Abstract 



Effective theories of quantum liquids (superconductors and superfluids of various types) 
are derived starting from microscopic models at the absolute zero of temperature. Spe- 
cial care is taken to assure Galilei invariance. The effective theories are employed 
to investigate the quantum numbers carried by the topological defects present in the 
phases with spontaneously broken symmetries. Due to topological terms induced by 
quantum fluctuations, these numbers are sometimes found to be fractional. The zero- 
temperature effective theories are further used to study the quantum critical behavior 
of the liquid-to-insulator transition which these systems undergo as the applied mag- 
netic field, the amount of impurities, or the charge carrier density varies. The classical, 
finite-temperature phase transitions to the normal state are discussed from the point of 
view of dual theories, where the defects of the original formulation become the ele- 
mentary excitations. A connection with bosonization is pointed out. 
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Preface 



The title of this report. Boulevard of Broken Symmetries, is chosen to indicate that the 
(condensed matter) systems we will be discussing in these pages have a spontaneously 
broken symmetry in common^ The notion of a spontaneously broken symmetry is one 
of the paradigms of modern physics. It plays a central role in today's understanding 
of many startling phenomena known in condensed matter physics, statistical physics, 
high-energy physics, and cosmology. A spontaneously broken symmetry indicates the 
presence of a global symmetry which is not apparent in the state the system is in. The 
symmetry is not lost, but implemented in a nontrivial way. If the symmetry involved 
is a continuous one (and the dimensionality is larger than two), the breaking is ac- 
companied by the occurrence of gapless modes. Due to their gaplessness, they are 
the dominant elements in an effective description of the system with broken symmetry 
valid at low energy and small momentum. 

This will be the subject of the first chapter The topic has a long history, what is rel- 
atively new is how the effective description can be reconciled with Galilei invariance — 
one of the symmetries governing the nonrelativistic world of condensed matter physics. 
More specific, we will derive the effective theories at the absolute zero of temperature 
of classical hydrodynamics, of a superconductor both in the weak-coupling as well as 
in the strong-coupling limit, of superfluid '^He, and of a bosonic superfluid. We also 
consider the behavior of a superconductor and the bosonic superfluid at finite tempera- 
ture. Throughout this report we employ the powerful apparatus of quantum field theory, 
mostly in the functional-integral formulation. In this approach, the effective theories 
are obtained by integrating out certain degrees of freedom contained in the microscopic 
description of the quantum system under consideration. 

The occurrence of gapless modes is not the only manifestation of a broken contin- 
uous symmetry. A closely related one is the appearance of topological defects. These 
objects, which are often of paramount importance to understand the physical properties 
of a system with a broken symmetry, can have peculiar quantum numbers associated 
with them. This will be the subject of Chap. ^ We shall consider a superfluid -^He 
film, one-dimensional metals, a two-dimensional model exhibiting an exotic mecha- 
nism which leads to superconductivity, and free electrons confined to move in a plane. 
Again using the functional-integral approach to quantum field theory, we shall see how 

*The true story behind this title is the following. Every late summer during the eighties, Amsterdam was 
visited by a traveling theater group. They built up their tents in a central park of the city surrounded by the 
Concert Building and the major museums (and also by the high school this author attended). The name of 
this enterprise, marking the end of the summer and its fantasies, was called "Boulevard of Broken Dreams" 
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so-called topological terms arise in the effective theories after the fermionic degrees of 
freedom contained in the microscopic models are integrated out. These terms are the 
cause of the peculiar quantum numbers which are expected to arise in these systems. 

The next chapter. Chap. ^ deals with so-called dual theories. Duality is one of 
the hottest topics in contemporary condensed matter physics. Quite generally, a dual 
description of a system refers to a description in terms of a different set of variables 
than the original set. The systems we will consider in this chapter — a superfluid *He 
film, a superconducting film, and a bulk superconductor, all at finite temperature — 
have vortex solutions in common. In the original formulation, these topological defects 
are described as singular objects. In the dual theories discussed in this chapter, these 
objects are instead described in a nonsingular fashion by field theory. The dual theories 
are employed to discuss the (equilibrium) phase transition these systems undergo as a 
function of temperature. It will also be pointed out in this chapter that the concept of 
duality in two dimensions is closely related to bosonization — a powerful computational 
tool in this reduced dimensionality. 

In the last chapter. Chap. Q we discuss quantum phase transitions. These are phase 
transitions, taking place close to the absolute zero of temperature, which are driven not 
by temperature, but by some other parameter in the system. In contrast to equilibrium 
transitions at finite temperature, time is important in quantum phase transitions. This 
naturally leads to the use of quantum field theory to describe these. Recent experi- 
ments on quantum phase transitions in various two-dimensional systems have raised 
some very interesting questions. It is generally believed that the concept of duality 
will play a decisive role in answering these. We will be concerned in this chapter with 
2nd-order quantum phase transitions. In addition to a diverging correlation length, 
which is common for 2nd-oder equilibrium transitions at finite temperature, these tran- 
sitions also have a diverging correlation time. We will discuss in detail the universality 
class defined by repulsively interacting bosons which undergo a so-called superfluid- 
to-Mott-insulating phase transition. We will consider the system both in the absence 
and in the presence of disorder. We then continue to describe the quantum phase tran- 
sition which a quantum-Hall liquid undergoes in the absence of disorder as the applied 
magnetic field varies. We shall discuss scaling theory and consider some relevant ex- 
periments. 
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Notation 



We adopt Feynman's notation and denote a spacetime point hy x = = {t, x), ^ = 
0, 1 , • • • , d, with d the number of space dimensions, while the energy fcg and momentum 
k of a particle will be denoted by fc = A:^ = {ko, k). The time derivative do = d/dt 
and the gradient V are sometimes combined in a single vector 9^ = (9o, — V). The 
tilde on is to alert the reader for the minus sign appearing in the spatial components 
of this vector. We define the scalar product k ■ x = k^x^ = k^Quvkv = k^t — k • x, 
with Qfj^v = diag(l, —1, • • • , —1) and use Einstein's summation convention. Because 
of the minus sign in the definition of the vector it follows that ^^^a^n = doao + V • a, 
with an arbitrary vector. 

Integrals over spacetime are denoted by 



When no integration limits are indicated, the integrals are assumed to be over all pos- 
sible values of the integration variables. 

Natural units ?i = c = fee = 1 are adopted throughout unless expUcitly stated. 




while those over energy and momentum by 
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Chapter 1 

Nonrelativistic Effective 
Theories 



In this chapter we shall derive effective theories governing the low-energy, small- 
momentum behavior of some nonrelativistic systems. The systems considered here 
have in common a spontaneously broken global symmetry. On account of Goldstone's 
theorem such a breakdown is always accompanied by gapless modes. If the origi- 
nal symmetry G is spontaneously broken to some subgroup H, the Goldstone modes 
parameterize the coset space G/H. Being gapless, these modes are the dominant exci- 
tations at low energy and small momentum and for this reason the relevant degrees of 
freedom to build effective theories from. 

In a Lorenz-invariant theory, massive elementary excitations have a spectrum of the 
form £'^(k) ^ h^+m^, with k the momentum and m the mass of the excitation. In this 
case, a gapless excitation is easily described by taking the limit m ^ in the massive 
theory. In a Galilei-invariant theory, however, it is a priori not clear how to describe 
a gapless mode. Here, elementary excitations of a free theory have a spectrum of the 
form i?(k) ~ k^/2m, so that the limit m ^ cannot be taken. Although elements 
of the solution to this problem were present in the literature (see, e.g., [|l} g]), it was 
not until a paper by Greitner, Wilczek, and Witten that the issue was settled. They 
considered the nonrelativistic Goldstone mode of a spontaneously broken global U(l) 
symmetry generated by the total particle number Their general considerations based 
on symmetry principles were powerful enough to fully determine the effective theory. 

In the following, we shall compute the effective theory describing the nonrela- 



tivistic Goldstone modes as they arise in an ideal classical fluid (Sec. 1.1), in a BCS 
superconductor (Sec. 1.3), in superfluid ^^He-a (Sec. 1.5), and in a weakly interacting 



Rose gas (S ec. |l.6| ). In addition, we investigate the strong-coupling limit of the pairing 



theory (Sec. 1 .4 ) and the high-temperature limit of a weakly interacting Rose gas (Sec. 



1.7) as well as that of the RCS theory (Sec. 1.8). 
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1.1 Effective Theory of Hydrodynamics 

The hydrodynamics of an ideal, classical fluid was already well understood in the 19th 
century. The case of isentropic flow, for which the entropy per unit mass is constant, 
is particularly simple. The pressure P is then a function of the mass density p only, 
and the flow is automatically a potential flow. A feature of such a fluid is that it sup- 
ports unattenuated sound waves, i.e., propagating density oscillations. The waves are 
unattenuated because viscosity and thermal conductivity, which usually serve to dis- 
sipate the energy of a propagating mode, are absent. More important to our present 
considerations is that sound waves are gapless. It would be gratifying to identify them 
as the Goldstone mode of a broken continuous symmetry since that would explain their 
gaplessness. We will argue in this section that such an identification is indeed possible. 

To describe the hydrodynamics of an isentropic fluid, we use Eckart's variational 
principle [Q] and start with the Lagrangian 

ipv2-pe + ,/)[aop + V- (pv)], (1.1) 

where v is the velocity field, p the mass density, and e the internal energy per unit 



mass. For isentropic flow e is a function of p alone. The first and second term in (1.1) 
represent the kinetic and potential energy density, respectively. The variable is a 
Lagrange multiplier introduced to impose the conservation of mass: 

<9op + V-(pv)=0; (1.2) 



its dimension is [(p] = m s . The variation of ( |1.1[ ) with respect to v yields the 
equation 

v = V(/). (1.3) 

It shows that, indeed, isentropic flow is automatically a potential flow and it also iden- 
tifies the Lagrange multiplier (j) as the velocity potential. (How to incorporate vortices 



in this variational approach will be discussed in the next section.) With ( |1.3[ ), the La- 
grangian ( |l.l[ ) becomes 

/: = -p[ao0 + i(v0)2 + e], (1.4) 

where we performed integrations by part. A second field equation can be obtained by 
varying the Lagrangian with respect to p. This yields the Bernoulli equation 

dQ(l) + ^{W4>f + h^0, (1.5) 

with h — d{pe)/dp the specific enthalpy. From this definition of h, one can easily 
derive the thermodynamic relation for isentropic flow 

\/h=-\/P, (1.6) 
P 



with P — p^de/dp the pressure. On taking the gradient of (1.5) and using (1.6) one 
obtains Euler's equation 

aov + iVv^ + -VP = (1.7) 

p 



(glAmstex 



1,1 Effective Theory of Hydrodynamics 



3 



governing the flow of the fluid. 

We next wish to investigate the symmetry content of the theory. Classical hydrody- 
namics has the following invariances: 



(i) Invariance under spacetime translations, 



C/j +e^i, with a constant vector. 



(ii) Invariance under global translations of the velocity potential, 
a a constant. 



+ a, with 



(iii) Invariance under Galilei boosts. 



t - 
do 



t' = t, 



d'o = do 



X 

u-V, V 



x' = X — ut; 
> V = V, 



(1.8) 



with u a constant velocity. 



According to Noether's theorem, symmetries imply conservation laws. Associated with 
the above invariances we obtain the following conservation laws: 

(i) dp,t^,y = 0, where t^i, i' — 0,1, 2, 3) is the energy-momentum tensor, with 
toi = Pi the momentum density and <oo = Ti- the Hamiltonian; 

(ii) 9^5^ = 0, with go the mass density and gi the mass current; 

(iii) dfigf_ij — 0, with goj = —goxj + tpj and gij — giXj — ttij the corresponding 
charge densities and currents. Physically, the conservation dGoi/dt = of the 
charges Go; — goi means that the center of mass of the fluid, X = xp/M, 
with M = J^p the total mass, moves with constant velocity. 



M— = 
dt 



Here, the right-hand side denotes the total momentum of the fluid. 



(1.9) 



From the Lagrangian (1.4) we obtain as explicit form for the various charge densities 
and currents [|]|: 



90 

Pji^ toj) 
tij 

W(= too) 



dC 

ddocj) 
dC 

m4> 

dC 



P 



ddocj) 



£5. 



djcf) = pvj 
dC 



ddoq 
dC 



djcj) = PSij + pvtVj 
do4' — ^ = + pe 



dd, 

dC 



do4> = {n + p)v. 



(1.10) 

(1.11) 

(1.12) 
(1.13) 
(1.14) 
(1.15) 
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Time derivatives da<j) have been eliminated through the field equation (1.5), so that, for 
example, C in the last equation is replaced with 



C ^ ph ~ pe 



d_ 

dp 



1 (pe) = P. 



(1.16) 



A few remarks are in order First, the Hamiltonian H is the sum of the kinetic and 
potential energy density, as required. Second, the equivalence of the mass current g 
and the momentum density p, which is a hallmark of Galilei invariance, is satisfied by 
the theory. Finally, the set of equations (1.1C)-(1.14) constitutes all the equations of 
hydrodynamics. This brings us to the conclusion that the Lagrangian (1.4) encodes all 
the relevant information for the description of an isentropic fluid. 

We next turn to the description of sound waves. We restrict ourselves to waves 
of small amplitude. These generate only small deviations in the mass density p and 
pressure P of the fluid at rest, so that we can expand the Lagrangian ( 1 .4 ) in powers of 
p~ p — p, with IpI << p: 



idocb + W)iP + p)-ep-hp~ Yh'-p" + O(p^) 



(1.17) 



The derivative (') is with respect to p and i s to be evaluated at p — p. Since for the 
system at rest (j> is constant, it follows from (1.5) that h ~ Q. If we denote the thermo- 
dynamic derivative dP/dp by c^, which has the dimension of a velocity squared, the 
coefficient of the quadratic term in p can be written as 



1 



(1.18) 



p p 

Apart from an irrelevant constant term (— ep) the Lagrangian becomes to this order 



We next eliminate p by substituting 



(1.19) 



P - 



(1.20) 



which follows from expanding the field equation (1.5). Physically, this equation reflects 
BernoulU's principle: in regions of rapid flow, the mass density p — p + p and therefore 
the pressure is low. It also shows that the expansion in p is one in derivatives. The 
higher-order terms that have been neglected correspond to higher-order derivatives. 
At low energy and small momentum, these additional terms can be ignored. After 
eliminating p, we obtain a Lagrangian governing the velocity potential [S]: 



2 ^ 



(1.21) 



This is precisely the effective theory of a nonrelativistic (Abelian) Goldstone mode 
obtained in Ref. [|]] using general arguments. Apart from an irrelevant constant it is 
identical to a proposal by Takahashi [gl which was also based on symmetry principles. 
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The field equation one obtains for the velocity potential from ( |1.21 1 is nonlinear: 



(1.22) 



The information contained in this equation cannot be more than the conservation of 
mass because was initially introduced in ( 1 . 1 ) as a Lagrange multiplier precisely to 
enforce this conservation law. Indeed, remembering that the combination p{da(f)+ jv^) 
denotes c§ times = p{x) — p, w ith p the constant mass density of t he flu id at 
rest, we see that ( 1.22 ) reproduces (L2) in this approximation. To simplify ( 1.22 ), we 
replace j5 in the first and last term with the full mass density p (which is justified to this 
order) to arrive at the known but unfamiliar field equation 



(1.23) 



of sound waves. If we ignore the nonlinear terms, it becomes the more familiar wave 
equation 

9o'0-cgv20 = O, (1.24) 

implying a gapless linear spectrum, and identifying c, which was introduced via the 
thermodynamic derivative dP/dp = c^, as the sound velocity. 

The combination do(j> + ^iVcj))^ appearing in the description of a nonrelativistic 
gapless field is dictated by Galilei invariance. To obtain the transformation property 
of the velocity potential under a Galilei boost (1.8) we note that since Vcf) is a ve- 
locity field, y(j){x) V'0'(a;') — V(/)(a;) — u, with x' = {t,x — ut). This gives as 
transformation rule for 6 



4>{x) ^ 4>'{x') ^ c^{x) ^xx-^ + f{t), 



(1.25) 



with f{t) a yet undetermined function of time. To determine fit) we note that the 
factor —dQcj) in the Lagrangian (1.4) is the chemical potential per unit mass. Indeed, 



using the standard definition p, ~ dTi / dp, we find 

p = ^v^ + h = -do(j), 



(1.26) 



where in the second equality we used the field equation (1.5). This identification fixes 
the transformation rule of —do<j>: 



doH^) -d'o(f>'{x') = -do^{x) - u • v{x) + y 



(1.27) 



and in combination with (ITU), yields for f{t) 



(1.28) 



ao/(i) = iu2, or /W-iu^i 

up to an irrelevant constant. It is easily checked that both the combination do4> + 
i(V(/))^ appearing in the effective theory (1.21 ) as well as the field equation (1.23) are 



invariant under Galilei boosts. So, contrary to what is sometimes stated in the literature 
[ 8[], s ound waves are invariant under Galilei boosts. The linearized wave equation 
( 1 .24 ) is, of course, not invariant because essential nonlinear terms are omitted. 
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Greitner, Wilczek, and Witten arrived at the effective Lagrangian (1.21) by requir- 



ing that the effective theory of sound waves can only be constructed with the help of 
the Galilei-invariant combination do<j>+ ^{Wcj))^, and that it should give the dispersion 
relation i?2(k) 



From the effective Lagrangian (1.21) one can again calculate the various Noether 
charge densities and currents 1^]. They are, as might be expected, of the same form as 



the exact expressions (1.10)-(1.14),but now with the approximations 



(1.29) 



as follows from Eq. ( |1.20[ ), 



n 



2 



2p 



and 



P 



~p{do(t)+W) =i clip- P)- 



(1.30) 



(1.31) 



This last equation is consistent with the expression one obtains from directly expanding 
the pressure: P{p) = P + pP' since P = and P' = Cq. 

To recapitulate, we have derived an effective theory describing a gapless mode 
starting from the Lagrangian (1.4) which entails the complete hydrodynamics of an 
isentropic fluid. We now wish to argue that this gapless mode is a Goldstone mode 
associated with a spontaneously broken symmetry. A first indication in favor of such 
an interpretation follows because the Hamiltonian H displays a property typical for a 
system with broken symmetry, namely that it is a function not of the velocity potential 
itself, but of the gradient of the field. The energy is minimal if is uniform throughout 
the sample, i.e., there is rigidity 

Usually the broken symmetry can be identified by the general property that the 
Goldstone mode is translated under the broken symmetry operations. (There may be 
other effects too, but the translation is always present.) Here, this general characteristic 
does not uniquely identify the broken symmetry because is translated under two 



symmetry operations. According to the transformation rule ( 1.25 ) with f{t) given in 
(1.28), we have that under a Galilei boost 



6^(j){x) := (t)'{x) - (t){x) = -u • X + tu • V0(a;), 



(1.32) 



where we took the transformation parameter u infinitesimal small so that quadratic and 
higher powers in u may be ignored. The first term at the right-hand side shows that 
the velocity potential is translated under a Galilei boost. The second symmetry under 
which (j) is translated is generated by the total particle number, or, equivalently, by the 
total mass M = p. To see this we first compute from (1.4) the conjugate momentum 
7r0 of (f), 

implying that and p are canonically conjugate [|lO[]: 



Wt,x),p(t,x')} = -<5(x-x'), 



(1.34) 
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where { , } denotes the Poisson bracket. We then use a central result of classical field 
theory stating that the charge Q of a continuous symmetry is the generator of the cor- 
responding transformations of the fields, xi^) ^ x'i^)- More specifically, for an 
infinitesimal transformation 5^x{x) — x'i^) ^ x{x) one has 

S^xix)=-a{x{x),Q}, (1.35) 



with a the transformation parameter Equation ( 1 .34 ) thus implies that under the global 



U(l) symmetry generated by the total mass, </> is indeed translated 

5^'^{x) = -a{^{x), M} = a. (1.36) 

The point is that the generators of Galilei boosts Goj — J^{—Xjgo + tpj) also contain 
the mass density p = go- It is therefore impossible to distinguish a broken Galilei 
invariance from a broken mass symmetry by considering the algebra alone. 

Let us at this point pause for a moment and consider the case of a superfluid. It is 
well established that in the normal-to-superfluid phase transition, the global U(l) sym- 
metry generated by the total mass is spontaneously broken. This comes about because 
in a superfluid, which is a quantum system, many particles Bose-Einstein condense in 
a single quantum state. The coherence allows us to describe the system by a complex 
field 

4>{x) = Vp(s)/me""'^(")/'' (1.37) 

normalized such that yields the particle number density p/m of the condensate. 
To underscore the quantum nature of a superfluid, Planck's constant has been made 
explicit. The field ip, which turns out to describe the Goldstone mode of the broken 
global U(l) symmetry, is a phase field and therefore compact. The field tp is frequently 
referred to as the condensate wavefunction, even in the modern literature. In our view, 
this is somewhat misleading. As has been stressed by Feynman [pi]], xjj is a classical 
field describing the coherent behavior of many condensed particles in the same way 
as the classical field of electrodynamics describes the behavior of many photons in a 
single state. For these classical fields there is no probability interpretation as is required 
for wavefunctions JI^. 

It has been argued by Feynman jrT[ ] that the i/'-field is governed by a nonrelativistic 
I ■01 ^-theory defined by the Lagrangian 

C^, = ihrdo^ - ^^IVV-P - ^(m|0p - p)\ (1.38) 
2m 2p 

The potential energy has its minimum along a circle away from the origin at = 
p/m, implying a spontaneous breakdown of the global U(l) symmetry. The Lagrangian 
(1.38) with (1.37) reduces to the one given in ( 1.19| ) when the term — fi.^(Vp)^/8m^p 



is ignored. Using the expression for the pressure, cf. (1.16) 
P 



dp ^ ddip 



(pe), (1.39) 



we find that it gives the contribution — ?i^(V^p)/4m^ to the pressure — the so-called 
quantum pressure. The reason for calling it this way is that it is the only place where 
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Planck's constant appears in the equations. To the order in which we are working, it 



is consistent to ignore this term. Because the field equation for V' derived from (1.38) 
has the form of a nonlinear Schrodinger equation, we will refer to ^/i as a Schrodinger 
field. We trust however that the reader realizes that it is a classical field unrelated to a 
Schrodinger wavefunction. 

Let us compare the transformation properties of the Schrodinger field with that of 
the velocity potential cj) of an isentropic fluid. Under a Galilei boost, ip{x) transforms 
as[|l| 

Tp{x) ip'{x') = exp[i{-u ■ X + ^uh)m/h] iP{x). (1.40) 

With rrvp/fi denoting the phase of the Schrodinger field, we see that ip transforms in 
the same way as does the velocity potential. 

Using that the canonical conjugate of the T/i-field is tt^ = ihi/;*, we easily derive 
the Poisson bracket 

{ipix), M} = -i{m/h)^{x). (1.41) 

This shows that the total mass M generates phase transformations on the T/i-field: 
'(p{x) tp'{x) = exp{ia'm/'h)ip{x). The phase ip of the Schrodinger field is conse- 
quently translated under the symmetry, just like the velocity potential cj). This transfor- 
mation property identifies Lp as the Goldstone mode of the broken U(l) mass symmetry. 



(1.34) 



The Poisson bracket (1.41) also implies that (p and p are canonical conjugate M\, cf. 



Mt,x),p(t,x')} = -5(x-x'). (1.42) 

A similar relation holds for superconductors. On quantizing, the Poisson bracket is 
replaced by a commutator. The Heisenberg uncertainty relation that results for the 
conjugate pair has recently been demonstrated experimentally [jTsj]. 

As remarked above, a necessary condition for the spontaneous breakdown of the 
global U(l) symmetry is the presence of a condensate. Such an intrinsic quantum 
phenomenon, requiring many particles in a single state, has no analog in a classical 
setting. Hence, the U(l) symmetry cannot be broken in classical hydrodynamics. This 
leaves us with the second possibility, namely that of a spontaneously broken Galilei 
invariance. The breakdown is a result of the presence of a finite mass density. This can 
be inferred [Q] from considering the transformation of the momentum density p(a;) 
under a Galilei boost: 

S^Pi{x) = ~Uj{pi{x), Goj{t)} = -u.ip{x) + tu ■ \7p,{x), (1.43) 
or with 6^ = "j'^F 

Sfpiix) = -p{x)Sj, + tdjPiix). (1.44) 

If the mass density p is finite, the right-hand side is nonzero, which is a symmetry- 
breaking condition. 



1.2 Including Vortices 

There is an essential difference between the spontaneous breakdown of the Galilei sym- 
metry and that of the global U(l) symmetry. Although both symmetries are Abelian, 
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the latter is a compact symmetry, whereas the Galilei group is noncompact. More 
specifically, the transformation parameter a of the U(l) group has a finite domain 
(0 < a < 27r), while the domain of u, the transformation parameter of the Galilei 
group, is infinite. As a result, the velocity potential of classical hydrodynamics cannot 
be represented as the phase of a complex field. An immediate physical manifestation 
of the difference is that a system with broken U(l) invariance supports topologically 
stable vortices, whereas a system with broken Galilei invariance does not. This is not 
to say that vortices are absent in the latter case, it merely states that their stability is 
not guaranteed by topological conservation laws. Closely connected to this is that the 
circulation is not quantized in classical hydrodynamics, which is known to exist in su- 
perfluids. Yet, the circulation is conserved also in isentropic fluids. This is again not 
for topological, but for dynamical reasons, the conservation being proven by invoking 



Euler's equation ( 1 .7 ) as was first done for an ideal, incompressible fluid by Helmholtz 
[ ]l4] ] and generalized to a compressible fluid by Thomson [|l5|]. 

The easiest way to observe vortices in a classical fluid is to punch a hole in the 
bottom of the vessel containing the fluid. As the fluid pours out a vortex is formed in 
the remaining fluid — a phenomenon daily observed by people unplugging a sinkhole. 
Often, as happens in, for example, superfluid '*He, the core of a vortex is in the normal 
state so that the symmetry is restored there. In the present context this would mean that 
inside the vortex core, the fluid mass density p is zero. This is indeed what is observed: 
the vortex core consists of air, therefore no fluid is present and p = there. 

In the eye of a tropical cyclone — another example of a vortex, nature does its best to 
restore the Galilei symmetry, record low atmospheric pressures being measured there. 
(A complete restoration would imply the absence of air corresponding to zero pres- 
sure.) 

It is customary to incorporate vortices in a potential flow via the introduction of 
so-called Clebsch potentials [|T^. We will not follow this route, but instead use the 
powerful principle of defect gauge symmetry developed by Kleinert [|l7l [T^ll^. In 
this approach, one introduces a so-called vortex gauge field 0^ = (0^ i ) in the 
Lagrangian via minimally coupling to the Goldstone field; 

~^^,,l>^~^^,<i> + ^l, (1.45) 

with — (9o, — V) and 

V X 0P = -2a;, (1.46) 

so that V X V = 2a; yields (twice) the vorticity a; of the vortex. The combination 
9^0 + is invariant under the local gauge transformation 

<i>(x) ^ + a(x); d^aix), (1.47) 

with 0^ playing the role of a gauge field. The left-hand side of ( 1.46| ) may be thought of 



as defining the "magnetic field" associated with the vortex gauge field = V x ^ 

For illustrative purposes, let us consider an external, static vortex with circulation F 
located along a line L, which may be closed or infinitely long [|l^ . Then,a; = ^rd{L), 
where S {L) is a delta function on the line L, 

5,{L)^ J^dy,S{^-y). (1.48) 
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This model with a static, external vortex may be thought of as describing the steady 
flow in the presence of a vortex pinned to a fixed impurity. The field equation for <j> 



obtained after the substitution (1.45) reads 



dofi + cgV • V = dov^ + iv • Vv^ - V (1.49) 

with fi = —{dQ(f) + 4>q) the chemical potential and v = Vi/) — <p^ the velocity of the 
flow in the presence of the vortex. The last term gives a coupling of the velocity field 
to the "electric field" associated with (j)^, 

eP = - ao^P. (1.50) 



Note that the field equation ( 1.49 ) is invariant under local vortex gauge transformations. 
Ignoring the higher-order terms and choosing the gauge ipQ — 0, we obtain as equation 
for the flow in the presence of a static vortex: 

V-v = 0, or V • (V0 - ^'') = 0, (1.51) 

which is solved by 

0(x) = - /G(x-y)V-(/)P(y). (1.52) 
Here, G(x) is the Green function of the Laplace operator 

Straightforward manipulations then yield the well-known Biot-Savart law for the ve- 
locity field in the presence of a static vortex [|l4, 17 1 



r 

47r Jl |x - y|' 



vv(x) = — / dy X ' (1.54) 



where the integration is along the vortex. This exemplifies the viability of the vortex 
gauge principle as an alternative to describe vortices in a potential flow. 

Let us continue to study the dynamics of vortices — a subject that has recently re- 
ceived considerable attention in the literature [p0|]. Because the vortex motion is deter- 
mined by the flow itself, the vortex can no longer be considered as external. We shall 



see that the nonlinear part of the field equation (1.49) becomes relevant here. 

In the absence of external forces, the vortex moves with a constant velocity, say. 
The fl ow in the presence of a moving vortex can be obtained from the static solution 



( |1.54| ) by replacing the coordinate x with x — v^t. This implies that 

9oVv(x - WLt) = -VL • Vvv(x - VLt). (1.55) 

Since the solution Vv is curl-free outside the vortex core, the right-hand side may there 
be written as — V(vi • Vv). To study sound waves in the presence of a moving vortex, 
we write the velocity field as v(x) = Vv(x — v^t) + V(/>(x), with (f> describing small 
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variations around the moving vortex solution. Equation (1.55) then requires that we 



write for the chemical potential in (1.49) 



/i(a;) = vl • Vv(x - vl^) - 9o0(a;). (1.56) 
This leads to the linearized field equation [ pi} ] 

d^4>{x) - clV^ix) = -vv(x) • V9o[2^(t,x) - ^it,0)], (1.57) 



describing sound waves in the presence of a moving vortex. In deriving ( 1.57 ) we again 
used the gauge (pQ ~ 0, and approximated Vv(x — v^t) by Vv(x). To linear order in 
this is allowed since the vortex, being driven by the sound wave, has a velocity 

Vi(t) = W^{t, x = VLf) « V^(i, 0). (1.58) 

We also neglected a term quadratic in Vv which is justified because this velocity is 
much smaller than the sound velocity outside the vortex core [|T]] . The first term at 



the right-hand side in (1.57) stems from the nonlinear term cJqv in the general field 



equat ion ( |1.49| ). Thus the nonlinearity of sound waves becomes detectable. Equation 



(1.57) can be used as a basis to study the scattering of phonons of a free moving vortex 



2_ 

So far we have contrasted the spontaneous breakdown of the Galilei invariance 
(caused by a finite mass density) and that of the global U(l) symmetry (caused by 
a nonzero condensate). A superfluid, however, has a finite mass density as well as a 
nonzero condensate. Both symmetries are therefore broken and we expect two different 
Goldstone modes to be present. This is indeed what is observed in superfluid "^He. The 
system supports besides first sound, which are the usual density waves associated with 
the spontaneously broken Galilei invariance, also second sound, or entropy waves. The 
latter mode depends crucially on the presence of the condensate and is the Goldstone 
mode associated with the spontaneously broken global U(l) symmetry generated by 
the total mass. At the transition point, the second sound velocity vanishes, whereas 
the first sound velocity remains finite. This is as expected since only the condensate 
vanishes at the superfluid phase transition; the total mass density remains finite. 



1.3 Anderson-Bogoliubov Mode 

In this section we study the so-called Anderson-Bogoliubov mode ^ of a neutral 
superconductor. This mode is known to be the gapless Goldstone mode associated with 
the spontaneously broken global U(l) symmetry generated by the total mass. Given the 
general arguments of Ref. [|]| we expect that the effective theory of a superconductor 
is exactly of the form ( 1.2l[ ) we obtained for sound waves in classical hydrodynamics. 



but with replaced by a compact field. 

Our starting point is the famous microscopic model of Bardeen, Cooper, and Schri- 
effer (BCS) defined by the Lagrangian [|^ 

c = r^[ido - e(-«v)]VT + rd^do - a-i^Mi - Ao^* 

:= Ca + Cu (1.59) 
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where ll\ = — AqV'* ■0| i'x V'T ^ contact interaction term, representing the effective, 
phonon mediated, attraction between electrons with coupling constant Aq < 0, and Lq 
is the remainder. In ( 1.59| ), the field V't(I) anticommuting field describing the 
electrons with mass m and spin up (down); ^(— iV) = e(— iV) — /iq, with e(— iV) = 
—V'^ /2m, is the kinetic energy operator with the chemical potential /iq subtracted. 



The Lagrangian (1.59) is invariant under global U(l) transformations. Under such 



a transformation, the electron fields pick up an additional phase factor 



(1.60) 



with a =1, [ and a a constant. Notwithstanding its simple form, the microscopic 
model ( |L59| ) is a good starting point to describe BCS superconductors. The reason 
is that the interaction term allows for the formation of Cooper pairs which below a 
critical temperature condense. This results in a nonzero expectation value of the field 
A describing the Cooper pairs, and a spontaneous breakdown of the global U(l) sym- 
metry. This in turn gives rise to the gapless Anderson-Bogoliubov mode which — after 
incorporating the electromagnetic field — lies at the root of most startling properties of 
superconductors [p5[|. 

To obtain the effective theory of this mode, the fermionic degrees of freedom have 
to be integrated out. To this end we introduce Nambu's notation and rewrite the La- 
grangian ( 1 .59 ) in terms of a two-component field 



In this notation, Cq becomes 



V'T 



i9o - C(-«V) 






ido + i{-iV) 



(1.61) 



(1.62) 



where we explicitly employed the anticommuting character of the electron fields and 
neglected terms which are a total derivative. The partition function. 



Z 



exp 



(1.63) 



must for our purpose be written in a form bilinear in the electron fields. This is achieved 
by rewriting the quartic interaction term as a functional integral over auxiliary fields A 
and A*: 



(1.64) 



exp ^-lAo j -01 V'T ^[ V'T 



DA*DAgxp 



where, as always, an overall normalization factor is omitted. Classically, A merely 
abbreviates the product of two electron fields 



A = AoV'iV'T- 



(1.65) 
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It would therefore be more appropriate to give A two spin labels A|| . Since and ipi 
are anticommuting fields, A is antisymmetric in these indices. Physically, it describes 
the Cooper pairs of the superconducting state. 

By employing ( 1 .64 ), we can cast the partition function in the desired bihnear form: 

(1.66) 



Z = J DV'^DV' J DA*DA exp (^y J A* A 



X exp 



-A 



Changing the order of integration and performing the Gaussian integral over the Grass- 
mann fields, we obtain 



Z 



DA*DA exp iSes[A*,A] + 



An 



A*A , 



(1.67) 



where S'ofr is the one-loop effective action which, using the identity Det(A) = exp[Tr 
ln(A)], can be cast in the form 



5eff[A*,A] = -zTrln 



Po - e(p) -A 

-A* po + e(p) 



(1.68) 



where po = ido and ^(p) e(p) — /iQ, with e(p) = p^/2m. The trace Tr appearing 
here needs some explanation. Explicitly, it is defined as 



S'ofr = -iTr \n[K{p,x)] = -itrln K{p, x)S{x - y)\^^^ 



(1.69) 

where the trace tr is the usual one over discrete indices. We abbreviated the matrix 
appearing in (1.68) by K{p, x) so as to cover the entire class of actions of the form 



S 



{x)K{p, x)'ip{x) 



(1.70) 



The delta function in ( 1.69 ) arises because K{p, x) is obtained as a second functional 
derivative of the action 



6^S 



Stp^{x) 5ip(x) 



K{p, x) 5{x - y)\ 



y=x' 



(1.71) 



each of which gives a delta function. Since the action has only one integral over 
spacetime, one del ta fun ction remains. Because it is diagonal, it may be taken out of 
the logarithm and ( 1.69 ) can be written as 



'eff 



-itv I lii[K{p,x)]S{x-y)\y^^ 



-itr 



ik-x 



e'*-^ In [if (p, x)] 6-^*=- 



(1.72) 
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In the last step, we used the integral representation of the delta function: 

S{x) = [ e-''^-^, (1.73) 



k 



shifted the exponential function exp(iA; • y) to the left, which is justified because the 
derivative s^ doe s not operate on it, and, finally, set equal to x^. We thus see that the 



trace Tr in ( 1.69 ) stands for the trace over discrete indices as well as the integration over 
spacetime and over energy and momentum. The integral J^, arises because the effective 
action calculated here is a one-loop result with fc^ the loop energy and momentum. 



The integrals in (1.72) cannot in general be evaluated in closed form because the 
logarithm contains energy-momentum operators and spacetime-dependent functions in 
a mixed order To disentangle the integrals resort has to be taken to a derivative ex- 
pansion in which the logarithm is expanded in a Taylor series. Each term contains 
powers of the energy-momentum operator which acts on every spacetime-dependent 
function to its right. All these operators are shifted to the left by repeatedly applying 
the identity 

f{x)pf,g{x) [pf, - id^)fix)g{x), (1.74) 

where f{x) and g{x) are arbitrary functions of spacetime and the derivative 9^ = 
{do, —V) acts only on the next object to the right. One then integrates by parts, so 
that all the p^'s act to the left where only a factor exp(i/c • x) stands. Ignoring total 
derivatives and taking into account the minus signs that arise when integrating by parts, 
one sees that all occurrences of (an operator) are replaced with fc^ (an integration 
variable). The exponential function cxp{—ik ■ x) can at this stage be moved to the left 
where it is annihilated by the function exp(j/e • x). The energy-momentum integration 
can now in principle be carried out and the effective action be cast in the form of an 
integral over a local density CcS- 

Scs ^ f C,s. (1.75) 



This is in a nutshell how the derivative expansion works. 

In the mean-field approximation, the functional integral ( 1.67 ) is approximated by 
the saddle point: 

Z = exp (^tS,s[A:^^f, A^f] + ^J^ , (1.76) 

where Amf is the solution of mean-field equation 

= --A. (1.77) 



If we assume the system to be spacetime independent so that A,nf (a;) = A, Eq. ( 1.77 ) 
yields the celebrated BCS gap equation: 

1 f 1 



1 f 1 



2 7k E{\^y 



(1.78) 
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where rj is an infinitesimal positive constant that is to be set to zero at the end of the 
calculation, and 



E{k) = ^eO^) + |AP (1.79) 
is the spectrum of the elementary fermionic excitations. If this equation yields a solu- 



tion with A 7^ 0, the global U(l) symmetry (1.60) is spontaneously broken since 



A->e2'"A^A (1.80) 

under this transformation. The factor 2 in the exponential function arises because A, 
describing the Cooper pairs, is built from two electron fields. It satisfies Landau's 
definition of an order parameter as its value is zero in the symmetric, disordered state 
and nonzero in the state with broken symmetry. It directly measures whether the U(l) 
symmetry is spontaneously broken. 



In the case of a spacetime-independent system, the effective action ( 1 .68 ) is readily 
evaluated. Writing 

po - e(p) -A ^^fpo" ap) 

-A* PO + e(p) \ PO + ^(P) 



(1.81) 



A 
A* 

and expanding the second logarithm in a Taylor series, we recognize the form 

up to an irrelevant constant. The integral over the loop energy fcp gives for the corre- 
sponding effective Lagrangian 

£eff = ^ [i?(k) - C(k)] . (1.83) 

To this one-loop result we have to add the tree term |Ap/Ao. Expanding -E(k) in a 
Taylor series, we see that the effective Lagrangian also contains a term quadratic in A. 
This term amounts to a renormalization of the couphng constant; we find to this order 
for the renormalized coupling constant A: 

1 1 1 /■ 1 

(1.84) 



A Ao 2Ale(k)r 

The integral at the right-hand side diverges, to regularize it we introduce a momentum 
cutoff A. In this way we obtain 

A Aq zvr^ 
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where we omitted the (irrelevant) finite part of the integral. It should be remembered 
that the bare coupling constant Ao is negative, so that there is an attractive interaction 
between the fermions. We can distinguish two limits. One, if the bare coupling constant 
is taken to zero, Aq 0^, which is the famous weak-coupling BCS limit. Second, the 
limit where the bare coupling is taken to minus infinity Aq — oo. This is the strong- 
coupling limit, where the attractive interaction is such that the fermions form tightly 
bound pairs [|27|]. These composite bosons have a weak repulsive interaction and can 
undergo a Bose-Einstein condensation (see succeeding section). 

To this order there is no renormalization of the chemical potential, so that we can 
write /i = /iQ. 

Since there are two unknowns contained in the theory, A and /i, we need a second 
equation to determine these variables in the mean-field approximation. To find the 
second equation we note that the average fermion number N, which is obtained by 



differentiating the effective action (1.68) with respect to /i 



N^^, (1.86) 
is fixed. If the system is spacetime independent, this reduces to 

n = -itT f Go{k)T3, (1.87) 



k 



where n = N/V, with V the volume of the system, is the constant fermion number 
density, T3 is the diagonal Pauli matrix in Nambu space. 



and G'o(fc) is the Feynman propagator. 



Go(fc) = ( , ) (1-89) 



A* ko + e(k) 

1 / fco e*'"'^ + C(k) A 

Here, t] is an infinitesimal positive constant that is to be set to zero at the end of the 
calculation. The exponential functions in the diagonal elements of the propagator are 
an additional convergence factor needed in nonrelativistic theories [GSh. If the integral 



over the loop energy fco in the particle number equation (1.87) is carried out, it takes 
the familiar form 



The two equations (1.78) and (1.87) determine A and /i. They are usually evaluated in 
the weak-coupling BCS limit. However, as was first pointed out by Leggett [p7|], they 
can also be easily solved in the strong-coupling limit (see succeeding section), where 
the fermions are tightly bound in pairs. More recently, also the crossover between 
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the weak-coupling BCS limit and the strong-coupling composite boson limit has been 
studied in detail ||, |^ |l|, ||]. 

We continue by writing the order parameter Amt as 



(1.91) 



where A is a spacetime-independent solution of the mean-field equation (1.77). This 
approximation, where the phase of the order parameter is allowed to vary in spacetime 
while the modulus is kept fixed, is called the London limit. The phase field ip{x), which 
is a background field, physically represents the Goldstone mode of the spontaneously 
broken U(l) symmetry. Our objective is to derive the effective theory describing this 
Goldstone mode. 

To this end we decompose the Grassmann field as, cf. p2 



(1.92) 



and substitute the specific form (1.91 ) of the order parameter in the partition function 



(1.66). Instead of the effective action (1.68|) we now obtain 

-iTrln 



Sctf — 



Po - doip -_^(p + V^) -A 

-A* Po + do^ + e(p - Vy^) 



(1.93) 



where the derivative dfj^ip of the Goldstone field plays the role of an Abelian background 
gauge field. We next write this effective action in the equivalent form 

oo 

^eff = -*Trln [G^^ (1 - Go A)] - *Tr^ - (Go A)' , (1.94) 

e=i 

where in the last step we ignored an irrelevant constant and the matrix A is given by 



A(p) = Ut3 + Ip . V(p + -^VV, 
m Zm 



(1.95) 



with U the Galilei-invariant combination 



(1.96) 



We shall consider only the first two terms in the series at the right-hand side of ( 1 .94 ), 
and ignore for the moment higher than first derivatives on the Goldstone field (p. The 
first term yields 



On account of ( 1 .87 ), this can be written as 



(1.97) 



(1.98) 



where n — fcp/37r^ is the fermion number density. We see that the phase of the order 
parameter appears only in the Galilei-invariant combination ( 1 .96 ). 
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We continue with the second term in the series (1.94). We shall first restrict our- 



selves to the weak-coupling BCS limit, where the chemical potential is well approxi- 
mated by the Fermi energy, jj. w kp/2m. In this limit, we can make the approximation 



47r 



de « v{0) 



An 



(1.99) 



where i'{0) = mkp/2TT^ is the density of states per spin degree of freedom at the Fermi 
level and / d^fc denotes the integral over the solid angle. In the last step, we extended 
the range of the ^ = e — /i integration from [—^,oo) to (—00,00) which is justified by 
the rapid convergence of the integrals involved. After some algebra, using the integrals 



r(0 

£;2;+i ^ r(z + i)|A|2'^ 



1 



with / 



1,2,.. ., and T{1) the gamma function, we obtain 

1 



^(2) 



KO) 



n 2 



2m 



(V^)^ 



(1.100) 



(1.101) 



(1.102) 



which is again invariant under Galilei transformations. 



Equations (1.98) and (1.102) give the effective theory of the Goldstone mode asso- 



ciated with the spontaneous breakdown of the global U(l) symmetry [E 



^cff 



2m 



KO) 



2m 



(1.103) 



As expected, it is of the same form as the effective theory we obtained in Eq. (1.21 ) of 
the previous section describing a hydrodynamic sound wave. After rescaling (p mip, 
the two forms become identical. The only difference being that the present Goldstone 
field (pis a compact field. For the velocity c of the Anderson-Bogoliubov mode we find 



(1.104) 



where vp = kp/m is the Fermi velocity. The effective theory (1.103) has been red- 
erived by various other authors ^ 

Graphically, the effective theory is represented by the diagrams depicted in Fig. |1.1| , 
where a line with a shaded bubble inserted stands for i times the full Green function G 
and the black bubble denotes i times the full interaction F of the Xo--fields, introduced 
in ( 1.92 ), with the background field U which is denoted by a wiggly line. Both G and 
r are 2 X 2 matrices. The full interaction is obtained from the inverse Green function 
by differentiation with respect to the chemical potential. 



F = -- 



dG- 



(1.105) 
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Figure 1.1: Graphical representation of the effective theory (1.103). The symbols are 
explained in the text. 



This follows because U, as defined in (1.96), appears in the theory only in the combi- 
nation n — U. To lowest order, the interaction is given by —t^. Employing Eq. (1.87) 



with Go replaced by the full Green function, we conclude that the first diagram indeed 
corresponds to the first part of the effective theory ( 1.103| ). The second diagram with- 
out the wiggly lines denotes i times the (0 0)-component of the full polarization tensor, 
IIoo, at zero energy transfer and low momentum q, 

i limnoo(0,q) = limtr / rgGT G (fco, k + q), (1.106) 

q— q^O Jj, 

where the minus sign associated with the fermion loop is included. To see that this 
represents the second part of the effective theory we invoke an argument due to Gavoret 
and Nozieres [|8]| . By virtue of relation ( 1 . 1 05 ) between the full Green function G and 
the full interaction F, the (0 0)-component of the polarization tensor can be cast in the 
form 



lim noo(0,q) 

q-^O 



ilim tr / t^G — — G(fco,k + q) 
q^" Jk dfi 

. d 



lim tr / T3G{ko,k + q) 
Jk 



dn 
djjL 



Ik 

1 d'^n 



(1.107) 



where Q. is th e therm odynamic potential and V the volume of the system. The right- 
hand side of (1.107) is n^n, with k the compressibility. Because it is related to the 
macroscopic sound velocity c via 



(1.108) 



we conclude that the (0 0)-component of the full polarization tensor satisfies the so- 
called compressibility sum rule of statistical physics [BSD 



1 9^0 

lim noo(0,q) = -TJIT^ 



It now follows immediately that t he second diagram in Fig. 
part of the effective theory (1.103). 



1.1 



(1.109) 
represents the second 
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Both the particle number density and the sound velocity can thus also be obtained 
from the thermodynamic potential by differentiating it with respect to the chemical 
potential; 

idn 1 imd^n 
^^^T7ir'^ ~2=^T7~ir^- (1.110) 

In the weak-coupling BCS limit, the spectrum of the elementary fermionic excitations 
does not differ significantly from the free-particle spectrum. The thermodynamic po- 
tential may therefore be approximated by that of a free fermion gas. The partition 



function (1.63) with C the free Lagrangian ( 1 .62 ) reads 



Zfrco = exp y-i J Vfrocj , (1-1 H) 

where the spacetime integral over the potential Vh-cc is given by minus the effective 



action (1.68) with A set to zero. We used a potential rather than an action in (1.111) 
because the free fermion system is homogeneous in spacetime. Carrying out the inte- 
grals over the loop energy and momentum, we obtain for the potential 

Vfree = -^m^/V/'. (1.112) 

This physically represents the ground-state energy per unit volume. The thermody- 
namic potential (at the absolute zero of temperature) is obtained by integrating V over 
space, 

n = [v. (1.113) 



With the help of ( |l.llO| ), the velocity of the Anderson-Bogoliubov mode can now be 



calculated; this yields again (1.104) 



One can of course continue and calculate higher-order terms of the effective the- 
ory. Since these involve always higher-order derivatives, they are irrelevant at low 
energy and small momentum. To identify the expansion parameter, let us compute the 
quadratic terms in the Goldstone fields, involving fourth-order derivatives. A somewhat 
tedious but straightforward calculation yields in the weak-coupling BCS limit: 

= ^ [(^oV)^ - |«^(9oV)(VV) + i4(VV)^] • (1.114) 

This equation shows us that the higher-order terms have either an extra factor 9q / | 
or (wp V)^/| Ap, where we recall that the BCS correlation length at zero temperature 
is given by 

^o = ~- (1.115) 
7r|A| 

It sets the scale over which the modulus of the order parameter varies. At low energy 



and small momentum these additional terms can indeed be ignored. In deriving (1.114) 
we integrated various times by parts and neglected the ensuing total derivatives. With 
the additional terms included, the energy spectrum of the Anderson-Bogoliubov mode 
becomes 

E^(k) = i«2k2 (1 - ^^^Cok") , (1.116) 
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in accor dance with Refs. [|39[|^ |l|. The minus sign in ( |l.79[ ) shows the stability of 
this Goldstone mode against decaying. 

When one considers even higher energies and momenta, one also has to account 
for variations in the modulus of the order parameter so that the London limit, where 
the modulus is taken to be constant, is no longer applicable. 



1.4 Composite Boson Limit 

In this section we shall investigate the strong-coupling limit of the pairing theory. In 
this limit, the attractive interaction between the fermions is such that they form tightly 
bound pairs of mass 2m. To explicate this hmit in arbitrary space dimension d, we swap 
the bare coupling constant for a more convenient parameter — the binding energy of 
a fermion pair in vacuum [^. Both parameters characterize the strength of the contact 
interaction. To see the connection between the two, let us consider the Schrodinger 
equation for the problem at hand. In reduced coordinates, it reads 



V2 

+ Ao(5(x) 

m 



(1.117) 



where the reduced mass is m/2 and the delta-function potential, with Ao < 0, repre- 
sents the attractive contact interaction £[ in (1.59). We stress that this is a two-particle 



problem in vacuum; it is not the famous Cooper problem of two interacting fermions 
on top of a filled Fermi sea. The equation is most easily solved by Fourier transforming 
it. This yields the bound-state equation 



Ao 



or 



1 

a;^ 



1 



(1.118) 



(1.119) 



/k k2/m + ea' 

This equation allows us to swap the c oupling constant for the binding energy ea- When 
substituted in the gap equation ( 1.78 ), the latter becomes 



1 



1 



1 



E{k)- 



By inspection, it is easily seen that this equation has a solution given by 



0, 



Mo 



(1.120) 



(1.121) 



where it should be noted that the chemical potential is negative here. This is the strong- 
coupling limit. To appreciate the physical significance of the specific value found for 
the chemical potential in this limit, we note that the spectrum i?b (q) of the two-fermion 
bound state measured relative to the pair chemical potential 2/io reads 



4m 



(1.122) 
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The negative value for fiQ found in ( 1 . 1 2 1 ) is precisely the condition for a Bose-Einstein 
condensation of the composite bosons in the q = state. 



To investigate this limit further, we consider the effective action (1.68) and expand 

(1.123) 



A(a;) around a constant value A satisfying the gap equation (1.78), 

A(x) = A + A(a;). 
We obtain in this way, 



oc ^ 

1=1 



A 

A* 



(1.124) 



where Go is given in ( 1.89 ). We are interested in terms quadratic in A. Employing the 
derivative expansion outlined on the pages p^-[l4], we find 



1 



1 



2 pl-E^p)ipo + qo)^-E\p-<i) 

X {A^ a* a* + bo + e(p)]bo + go - e(p - q)] AA* 

.2 r r 



(1.125) 



+A*'aA + [po - e(p)]bo +qo+ e(p - q)]A*A}, 

where — id^. It is to be recalled here that the derivative operates on everything 
to its right, while operates only on the first object to its right. Let us for a moment 
ignore the derivatives in this expressi on. A fter carrying out the integral over the loop 
energy and using the gap equation ( 1.78 ), we then obtain 

^^'^0) = -^ / 747TfA'A*2 + A*'A2 + 2|AnA|2). (1.126) 
8 A -E^ (k) V J 

In the composite boson limit A ^ 0, so that the spectrum ( 1.79 ) of the elementary 
fermionic excitations can be approximated by 



i?(k) « 6(k) 



The remaining integrals in ( 1.126 ) become elementary in this limit, 

^ _ 4r(3-d/2) ,/2-3 



£;3(k) 



(1.127) 



(1.128) 



We next consider the terms involving derivatives in (1.125). Following Ref. 
we set A to zero here. The integral over the loop energy is easily carried out, with the 
result 



C^'\q) 



1 



1 



AA* 



2 7k —go — k^/m + 2^o — /Am 



A*A. 



(1.129) 
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The integral over the loop momentum k will be carried out using the dimensional- 
regularized integral 



1 



r(i-d/2) 1 



(k2 + A'P) (47r)'i/2 (Af2)i-''/2 



(1.130) 



to suppress irrelevant ultraviolet divergences. To illustrate the power of dimensional 
regularization, let us consider the case d = 3 in detail. Introducing a momentum 
cutoff, we find in the large-A limit 



1 



(k2 + A/2) 



27r2 47r 



(1.131) 



From (1.130) however only the finite part emerges. This exemplifies that terms diverg- 



ing with a strictly positive power of the momentum cutoff are suppressed in dimen- 
sional regularization. These contributions, which come from the ultraviolet region. 



cannot physically be very relevant because the simple BCS model (1.59) stops being 



valid here and new theories are required. It is a virtue of dimensional regularization 
that these irrelevant divergences are suppressed. 
We thus obtain in the strong-coupling limit 



1 



k go - k2/?7i - Ea - q2/4m 



(1.132) 



%-^m'^/2(-go + + q2/4m)^/2-\ 



or expanded in derivatives 
1 



k go - k2/?7i - - q2/4r 



(1.133) 



r(i - d/2) 



r(2-d/2) 



go 



q 

4m 



The first term at the right-hand side yields as contribution to the effective theory 



.(2) _ r(i 



d/2) 



{AttYP 



i'^/2e:?/2-i|A|2. 



(1.134) 



To this we have to add the contribution |A|2/Ao coming from the tree potential, i.e., 
the last term in the partition function (1.67). But this combination is no other than 
the one needed to define the renormalized coupling constant via (1.84), which in the 
strong-coupling limit reads using dimensional regularization 



1 
A 



1 ni~d/2) 



(1.135) 



In other words, the contribution ( 1.134| ) can be combined with the tree contribution to 
yield the term | A|2/A. Expanding the square root in (1.132) in powers of the derivative 
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Qfj, using the value (1.121) for the chemical potential, and pasting the pieces together, 
we obtain for the terms quadratic in A [p9|. 



(2) i r(2-d/2) rf/2-2 Mn(a) ^ ^ 

where Mo{q) is the 2x2 matrix. 



(1.136) 



/ go-qV4m-(2-d/2)|A|Ve, 
i -(2-d/2)A*'/ea 



(1.137) 



-(2-d/2)AVe, 
-(7o-qV4™-(2-d/2)|A|Ve„ 



As we will see shortly in Sec. |l.q , this describes the super fluid st ate of a weakly inter- 
acting composite boson system. On comparing with Eq. (1.170) below, we conclude 
that the composite bosons have — as expected — a mass rrib = 2m twice the fermion 
mass TO, and a small chemical potential 



/io,b = (2 - d/2) 



\A\' 



(1.138) 



From ( 1.137 ) one easily extracts the so-ca lled Bogoliubov spectrum and the velocity 
Co of the sound mode it describes (see Sec. 1.6), 



TOb mea 
Also the number density flo.b of condensed composite bosons, 

as well as the weak repulsive interaction Ao.b between the composite bosons, 

l-d/4 el-'/' 



Ao,b = (471)'*/= 



r(2-d/2) TO'*/^ 



(1.139) 



(1.140) 



(1.141) 



follow immediately. We in this way have explicitly demonstrated that the BCS theory 
in the composite boson limit maps onto the Bogoliubov theory. 

In concluding this section, we remark that in d = 2 various integrals we encoun- 
tered become elementary for arbitrary values of A. For example, the gap equation 
(1.120) reads explicitly in d = 2 



Ai2 + |A|2-^, 



while the particle number equation ( 1.9C ) becomes 



27r 



lAI 



(1.142) 



(1.143) 



(glAmstex 



1.5 Superfluid ^He-a 



25 



Since in two dimensions, 

fi^M^^ep, (1.144) 

with kp and ep — kp/2m the Fermi momentum and energy, the two equations can be 
combined to yield [ p4| 

ea v//i2 + |A|2 - n 

_a ^ 9 V ^_^L_L tl, (1.145) 

ep V/^2^|A|2+/i 

The composite boson limit we have been discussing in this section is easily retrieved 
from these more general equations. Also note that in this limit, h — 2no,b- 



1.5 Superfluid ^He-a 

In this section we extend the previous analysis to a more complicated system, namely 
that of superfluid "^He [Q. Whereas the Cooper-pairing in a BCS superconductor is in 
a spin-singlet state, the pairing in ''He is in a triplet state. This means that besides the 
global U(l)^ group generated by the particle number N also the 80(3)"^ spin rotation 
group has to be considered. Moreover, since on account of the exclusion principle 
the pairs must have an odd angular momentum, also the space rotation group must be 
included. In "'He, it is generally excepted that the pairing is in the L = 1 state. The 
relevant symmetry group is now much larger than that of a BCS system which leads 
to a much richer structure of possible superfluid phases [ p4] , p3| ]. Since in most phases 
the three symmetries become intertwined, superfluid ''He displays surprising physical 
properties [^. 

We will consider a '^He film. The specific superfluid state we will study is the two- 
dimensional analog of the three-dimensional '^He-A phase, the so-called '^He-a phase, 
first considered by Stein and Cross (see also Ref. [Q). This state is characterized 
by the magnetic quantum numbers ms — 0, niL = — 1, where the orbital magnetic 
quantum number may be considered as the projection of the orbital angular mo- 
mentum on the nonexisting z-axis. The corresponding symmetry breakdown is 

S0(3)^ X S0(2)^ X U(l)^ D S0(2)^ x U(l)^+^, (1.146) 

where the generator of the residual symmetry group U(l)^+^ is the sum of the gener- 
ators of the group SO(2)^ of space rotations and the group U(l)^ of phase transforma- 
tions. We see that in this state, the spin rotation group SO(3)'^ is spontaneously broken 
to SO(2)'^, as in ferro- and antiferromagnets. Since ms = 0, the superfluid ■^He-a 
phase is, in addition, an antiferromagnet. We choose to represent the order parameter 
by the matrix Af with spin index a — 1,2,3 referring to the x, y or z component in 
spin space, and with orbital index i — 1, 2 referring to the x or y component in real 
space. (We will always denote vector spin indices — such as a — as superscript to dis- 
tinguish them from ordinary space indices.) In this basis, the ■^He-a order parameter is 
given by [|8|] 

■ 

4" = A I I , (1.147) 

1 ~i 
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where A is a constant whose physical significance will become clear when we proceed. 
Without loss of generality, we have chosen A to be real. 



It follows from (1.146) that the coset space is three dimensional, implying three 



Goldstone modes. To make these degrees of freedom explicit, we write the order pa- 
rameter in the more general form 



= Ad" (e^ - ie^), e 



2itf> 



(1.148) 



where and are two orthonormal vectors in real space, is a unit vector in spin 
space, parameterizing the coset SO(3)'^/SO(2)'^ ~ S'^, and 2lp is the phase of the order 
parameter The residual U(l)^+^ symmetry is reflected in the invariance of the order 
parameter (1.148) under the combined action of a U( 1 ) transformation — > + ^ a and 
a rotation of the e-frame through a about the z axis: (e^— ie^) (e^— ze^) exp(— ia). 
The transformation parameter a may be taken to be spacetime dependent. The specific 
form (1.147) is recovered by choosing rf" in the z direction of spin space, and and 



e in the x and y direction of real space, respectively, and by taking = 0. 



Superfluid ■'He can be modeled by the Lagrangian [ [39| ] 



(1.149) 



where the last term, with Ao a positive c ouplin g constant, is an attractive interaction 
term appropriate for P-wave pairing. In ( 1.149| ), %l^a are Grassmann fields describing 
the '^He atoms with spin a =t, i, and iV) is as defined below ( 1.59| ). Finally, 

di = [di— di)/2kF, with i = 1,2. We again linearize the theory by introducing 
auxiliary fields A and A* using the functional identity 



exp 



DA*DAcxp( -- 



a;; [M-^dM + irA-^Wr] Ko 




(1.150) 



The Euler-Lagrange equation for A shows that classically it merely stands for a product 
of two Grassmann fields: 

A^, = Ao^.(-ia,)V'r. (1.151) 

As may be inferred from this equation, A is symmetric in its spin indices, i.e., it is a 
spin triplet. Moreover, the index i indicates that the field is a vector under spatial rota- 
tions. The partition function Z of the theory can now be represented by the functional 
integral: 



z= j D^-tD* j DA*DA exp j a;;a;^ 



(1.152) 



X exp 



-ig{p,/fcF,A"} 



Pa + ?(P) 



5- 
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where { , } denotes the anticommutator. We recall that the operator p^j_=i 9^ acts on 
everything to its right. The derivatives 9^ A' andc},;A* ' contained in (1.152) arise after 



a partial integration in the corresponding terms in (1.150). The tensor g in (1.152) is 
the metric spinor 

9=( ? I), (1.153) 



-1 



and stands for the Nambu multiplet 



_ f ^1 

9^1 



(1.154) 



The reason for introducing the spinor g in this multiplet is that the two symmetric 
matrices A* (i = 1,2) may be expressed as a linear combination of the symmetric 
matrices cr"(?, with a" {a = 1, 2, 3) the Pauli matrices: 



a;, = A^Kg). 
where Af are the expansion coefficients, so that 



(1.155) 



(1.156) 



In this way, th e metric spinor g disappears from the argument of the exp onential func- 
tion in (1.152). This will simplify the calculation. With the value (1.148) for the matrix 
Af, the Lagrangian appearing in the partition function ( 1.152| ) can be cast in the con- 
cise form 



Po - ^(p)t3 - A(p) - I 



(1.157) 



where (z = 1, 2) are Pauli matrices in Nambu space which should not be confused 
with the Pauli matrices a" in spin space. The vertex A(p) is the same as the one 
we encountered before in the context of the BCS theory [see Eq. (1.95)]. In deriving 
(1.157) we have given the e-frame the standard orientation. 

We proceed in the same way as in the previous section, and decompose the Grass- 
mann fields as in Eq. (1.92). The partition function (1.152) then becomes after carrying 
out the Gaussian integration over the fermionic degrees of freedom 

Z = J DA*DA exp (^iS^^* , A] + ^A;;a;,^ , (1.158) 



with ^cff the effective action. 



A 



Trlog Po - e(p)T-3 - A(p) - — r,{p„ d^}a 



2kv 



(1.159) 



A general orientation of the e-frame, which may vary in spacetime, is accounted for 
by requiring the action to be invariant under the residual U(l)^+^ symmetry. This is 
achieved by the replacement 



B^ip B^ip - • 9^e^ 



(1.160) 
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Since in this section we are not interested in the Goldstone modes of the spontaneously 
broken spin rotation group, we keep the spin vector fixed in the z direction. 

We continue as in the previous section and expand the effective action in a Taylor 
series: 



5eff = i^Tr^i[Go(p)A(p)]^ 



(1.161) 



where the Feynman propagator now reads 



Go(fc) 



K - C(k)T 



3 ^, 3' 



kl - E^{k)+ir] 



A. 



(1.162) 



Here, E(k) is the spectrum of the elementary fermionic excitations 



i?2(k) = e'(k) 



(1.163) 



The constant A introduced in ( 1 . 148 ) is — apart from a factor fcp — seen to be the energy 



gap of these excitations at the Fermi circle. 

The exp ression (1.161 ) can be evaluated in the weak-coupling limit along the lines 
of Sec. L3 to obtain the effective th eory o f the Abelian Goldstone mode [p^. We 
arrive in this way again at the theory ( 1.103 ) we calculated from the BCS theory with 
1/(0) = m/2TT now the two-dimensional density of states per spin degree of freedom 
and n — fc|./27r the two-dimensional fermion number density. It describes a sound 
wave traveling with the speed / \/2. It is remarkable that to this order the additional 
contributions arising from the derivatives in the interaction term of the Lagrangian 
( 1.149 ) cancel so that the same effective theory is obtained as in the previous section 
where these derivatives were absent. 

In the next chapter we shall consider the effective theory describing the antiferro- 
magnetic spin waves associated with the spontaneous breakdown of the spin rotation 
group SO(3)'^ DS0(2)'^ in this particular superfluid '^He phase. Since this symmetry 
is an internal symmetry unrelated to Galilei invariance, there is no reason to expect a 
theory that is invariant under Galilei transformations. 



1.6 Weakly Interacting Bose Gas 

The model commonly used to describe a weakly interacting Bose gas (for a general 
introduction, see the textbooks [Bfl B2I Dl]|), is defined by the Lagrangian [521 



(1.164) 



where e(— iV) — — V^/2to and /xq is the chemical potential. The self-coupling is taken 
to be positive, Aq > 0, so that the contact interaction is repulsive. We will treat this 
model perturbatively in a loop expansion. For this to be applicable the (renormalized) 
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V 




Im(</)) 



Figure 1.2: Graphical representation of the potential energy (1.165). 



coupling constant must be small. Superfluid ^He is a strongly interacting system where 
this premise does not hold. Fortunately, after almost two decades of experimental 
efforts by various groups in the US as well as in Europe, Bose-Einstein condensation 
(BEG) has recently been found in weakly interacting Bose gases. A group in Boulder, 
Colorado were the first to produce such a BEG in a cloud of *^Rb atoms trapped in a 
magnetic field [^3||. The condensate formed at a temperature around 170 nK and was 
comprised of up to 2000 atoms; it was stable for about 15 s. A second unequivocal 
observation of BEG was subsequently established in a system of Na atoms, with a 
condensate containing two orders of magnitude more particles than in the Colorado 
experiments [|^. Various other systems of bosonic alkali atoms are presently under 
study. 

At zero temperature, the global U(l) symmetry is spontaneously broken by a non- 
trivial ground state. This can be easily seen by considering the shape of the potential 
energy 

V = -A*ol0l' + Ao|0^ (1.165) 



depicted in Fig. 1.2 which is seen to have a minimum away from the origin = 0. To 



account for this, we shift by a (complex) constant <j) and write 

0(a;) =e^'^(^) [0 + 0(a;)]. (1.166) 

The scalar field (p{x) is a background field representing the Goldstone mode of the 
spontaneously broken global U(l) symmetry. At zero temperature, the constant value 

Z Ao 
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minimizes the potential energy. It physically represents the number density of particles 
contained in the condensate for the total particle number density is given by 

n{x) = \(t)ix)\^. (1.168) 

Because is a constant, the condensate is a uniform, zero-momentum state. That is, 
the particles residing in the ground state are in the k — mode. In terms of the new 



variables, the quadratic terms of the Lagrangian ( 1.164 ) may be cast in the matrix form 



Co - i|'tAfo(p)l', $ = t , (1.169) 



with 



Moip) = (1.170) 



2 



-2Ao0* -po-e(p) + Mo-C^-4Ao|(/.| 



2 



where U is the Gahlei-invariant combination ( 1.90 . In writing this we have omitted a 
term V^tp containing two derivatives which is irrelevant at low momentum and also a 
term of the form V(p • j, where j is the Noether current associated with the global U(l) 
symmetry, 

Ztm 

This term, which after a partial integration becomes —ipW ■ j, is also irrelevant at low 
energy and small momentum because in a first approximation the particle number den- 
sity is constant, so that the classical current satisfies the condition 

V-j = 0. (1.172) 

The Feynman propagator of the theory is easily extracted by inverting the matrix Mq 
with the background field U set to zero. This yields upon using the mean-field value 



( I1.167D for (j) 

- fcg-i.^(k)+., ^'-'"'^ 

^ / fcoe^*=°" + e(k) + 2Ao|^P -2Xo^^ \ 

\ -2Ao^*' -ko e-'^^on + + 2Xo\4>\^ J ' 

with 7] a small positive constant that is to be set to zero after the integration has been 
performed, and -^(k) the single-particle Bogoliubov spectrum | 



Eik) = Ve'(k) + 2A*oe(k) 



e2(k)+4Ao|<^|2e(k). (1.174) 

The most notable feature of this spectrum is that it is gapless, behaving for small mo- 
mentum as 

E{k)^uoH (1.175) 
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with uq = -y/ ^la/m a velocity which is sometimes referred to as the microscopic sound 
velocity. It was first shown by Beliaev that the gaplessness of the single-particle 
spectrum persists at the one-loop order. This was subsequently proven to hold to all 
orders in perturbation theory by Hugenholtz and Pines [p7|]. For large momentum, the 
Bogoliubov spectrum takes a form 

S(k) - e(k) + (1.176) 

typical for a nonrelativistic particle with mass m moving in a medium. To highlight the 
condensate we have chosen here the second form in ( |1.174| ) where /iq is replaced with 
2Ao|0p. 

Since gapless modes in general require a justification for there existence, we expect 
the gaplessness of the single-particle spectrum to be a result of Goldstone's theorem. 
This is corroborated by the relativistic version of the theory. There, one finds two 
spectra, one corresponding to a massive Higgs particle which in the nonrelativistic 
limit becomes too heavy and decouples from the theory, and one corresponding to 
the Goldstone mode of the spontaneously broken global U(l) symmetry [^]. The 
latter reduces in the nonrelativistic limit to the Bogoliubov spectrum. Also, when the 
theory is coupled to an electromagnetic field, one finds that the single-particle spectrum 
acquires an energy gap. This is what one expects to happen with the spectrum of 
a Goldstone mode when the Higgs mechanism is operating. The equivalence of the 
single-particle excitation and the collective density fluctuation has been proven to all 
orders in perturbation by Gavoret and Nozieres [p^]. 

Given this observation we immediately infer that at low energy and small momen- 
tum the superfluid phase is described by the theory we encountered before 



2m 



(1.177) 



of a nonrelativistic sound wave, with the real scalar field Lp representing the Goldstone 
mode of the spontaneously broken global U(l) symmetry. The alert reader might be 
worrying about an apparent mismatch in the number of degrees of freedom in the nor- 
mal and the superfluid phase. Whereas the normal phase is described by a complex 
(/)-field, the superfluid phase is described by a real scalar field ip. The resolution of 
this paradox lies in the spectrum of the modes [^]. In the normal phase, the spec- 
trum -E'(k) = k'^/2TO is linear in the energy, so that only positive energies appear 
in the Fourier decomposition, and one needs — as is well know from standard quan- 
tum mechanics — a complex field to describe a single particle. In the superfluid phase, 
where the spectrum is i?^(k) ~ c^k^, the counting goes differently. The Fourier de- 
composition now contains positive as well as negative energies and a single real field 
suffice to describe this mode. In other words, although the number of fields is different, 
the number of degrees of freedom is the same in both phases. 



The effective theory is graphically again represented by Fig. 1.1. To lowest order, 



the inverse propagator is given by the matrix Mq in ( 1 . 1 70 ), so that the vertex describ- 
ing the interaction between the U and the ^-fields is minus the unit matrix. In terms of 
the full Green function G, the particle number density now reads 

n = - tr / G{k). (1.178) 



2 
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The (0 0)-component of the full polarization tensor, IIoo, at zero energy transfer and 
low momentum q, 

i limnoo(0,q) = limtr / GTGCfcck + q) (1.179) 

contains a symmetry factor i which is absent in fermionic theories. Following the 
same step s as in the case of superconductors, we again arrive at the compressibility 
sum rule (|l.l09|). 



The diagrams of Fig. |L1| can also be evaluated in a loop expansion (thereby inte- 
grating out the complex scalar field 0) to obtain explicit expressions for the particle 
number density n and the sound velocity c to any given order [|o|]. In doing so, one 
encounters — apart from ultraviolet divergences which can be renormalized away — also 
infrared divergences because the Bogoliubov spectrum is gapless. When however all 
on-loop contributions are added together, these divergences are seen to cancel [|o|]. 
We will in this report not proceed in this way, but instead compute the thermodynamic 
potential fi from which both n and c are obtained by differentiating with respect to the 
chemical potential [see (1.110)]. 

In the approximation ( 1 . 1 69 ) of ignoring higher than second order in the fields, the 
integration over the $ field is Gaussian. Carrying out this integral, we obtain for the 
partition function 



Z = e 



D$ exp [i / Cq] = e 



'/.^°Det-i/2[Mo(p)], 



(1.180) 



where Vo is the potential ( |1.165| ) with (f> replaced by (f>. 



Vo 



4Ao' 



(1.181) 



Writing 



Z = exp 



« / (Vo + Veff ) 



(1.182) 



we find from (1.180) that the effective potential is given to this order by 



Vcff = -7rtr / ln[Mo(fc)]. 

^ Jk 



(1.183) 



The use of a potential rather than an action here is to indicate that we are working with 
a spacetime-independent condensate, so that also U — 0. The easiest way to evaluate 
the loop integral over fc^ is to assume — without loss of generality — that (p is real and 
to first differentiate the expression with respect to /io 



Atr / In P° ~ - ^" 
o^J,o Jk \ -Mo 



-Mo 

6(k) 



k 1^0 



Mo 
e(k) 



(1.184) 
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where i?(k) is the BogoHubov spectrum (1.174). The integral over fco can be carried 
out with the help of a contour integration, yielding 



6(k) 



i?(k)- 



(1.185) 



This in turn is easily integrated with respect to ^q. 
obtain for the effective potential 



£;(k). 



Putting the pieces together, we 



(1.186) 



The remaining integral over the loop momentum k is diverging in the ultraviolet. We 
regularize it by introducing a momentum cutoff A. The model (1.164) is only valid 
far below the cutoff. Above it, the model breaks down and new physics starts. In the 
large-A limit, we arrive at 



Veff 



1 



127r- 



1 

4^ 



Am/iQ 



157r 



3/2 5/2 
2™ ' Mo ' 



(1.187) 



ignoring an irrelevant term (cx A'^) independent of the chemical potential. Equation 
(1.187) contains two ultraviolet-diverging terms. It should be realized that they arise 
from a region where the model (1.164) is not applicable. They cannot therefore be 
of deep significance. As a consequence of the uncertainty principle, stating that large 
momenta correspond to small distances, terms arising from the ultraviolet region are 
always local and can be absorbed by redefining the parameters appearing in the La- 
grangian [pT|. Since /io = 2Ao|0p, the two ultraviolet-diverging terms in (1.187) can 
be absorbed by introducing the renormalized parameters 



1 

67r' 



AoA3 



(1.188) 
(1.189) 



Because the diverging terms are — at least to this order — of a form already present in the 
original Lagrangian, the theory is called "renormalizable". The renormalized parame- 
ters are the physical ones that are to be identified with those measured in experiment. 
In this way, we see that the contributions to the loop integral stemming from the ultra- 
violet region are of no importance. What remains is the finite part 



Veff 



157r- 



(1.190) 



This result could have been obtained directly without renormalization if we, instead 
of introducing a momentum cutoff to regularize the integrals, had employed analytic 
regularization. In such a regularization scheme, where for example the integrals are 
analytically continued to arbitrary real values of the space dimension d, divergences 
proportional to powers of the cutoff never show up as was demonstrated in Sec. 1.4. 



Only logarithmic divergences appear as 1/e poles, where e tends to zero when the 
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parameter in which the analytic continuation is carried out is given its physical value. 
These logarithmic divergences \n{A/E), with E an energy scale, are relevant also in 
the infrared because for fixed cutoff ln{A/E) — oo when E is taken to zero. 

In so-called "nonrenormalizable" theories, the terms which are diverging in the 
ultraviolet are still local but not of a form present in the original Lagrangian. Whereas 
in former days such theories were rejected because there supposed lack of predictive 
power, the modern view is that there are no fundamental theories and that there is no 
basic difference between renormalizable and nonrenormalizable theories [62\. Even a 



renormalizable theory like ( 1.164 ) should be extended to include all higher-order terms 
such as a |(/)|^-term which are allowed by symmetry. These additional terms render the 
theory "nonrenormalizable". This does not however change the predictive power of the 
theory. The point is that when describing the physics at an energy scale E far below 
the cutoff, the higher-order terms are suppressed by powers of E/ A, as follows from 
dimensional analysis. Therefore, far below the cutoff, the nonrenormalizable terms are 
negligible. 

The thermodynamic potential 51 becomes to the order in which we are working 



n= iVo + Veff ) . 



(1-191) 



We are now in a position to determine the particle number density and the sound ve- 
locity using ( |1.11CI| ). We find 



la 
2 A 



8A 
3^ 



(1.192) 



and 



TO 



4A 
3^ 



(1.193) 



where in the last equation we made an expansion in the (renormalized) coupling con- 
stant A which was assumed to be small. These equations reveal that the expansion is 
more precisely one in terms of the dimensionless parameter Am'^/^/x^/^, or reintroduc- 
ing Planck's constant Xm^^"^ ^i^/^ /h'^ . 

Up to this point, we have considered the chemical potential to be the independent 
parameter, thereby assuming the presence of a reservoir that can freely exchange par- 
ticles with the system under study. It can thus contain any number of particles, only 
the average number is fixed by external conditions. From the experimental point of 
view it is, however, often more realistic to consider the particle number fixed. If this is 
the case, the particle number density n should be considered as independent variable 
and the chemical potential should be expressed in terms of it. This can be achieved by 



inverting relation (1.192): 



4:7:na 



l + ^(2nm3A3)V2 



TO 



32 ffia-^ 



(1.194) 
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where in the last step we employed the relation between the (renormalized) coupling 
constant A and the S-wave scattering length a [B3l S2p, 



A = 



2na 



m 



(1.195) 



This relation follows from comparing the differential cross sections for scattering of 
a slowly moving boson from a hard-core sphere of radius a and from a delta function 
potential with strength A in the Born approximation. For the sound velocity ( 1.193 ) 
expressed in terms of the particle number density we find 



2nX 
m 

Anha 



1 + (2nm^X^ 



1 + 16 



V 71" 



1/2- 



(1.196) 



It is important to note that is linear in the coupling constant. Without the interpar- 
ticle interaction characterized by A, the sound velocity would be zero. Moreover, the 
interaction must be repulsive for the system to support sound waves. Let us in this 
connection mention that a third experimental group has reported evidence of BEC in 
a system of ^Li atoms [|64|]. This would be a surprising result as ^Li atoms have — in 
contrast to ^^Rb and Na atoms — an attractive interaction and thus a negative scattering 

length. 

The effective theory ( |1.177| ) can also be put in a equivalent form 



£eff = ~nU{x) + ]u{x)^Uix), 
4 Ao 



(1.197) 



which can be easily generalized to systems with long-ranged interactions. A case of 
particular interest to us is the Coulomb potential 



y(x) 



(1.198) 



whose Fourier transform in d space dimensions reads 

F(k) = 2'^-V('^-i)/2r[i(d-i)] 



(1.199) 



The simple contact interaction Li — — Aq in (1.164) gets now replaced by 



Li = — 
2 



|0(t,x)|V(x-y)|0(t,y)|^ 



(1.200) 



The rationale for using the three-dimensional Coulomb potential even when consider- 
ing charges confined to move in a lower dimensional space is that the electromagnetic 
interaction remains three-dimensional. The effective theory (1.197) now becomes in 
the Fourier representation 



£eft = -nU{k) + l[/(fco, k)-^C/(fco, -k) 



(1.201) 
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and leads to the dispersion relation 



E\k) = 2'^^('^-i)/2r - 1)] ^|k| 



3-d 



(1.202) 



For d = 3, this yields the famous plasma mode with an energy gap given by the plasma 
frequency cjp — Anne^/m. 

To appreciate under which circumstances the Coulomb interaction becomes impor- 
tant, we note that for electronic systems l/|x| ^ /cp for dimensional reasons and the 
fermion number density n ~ /cp, where fcp is the Fermi momentum. The ratio of the 
Coulomb interaction energy ec to the Fermi energy ep = kp/2m is therefore pro- 
portional to n"^/'*. This means that the lower the electron number density, the more 
important the Coulomb interaction becomes. 

For later reference, we close this section by calculating th e fracti on of particles 
residing in the condensate. In deriving the Bogoliubov spectrum ( 1.1 74 ), we set |0p = 
/io/2 Ao . For our present consideration we have to keep <f> as independent variable. The 
spectrum of the elementary excitation expressed in terms of <j) is 



i?(k) = V kk)-Mo+4Ao|^p 



(1.203) 



It reduces to the Bogoliubov spectrum when the mean-field value (1.167) for (/> is in- 
serted. Equation (1.186) for the effective potential is still valid, and so is (1.191). We 
thus obtain for the particle number density 



2Ao 2dfioA 



(1.204) 



|0|2=^o/2Ao 



where the mean-field value for (h is to be substituted after the differentiation with re- 



spect to the chemical potential has been carried out. The integral in (1.204) gives rise 



to an ultraviolet-dive rging t erm which can can be cancelled by go ing over to the renor- 
malized parameters ( 1.189 ) and ( 1.188 ) in the first term of ( 1.204 ). In the second term, 
being a one-loop result, we may to this order simply replace the bare by the (one-loop) 
renormalized parameters. We find in this way [ p^ ] 



or for the so-called depletion of the condensate [ p6| , |67| ] 

,3\ 1/2 



n 
no 



1 



na 



3 V TT 



(1.205) 



(1.206) 



where 



no 



2A 



(1.207) 



is the number density of particles in the condensate. Equation (1.206) shows that even 



at zero temperature not all the particles reside in the condensate. Due to the inter- 
particle repulsion, particles are removed from the ground state and put in states of 
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finite momentum. It has been estimated that in superfluid He — a strongly interacting 
system — only about 8% of the particles condense in the zero-temperature state [^. 
However, it is well known (see next section) that at zero temperature all the particles 
nevertheless participate in the superfluid motion [|69|]. Apparently, the condensate drags 
the normal fluid along with it. 



1.7 BEC at Finite Temperature 

We continue to consider BEC in a weakly interacting Bose gas at finite temperature. 
Some finite-temperature considerations can be found in a textbook by Po pov [[70| ]. At 



the absolute zero of temperature, we saw that the Bogoliubov spectrum ( 1.174 ) van- 
ishes linearly as the momentum goes to zero. Since this property is a direct conse- 
quence of the spontaneously broken global U(l) symmetry, we expect it to hold at any 
temperature below the critical temperature Tc- We shall compute the sound velocity in 
the vicinity of Tc . 

To calculate the thermodynamic potential at finite temperature, we have to general- 
ize the zero- temperature quantum field theory used up to this point to finite temperature. 
For our purposes, this is achieved simply by going over to imaginary time t —it, 
with < r < 1 /T, and by replacing integrals over energy fco with summations over 
frequencies 

' g{ka)^iTY,9{^^n). (1.208) 

kg 



Ik 



where g is an arbitrary function, and ujn are the so-called Matsubara frequencies, 

2mrT for bosons ,^ 

(2n+l)7rr for fermions. 

These rules allow us to calculate the partition function Z at finite temperature. The 
thermodynamic potential is then obtained using the relation 

n = -T\n{Z). (1.210) 

As in the previous section, we evaluate the thermodynamic potential for the case 



where the condensate is spacetime independent. That is, ip introduced in (1.166) is as 



sumed to depend only on temperature and the scalar field ip representing the Goldstone 



mode is set to zero. It was shown in (1.191) that in the Bogoliubov approximation of 



ignoring terms higher than second order in the complex field (p, the thermodynamic 
potential is given by the effective potential. Applying the rules outlined above to the 



zero-temperature expression (1.183), we obtain 



^cff(7^) = ? II / ln[Mo(»c^n,k)], (1.211) 



with the matrix A/o given by ( 1 . 1 70 ) with [/ — 0. The mean-field value of the potential 
energy is given by 

Vo(r) = -fio\HT)\^ + MHT)\\ (1.212) 
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with |^(T) p denoting the number density of particle residing in the condensate. 

The summation over the Matsubara frequencies can be carried out using standard 
methods. In this way, one obtains for the effective potential 



Vcff(r) 



T/^Ml 



-iJ(k)/T 



(1.213) 



The first term at the right-hand side has been analyzed in the previous section and 
was shown to le ad to a renormalization of the parameters Aq and iiq in addition to 
the contribution ( 1.190 ) which we will ignore here. We shall study the second term 
in the vicinity of the critical point by expanding it in a high-temperature series. The 
expansion is justified only if Tc is in the high-temperature regime. It will turn out that 
this is indeed the case for the weak-coupling theory we are discussing. We obtain for 
the thermodynamic potential 



n 

V 



= V(T) 



(4A|^|^ 




(4A|^| 



el' 



2 \2 



(1.214) 



where the integration variable is q ~ |k|/v2mr. The potential is denoted by V{T) 
to indicate that we included the zero-temperature renormalization of the parameters in 



(1.212). The first integral appearing here is finite and yields 

dyy^ In (l _ e~^^) = -ir(|)C(|), 



(1.215) 



with r the gamma function and ( the Riemann zeta function. The last two integrals, 
however, are infrared divergent. We regularize these by analytically continue the fol- 
lowing equations to arbitrary values of a 



dx 



with the result 

n 

V 



dx x" 
~ (e^ - 1)2 

3/2 



r(a)C(«) 



r(a)[C(a-l)-C(«)], 



(1.216) 
(1.217) 



4A2|0|4+i(4A|^|2 



(1.218) 
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To this order in 1 /T the thermodynamic potential is of the Landau form, involving 
terms up to fourth order in the order parameter ^: 



3/2 



with 

a (^-independent term, and 



-C(|)T5/2 _ ^^(|)t3/2 _ 1^2^M )yl/2 



ai 



TO \ 3/2 



a2 = A 



1 - 12A 



(1.219) 
(1.220) 

(1.221) 
(1.222) 



We note, however, that our expansion is not one in terms of since the higher-order 
terms in 1/T contain, besides higher-order terms in 0, also and |0|^-terms. 

The critical temperature is determined by setting the coefficient ai to zero. An 
approximate solution to this equation is given by 



Tr = 



1 2C(i) 



2/3 - 



3 ^(3\^^ 



0(A2/3), 



(1.223) 



where we expanded in the coupling constant A. This justifies the high-temperature 
expansion we have been using because the leading term is of the order A^^/^ which is 
large for a we ak-cou pling theory. 

Equation ( 1.223 ) expresses the critical temperature in terms of /i, the chemical po- 
tential. As we mentioned in the previous section, from the experimental point of view 
it is sometimes more realistic to consider the particle number density as independent 
variable. To express Tc in t erms of the particle number density, we invert the equation 
for n obtained from ( 1.219 ) with cf> — 0, 



n{Tc 



1 dfl 



f mTc 



V 2n 



3/2 



(1.224) 



to obtain /i = When this is substituted in the condition ai = 0, the latter can be 

manipulated in the form 



(T,) - 4n(T,)A {^f" am'' - cm'/' = 0, 



(1.225) 



from which the critical temperature as function of the particle number density can be 
determined. Setting A equal to zero, we obtain the critical temperature 



2tt 

TO VC(|) 



2/3 



(1.226) 
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of a free Bose gas. To first nontrivial order in an expansion in A, we find 

4 ca: 



= Tn 



1 



3^^/3(1) 



n(Te)a3 

C(|) 



1/3- 



(1.227) 



where in the last equation we replaced the coupling constant A with the scattering 
length a using relation (1.195). This equation implies a (slight) increase of the critical 
temperature due to the weak interaction since (^(5) < 0. This is qualitatively different 
from the strongly interacting ''He system. A free gas with *He parameters at vapor 
pressure would have a critical temperature of about 3 K, whereas liquid ^He becomes 
superfluid at the lower temperature of 2.17 K. 



The value of the order parameter (j>{T) obtained from Eq. (1.219) reads 

1 ai 



i0(T)r = 



2^2' 



(1.228) 



Near the critical temperature this can be cast in the form 



3C( 



V 2n 



1 



T 



3n{Tc 



(1.229) 



Recalling that |0(T) p represents the particle number density no{T) of the condensate, 
we see that when the critical temperature is approached from below, the Bose-Einstein 
condensate is drained of particles, and that at Tc it vanishes altogether. 

We are now in a position to investigate the spectrum of the elementary excitations 



at finite temperature. If one simply substitutes the value (1.228) for (j){T) in the Bo- 
goliubov spectrum (1.203) one finds, contrary to what is expected, that it has an energy 
gap. The solution to this paradox, which also exists in the relativistic |(/)|'' -theory [[7T||, 
lays in the observation that, as at zero temperature [|63|], perturbation theory should 
be carr ied out consistently. This has not been don e up t o this point. Whereas Eq. 
( 1.228 ) inclu ded th ermal fluctuations, the spectrum ( 1.203 ) did not. To fix this, we go 
back to Eq. ( 1.221 ) and note that it reflects a change in the chemical potential due to 
thermal fluctuations. With this modification the finite-temperature spectrum becomes 
gapless — in accordance with a theorem due to Hohenberg and Martin [f7^]. 

Let us work this out in some detail near the critical temperature, where |0(r)p is 



given b y ( 1.229 ). Retaining only the leading term in the 1/T expansion, we infer from 
( 1.221 ) the following change in the chemical potential 

2tt J ' 



/i^/i = /i-4C(|)A 



(1.230) 



and 2A|(/)p = fl, where fi is given a tilde to indicate that it is dressed by thermal 
fluctuations. In this way, the finite-temperature spectrum becomes 



i?(k) = Ve'(k)+4A|0(T)|2e(k), 



(1.231) 
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which is indeed gapless. Since we included the zero-temperature renormalization of 



the parameters stemming from the first term at the right-hand side of (1.213), the 



renormalized parameters A and /i feature in these equations. The gaplessness of the 
finite-temperature spectrum was shown to be true in all orders of perturbation theory 
by Hohenberg and Martin [^]. The result we just obtained is completely analogous 



to the zero-temperature Bogoliubov spectrum (1.174). The expression for the sound 



velocity at finite temperature we extract from the spectrum reads 



c\T) 



-\m\' 

m 



(1.232) 



It vanishes when the temperature approaches Tc- This is in accord with the observation 
that the gapless Goldstone mode of the spontaneously broken global U(l) symmetry 
vanishes at the critical temperature. We note that also at finite temperature, the sound 
velocity squared is linear in the coupling constant — without the interparticle repulsion 
the sound velocity would be zero. 

At the end of the preceding section we mentioned the well-known fact that despite 
only a fraction of the particles at zero temperature reside in the zero-momentum state, 
all particles participate in the superfluid motion. To show this, let us assume the entire 
system moves with a velocity u relative to the laboratory system. As is known from 
standard hydrodynamics the time derivate in the frame following the motion of the fluid 
is 9o + u • V [see Eq. (1.8)]. If we insert this in the Lagrangian (1.164) of the near-ideal 
Bose gas, it becomes 



C = ri^do ~ e(-^V) + Mo - u • (-^V)]0 - Ao|</)|^ 



(1.233) 



where the extra term features the total momentum iV)0 of the system. The 

velocity u multiplying this is on the same footing as the chemical potential which 
multiplies the particle number Whereas /ip is associated with particle number 

conservation, u is related to the conservation of momentum. 

In the two-fluid picture, the condensate can move with a diffe rent ve locity Vg as the 
rest of the system. To bring this out we introduce new fields, cf. ( 1.166 ) 



(1.234) 



(l>{x) ^ (t)' (x) =e" 

in terms of which the Lagrangian becomes |^3| 

C = r [*ao-e("^V)+A^o-imv,•(v3-2u)-(u-v,)•HV)](/.~Ao|0|^ (1.235) 

where we dropped the primes on (p again. Both velocities appear in this expression. 
Apart from the change u — > u — Vg in the second last term, the field transformation 
resulted in a change of the chemical potential 



Mo Mcff := Mo - ^mvs • (v^ - 2u) 
where /ioff may be considered as an effective chemical potential. 



(1.236) 
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The equations for the BogoHubov spectrum and the thermodynamic potential are 
readily written down for the present case when we keep these two changes in mind. In 
particular, the effective potential (1.213) reads 



[£;(k) + (u-v,).k]/T 



(1.237) 



where E(k) is the Bogoliubov spectrum (1.203) with the replacement Eq. (1.236). 
The mean-field potential Vq (T) is given by ( 1.212| ) with the same replacement. The 
momentum density, or equivalently, the mass current g of the system is obtained in this 
approximation by differentiating the potential Vo{T) + VosiT) with respect to — u. We 
find, using the equation 



du 



that it is given by 



g — mnrVs 



kexp{[ii;(k) + (u-v,).k]/r}-i- 



(1.238) 



(1.239) 



The last term is the contribution stemming from the elementary excitations. In the zero- 
temperature limit, this term vanishes and g = mnvg. This equation, comprising the 
total particle number density n, shows that at zero temperature indeed all the particles 
are involved in the superflow, despite the fact that only a fraction of them resides in 
the condensate [E^. When the condensate moves with the same velocity as the rest 



(Vg — u), the last term in ( 1 .239 ) vanishes again, now by symmetry. 

Assuming that the difference between the normal and superfluid velocity is small, 
we may expand the last term in (1.239) to linear order in this difference to find 

gB(k)/T 



g = PVs 



1 

3T 



(u- V,), 



■ (e£(k)/T_i)2 

where p = mn is the total mass density of the fluid, and we used the general result 



(1.240) 



(1.241) 



with /(|k|) an arbitrary function depending only on the length of k. The last term stem- 
ming from the elementary excitations may be used to define the normal mass density 

Pn, 



Pn 



1 

3T 



=B(k)/T 



(eS(k)/T _ 1)2 



Writing 



P = Pb 



Pn 



for the total mass density, we may cast ( |1.240| ) in the form 

g = PsVs + p„Un. 



(1.242) 



(1.243) 



(1.244) 
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These last two equations are the basic equations of the two-fluid model [|69|]. The 
model was introduced by Tisza [ [74| | using idea's of London to give a phenomenolog- 
ical description of superfluid '*He. It not only successfully explained various startling 
experimental properties of the strongly interacting system, but also predicted new phe- 
nomena which were later confirmed by experiment 



1.8 BCS at Finite Temperature 



We in this section apply the high-temperature expansion to the BCS theory to show 
that it yields the usual Ginzburg-Landau theory [ [75[ ) . 

The one-loop effective potential for a uniform system at finite temperature is readily 
seen to take the form 



Veff = 



[ E{k) -2T [ \n 



1+C 



-B(k)/T 



(1.245) 



with i? (k) the BCS spectrum ( 1 .79) of the elementary fermionic excitations. The factor 
—2 in (1.245) arises because there are two fermion species, with spin a =1 ™d |, 



respectively. To the one-loop potential we have to add the tree contribution 



Vo 



Ao 



(1.246) 



The form ( 1 .245 ) for the effective potential is usually not taken as starting point to 



derive the Ginzburg-Landau theory. Normally, one first carries out the integration over 
the loop energy ^ and then performs the sum over the Matsubara frequencies [[76|]. 
But here we follow — in analogy with the previous calculation — the opposite route and 
have carried out the summation first and will now perform the integration over ^ (in the 
weak-coupling limit). 



Veff = 



-2Ti.(0) J ln(l + c-^/^ 



(1.247) 



where we ignored the first term in (1.245) corresponding to the one-loop quantum 
contribution to the effective potential which we studied in Sec. pTj] , save for the renor- 
malization of Ao it leads to. Our problem is seen to reduce to one in a single dimension. 
To cope with divergences that will appear we dimensional regularize the integral and 
consider the problem in 1 — e dimensions. Equation (1.247) then becomes 



Veff = -4 T^i^iO) Ay y-' In (l + g- V^^WT^^ , (1.248) 

with y the dimensionless variable y — (,/T. We kept the one-dimensional integration 
measure Ki, where 

Kd = T (1.249) 

(4^)d/2p(d)' 
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is the area of a unit sphere in d spatial dimensions divided by {2ttY. Because A 
appears only in the combination | Ap/T^, it follows that an expansion in this parameter 
is tantamount to one in 1/T. Carrying out this expansion, we arrive at 



eflF 



2u{0) / 
Jo 



dyy~ 



-2T^\n{l + e-y)+\A\^- ^ 



yey + 1 




(1.250) 



In analogy with what we did in the case of a weakly interacting Bose gas, we regular- 
ize the infrared-divergent integrals by analytically continue the following equations to 
arbitrary values of a 



f 



dx ln(l + e-") = — , 



dx x" 



r(a)(l 



^1— a 



)C(a) 



Jo 



dx 



/o X (e^-M)2 
In this way we obtain in the limit e — > 



X + 1 

= r(a) [(1 - 2i-")C(a) - (1 - 22-«)C(a - 1) 



(1.251) 
(1.252) 
(1.253) 



Veff 



KO)|A|^ 
A|4 



-+7 + 21n(2) + 2C'(0) 



-JC'(-2)K0) ^2 • 



To make contact with the standard approach we use the identities 



(1.254) 



C'(0) = -iln(27r), 



C'(-2) = -^C(3) 



and substitute 



e V T 



(1.255) 



(1.256) 



where the Debeye energy uju, being a measure of the inverse lattice spacing, is the 
physical ultraviolet cutoff and the temperature T is the relevant infrared scale. [This 
correspondence between the pole 1/e of dimensional regularization and the logarithm 
In (A) appearing in the regularization with a cutoff is commonly used in the context of 
quantum field theory.] The critical temperature is again determined by the condition 
that the coefficient of the term quadratic in the order parameter be zero, i.e.. 



■7-ln(i7r) 



1 

A=0' 



(1.257) 
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where we included the tree contribution ( 1 .246 ) with the bare coupling Aq replaced with 
the renormalized one (1.84). This yields the standard result for the critical temperature 

Te = -e'^WDe^/'^(°)^, (1.258) 

TT 

where it should be kept in mind that the coupling constant A is negative in the weak- 
coupling limit. Employing this expression for the critical temperature, we can ma- 
nipulate the effective potential (1.254) with the tree contribution (1.246) added in the 
canonical form 



V 



v{0) In 



T2 



(1.259) 



valid close to the critical temperature. The minimization of this potential with respect 
to A yields the well-known temperature dependence of the order parameter near the 
critical temperature 



my 



-T' 1 



(1.260) 



which should be compared with Eq. (1.229) we obtained for a superfluid. 
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Chapter 2 

Induced Quantum Numbers 



Solitons play an important role in condensed matter physics (for a general introduction 
see Ref. p% 78 1). One of the most famous solitons is the magnetic vortex in a super- 



conductor first discussed by Abrikosov [|79|]. In various models, solitons are found to 
have peculiar quantum numbers associated with them. We, in this chapter, exclusively 
study solitons arising in fermionic systems. The unusual quantum numbers induced 
by the solitons can then be computed by integrating out the fermionic degrees of free- 
dom. The effective theory thus obtained contains the same information as the original 
fermionic model at low energy and small momentum. In particular, approximate ex- 
pressions for the fermionic currents — so-called Goldstone-Wilczek currents [ |80| ] — can 
be extracted from it. To see how this is connected to the quantum numbers of solitons, it 
should be realized that the effective theory and the ensuing Goldstone-Wilczek currents 
are built from background fields, which in our definition include possible Goldstone 
fields. But solitons are precisely specified by these fields. Hence, when a specific field 
configuration describing a soliton is substituted in the Goldstone-Wilczek currents, the 
induced quantum numbers localized on the soliton can be evaluated. These quantum 
numbers arising in purely bosonic effective theories can be peculiar in that they can 
have fermionic, or even anyonic characteristics. We will encounter various examples 
of this along the way. 

We will also consider phase transitions where the induced quantum numbers local- 
ized on the solitons change from fermionic to bosonic. Such a transition, first consid- 
ered by Wen and Zee [^1]], is called a statistics-changing phase transition. 

2.1 Skyrmion Lattice 

In this section we investigate the effective theory describing the spin degrees of free- 



dom of superfluid '^He-a. In this superfluid state, as we remarked in Sec. 1.5, the 
SO'^(3) spin rotation group is spontaneously broken to the group SO'^(2) of spin ro- 
tations about the preferred spin axis. Since the state is characterized by the magnetic 
quantum number ms — 0, we expect it to display besides superfluid also antiferromag- 
netic behavior. The spin waves of an ordinary uniaxial antiferromagnet, which are the 
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Goldstone modes of the spontaneously broken SO'^(3) symmetry, are known to be de- 
scribed by the 0(3) nonHnear sigma model [p^. In 2+1 dimensions, the model allows 
for a topological term in the action — the so-called Hopf term [ pi] ] . 

It has been suggested by Dzyaloshinski, Polyakov, and Wiegmann [|3|] that the 
Hopf term with 6 — n should be included in the effective theory describing a spin-^ 
Heisenberg antiferromagnet in 2+1 dimensions. But subsequent microscopic calcula- 
tions showed that the topological term cannot be derived from the Heisenberg model 
It will turn out that in ''He-a the situation is different in that here the quantum- 
induced effective action does accommodate such a topological term. 

To de rive th e effective theory, we use the well-known fact that the unit spin vector 
d°'{x)m{ 1.148 ) can always be rotated in the third direction by introducing a spacetime- 
dependent 2x2 SU(2)-matrix s{x): 



We next introduce the decomposition 



sgx* 



(2.1) 



(2.2) 



in the Lagrangian (1.157). The theory then takes the form of an SU(2) gauge theory: 

£ = i^t fpo - So - Cfe - S,)t3 - ^{{v^ - B,)T,,a^}\ X, (2.3) 



2fcp 



where X — {x^ <7X*)^' with the superscript T indicating the transpose, and the 2x2 
matrix- valued field 

(2.4) 



plays the role of a gauge field. Remember that 9^ = {do, —V). Also recall that the 
Pauli matrices T3 and (i = 1, 2) operate in Nambu space, whereas the Pauli matrices 
cr" (a = 1, 2, 3) operate in spin space. After integrating out the fermionic degrees of 
freedom, we now obtain instead of (1.159) the effective one-loop action 



'S'eff 



Trlog po-So-e(p-B)T3 



A 

2k^ 



{pi - Bi,(j^}n 



Expanded in a Taylor series this becomes, apart from an irrelevant constant, 

oc 

^eff = ^Tr^-[Go(p)A(p)]^ 
where the vertex A(p) is now given by 

A(p) = Bo + :^[{Bf - 2p • B - *V • B]t, - ^Bfn, 
Zra kf 



(2.5) 



(2.6) 



(2.7) 



and Gq{p) is thepropagator(1.162). We have written B • B as (B)^ so as not to confuse 
it with the second spin component (a — 2) of the spatial vector B". 
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The first term S^,^ in the expansion (2.6) is readily shown to give 



c(l) _ 



4m 



tr / (B)= 



(2.8) 



The trace tr stands for the trace over the remaining sigma matrix indices, where one 
should bear in mind that the B^'s defined by (2.4) are elements of the SU(2) algebra. 
For the second term S^g in the expansion ( |2.6[ ) we obtain to lowest order in derivatives: 

Sil' = JkO) tr l^iBof ~ {a^Bof + f [(B)^ + ^{a^Bf] | , (2.9) 



where 7 stands for 



7 



1 



In 



cA 



(2.10) 



with A an ultraviolet energy cutoff and c an irrelevant numerical constant. The term 
proportional to (A//i)^ stems from the last term in the vertex (2.7). Since A << 
/i w ep in the weak-coupling limit we are considering, the factor (A//i)^ In ( cA/A ) is 
negligible compared to 1, and may be set to zero. Adding the contributions (2.8) and 
(2.9), and using the identity 



<a\2 



(2.11) 



which is easily derived from the definitions (2.1) and (2.4), we conclude that, apart 
from a possible topological term, the effective theory describing the antiferromagnetic 
spin waves is the 0(3) nonlinear sigma model | 

^cff = \m [{dod^ f - ^vUd^n^] . (2.12) 

It is gratifying to see that the B-terms obtained from the first and second term in the 
expansion of the effective action conspire to precisely generate the right combination 
tr[(B)2 - [a^Bf]. 

The 0(3) nonlinear sigma model can be equivalently represented by the CP^ model 



(2.13) 



adopting units such that the spin-wave velocity is unity vp / \/2 = 1. Here, the complex 
scalar z = {zi^Z2)'^, which is subject to the constraint z'^ z — 1, is defined by the 
equation 

(2.14) 



and is related to the 2x2 matrix s in the following way 



zi 

Z2 



(2.15) 



In terms of z, the field B,, reads 



Bl = -izta^z. 



(2.16) 
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Since (P defines z only up to a phase, the Lagrangian (2.13) has a local gauge symme- 
try: 

^(a;) ^ e*"(^)z(a;). (2.17) 

In the original variables, such a transformation corresponds to a spin rotation about 
the preferred spin axis. That is, the unbroken component of the global SO'^(3) spin 
rotation group becomes a local gauge symmetry in the effective theory. We will refer 
to this symmetry as local spin gauge symmetry. The reader is referred to Ref. [^sjl for 
a general discussion on this phenomenon. 

We next turn to the topological term. When written in terms of the spin gauge field 
B^, the Hopf term takes the form of a Chern-Simons term: 

OSes ^ ^ Jj^.xBld.Bl (2.18) 

where 6' is a real parameter. 

The question whether a Chern-Simons term is induced at the quantum level can 
be addressed in great generality. To the author's knowledge. So [^; Ishikawa and 
Matsuyama [ |87| ] were the first to reveal the underlying structure, showing that 9 is 
determined by the Feynman propagator G{k) of the underlying fermionic theory 

1 f { ^dQ-^ BQ-^ aG-^\ ^^^^^ 



They also showed that the right-hand side of ( 2.19| ) can acquire a topological meaning 



This work appears to have been largely ignored in the literature, and the results have 
later been rederived by others [|8l p9t 



For ■'He-a, where according to Eq. (1.162) the inverse propagator reads 



G^\k) = ko- e(k)r3 - -^k • Ta\ (2.20) 

Kp 



the 0-parameter takes the value [ P0[ | 

9 = sgn(ei X e^)TT. (2.21) 

Here, sgn denotes the signum, and it should be realized that the cross product of two 
vectors is a scalar in two space dimensions. The dependence of 9 on the orientation of 
the e-frame is crucial because it changes sign under both a parity transformation and 
time inversion. Under a parity transformation in 2h-1 dimensions, one spatial coordinate 
is reflected: 

{t,xi,X2) ^ {t,-xi,x2); (2.22) 
and an arbitrary vector Vp, transforms as: 

Vo,2it,Xi,X2) ^ Vo^2it,-Xi,X2), Vi{t,Xi,X2)^-Vi{t,-Xi,X2)- (2.23) 

The sign change in 9 offsets the one in the Chern-Simons term Ccs under these trans- 
formations, so that the combination 9 Ccs is invariant | |8^ , |90| ]. 
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The appearance of a Chern-Simons term in the effective theory of a nonrelativis- 
tic model might seem surprising. However, the two possible terms eijEgdiBj and 
CijBfdoBj one can write down are related by the requirement of invariance under the 
spin gauge transformation B^ — ^ B^ + d^a, in such a way that the two can be com- 
bined into a single Chern-Simons term. The presence of both terms can be checked 
explicitly by invoking the derivative expansion we used before. The specific form of 



the interaction, reflected in the last term of the vertex (2.7), as well as of the propagator 



(1.162), are crucial in obtaining them. 

It was pointed out by Wilczek and Zee [ ^ ] that the Hopf term imparts spin to the 
solitons of the (2H-l)-dimensional 0(3) nonlinear sigma model. These solitons, first 
discussed by Belavin and Polyakov [p2[], belong to the class of topological defects that 
are regular throughout coordinate space. That is, they do not possess a singular core as 
is often the case. They are characterized by the topological charge [ |77[ ] 

Q = -L /■ e,,e''P''d''dd^djd\ (2.24) 

where the spin indices a, /3, 7 run over 1, 2, 3, and the antisymmetric Levi-Civita sym- 
bol is defined such that e^^^ — 1. We repeat that vector spin indices are always repre- 
sented by superscripts. 

The natural language for the description and classification of defects in ordered 
systems is provided by homotopy theory. (For general accounts on homotopy theory 
see, for example, [^, |Q] or [p5|].) The way to characterize a general defect of 
dimension ed is to surround it by a spherical surface of dimension r such that | 



r = d - Ed - 1, (2.25) 

with d the dimension of the medium under consideration. The minus one at the right- 
hand side represents the distance from the defect to the surrounding hypersphere. In 
each point of this surface, the Goldstone fields define a map of the r-sphere S"^ onto 
the coset space G/H. Topologically stable defects arise when the contour in the 
coset space is not contractable. They are therefore classified by the homotopy groups 



For the case at hand, the topological charge ( |2.24| ) is the winding number of the 
map 

d"(x) : ^ (2.26) 

of compactified space into the internal two-sphere parameterized by d", where 
it is to be noted that the second homotopy group 712(8^) = Z. The corresponding 
topological current, which is conserved independently of the field equations, is 

J^. = ^€,,,^e"P-^d''~d,dP~dxd\ (2.27) 

OTT 

This current may alternatively be written in terms of the field B^ as: 

= ^G^, (2.28) 
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where the dual field is defined by 

= e^,.xd,Bl (2.29) 

The simplest soliton has Q = 1 and is called a skyrmion. Due to the Chern- 
Simons term, skyrmions acquire fractional spin 9/2tt and statistics 9 p7[ |. The 
values 6* = ±7r found in superfluid -^He-a imply that here the skyrmions have spin ^ 
and are fermions. The fractional statistics (for an introductory review, see Ref. [Q) 
imparted by the topological term to the skyrmion can be easily understood by rewriting 



the Chern-Simons term (2.18) as 



zvr 



where is the topological current (2.27). This shows that the effect of the Chern- 
Simons term is to couple the topological current to the spin gauge field B^, with a 
charge 9/2tt. In other words, besides carrying flux, skyrmions also carry charge. Now, 
if a skyrmion winds around another skyrmion, its wavefunction will acquire a phase 
exp(j6') through the Aharonov-Bohm effect, thus turning skyrmions into anyons with 
statistics parameter 6. 

A skyrmion configuration is given by 

(cos(/) sinM(r) \ 
sm(t> sinu(r) , (2.31) 
cosu(r) J 

where (r, (/>) are circle coordinates in the spatial plane, while u{r) is a function with 
u{r) = at the origin of the skyrmion (r = 0) and u{r) vr for r ^ oo. The 



corresponding spin gauge-field configuration Bf^, which can be found by solving (|2.1|) 



for s, reads 

Bl = + i(l - cosu)9^<^, (2.32) 

where a is an arbitrary gauge parameter originating from the circumstance that s is 
determined by d" only up to a phase factor exp(iQ;(T^). We fix the gauge by demanding 
that the "magnetic" flux 

$ = / Go (2.33) 

J X 

is regular everywhere. Since eijdidj(p is singular at the origin, we choose a = 0. This 
leads to 

Go = —-^smu. (2.34) 
2r ar 

The corresponding magnetic flux $(i?) through a disk of radius R is 

= 7r[l - cosw(i?)]. (2.35) 

With m(oo) = TT, it follows that the flux piercing the plane is 2tt. That is, each skyrmion 
carries one unit of flux. As pointed out by Huang, Koike, and Polonyi [p^, these 
conclusions may be nicely visualized by placing a Dirac monopole with unit magnetic 
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Dirac string 




monopole 



Figure 2.1: Compactified space represented by a two-sphere with a monopole at its 
origin. The Dirac string pierces the plane at the center of the skyrmion. 



charge outside of the system at the origin of a two-sphere representing the compactified 
spatial plane, and letting the Dirac string pierce the plane at infinity. Alternatively, one 
may extract a unit flux tube so that the Dirac string pierces the plane at the origin where 
the center of the skyrmion resides (see Fig. 2.1). 

Let us mention that we are working in the ordered phase of the 0(3) nonlinear 
sigma model. The (2H-l)-dimensional model £ = {df^d")^ /2g'^, with g the coupling 



constant, has a phase transition into a disordered, strong-coupling phase at [100| 



27r2 



(2.36) 



where A is an ultraviolet energy cutoff. Since in our case, according to (2.12), g — 



i^(0)/2 = m/47r, with m ^ GeV >> A, we are indeed in the ordered, weak-coupling 
phase. The two phases differ in the way the so-called flux symmetry generated by 
the topological current (2.28) is realized [101 1. In the weak-coupling phase the flux 
symmetry is unbroken, while in the strong-coupling phase it is realized in the so- 
called Kosterlitz-Thouless mode. In the latter phase, there is — at least in perturba- 
tion theory — algebraic long-range order due to the presence of a gapless Kosterlitz- 
Thouless boson. It has been pointed out by Huang, Koike, and Polonyi |^], that a 
condensation of the Dirac monopoles would invalidate this picture. If this happens, 
they argued, the flux symmetry becomes anomalous and the Kosterlitz-Thouless boson 
acquires an energy gap so that the algebraic long-range order of the ordered, strong- 
coupling phase is lost. These authors found numerical evidence supporting this sce- 



nario, but Bitar and Manousakis [ 102| criticized their numerical analysis and concluded 



that the monopoles do not spoil the perturbative picture. 

We next study the connection between the microscopic fermionic model and the 
effective theory, i.e., the nonlinear sigma model with a Chern-Simons term added. In 
particular, we will be interested in how the spin currents of the original model 



Jo 



1 



2m 



(2.37) 
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with Jq the spin density, are represented in the effective model. In ( 2.37 ), di—di—di, 
is the gradient operator acting to the right and left. Since the global SO'^(3) spin rota- 
tion group is spontaneously broken to SO'^(2) in '^He-a, the two components with 
a — \.2 are supercurrents, meaning that the corresponding charges are transported by 
the spin waves without dissipation. We shall demonstrate that the induced spin currents 
are represented by the isospin currents of the nonhnear sigma model. 



dL 



(2.38) 



and that due to the presence of the Chern-Simons term, the spin currents acquire an 
additional contribution. 



Let us rewrite the spin currents (2.37) in terms of x and the fields 



where ^jj and x related via ([2.2|). We find 

JO 

3 



1 



2™ 



xH<y''HK)x-x^B^x 



(2.39) 



where R"'^ is an orthogonal rotation matrix defined via 



(2.40) 



If we differentiate the Lagrangian ( |2.3| ) with respect to the gauge fields, we obtain 
1 dC 



2dBg 
1 dC 



2dB" 



1 

2m 



X^a'^{-td,)x-X^B'^X 



(2.41) 



where in the last equation with a = 3 we neglected the contribution arisin g from the 
interaction term. This is justified because A/fcp << 1- It follows from ( 2.39 ) and 
(^4l|)that 



2 dBf 



(2.42) 



If we now wish to calculate in the effective theory, we have to replace C in this equa- 
tion with the effective Lagrangian. Omitting the Chern-Simons term for the moment, 
we obtain in this way 

= -z.(0)i?"^i?^, (6-1,2). (2.43) 

We proceed to show that the right-hand side of this equation is the isospin current 
(2.38). To this end we write: 

= iw(0)tr sV"9^s + ii^(0)trs V"scr3B^c^^ (2.44) 
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where we used (2.1) and the definition of i?^. Then, employing (2.40) and performing 



the trace, we see that the induced spin currents are indeed given by the isospin currents 



of the nonlinear sigma model [ 103 1 



■a ^ ja 

J fl /J, 



(2.45) 



We next turn our attention to the anomalous contribution to the spin currents arising 
from the Chern-Simons term: 



47r2 



Tja3 



(2.46) 



In particular, we are interested in the spin carried by a skyrmion. The zero-component 
of the dual field describing a skyrmion configuration is given by Eq. (2.34). To find 
the spin, we calculate the projection of the spin density jg^ onto the spin-quantization 
axis. Usually this axis is fixed in spacetime, but for a skyrmion configuration the spin- 
quantization axis (i"(x) varies in space. The spin cr of a skyrmion is therefore given 

' d'^jU (2.47) 



(2.48) 



This quantity is most easily evaluated when use is made of the identity 



Q/3 



which is readil y de rived by multiplying the definition of the R matrix, Eq. ( 2.40 ), with 
and using (2.1). We find in this way that the spin of a skyrmion is related to its 
topological charge ( 2.28 ) via 

a = ^Q. (2.49) 

More precisely, a skyrmion with topological charge Q ~ 1 has a spin 6'/27r in accord 

with its statistics. 

This result can also be derived by directly employing ( 2.42 ) with ( 2.48| ): 



1 dC 



2dBl' 



(2.50) 



and replacing C with the Chern-Simons term. 

Equation (2.49 ) is the anomalous contribution to the total spin of the system. Hence, 
the difference in the number N^,i of spin-| and spin-^ ^^He atoms is given by the 
skyrmion number [ 103 ] 



N^-N^^ ±Q. 



(2.51) 



It follows that a finite number of skyrmions can be created by taking N-^ 7^ ^j,- If 
the difference \N-r — A'^x | is large enough the skyrmions, which are fermions due to the 



Chern-Simons term, may form a hexagonal lattice (see Fig. 2.2) analogous to the one 



formed by so-called Anderson-Toulouse-Chechetkin vortices |104] in rotating super- 
fluid '^He- A [ 1 05 1 . Experimentally, however, such a skyrmion lattice may be difficult to 



realize because an external magnetic field necessary to spin polarize the system locks 
the spin-quantization axis to a plane normal to the field. 
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Figure 2.2: Hexagonal spin-^ skyrmion lattice. The arrows denote the direction of the 
spin vector d" (x) . 



2.2 Peierls Instability 

In this section we shall discuss a model describing quasi one-dimensional metals. The 
effective low-energy, small-momentum continuum theory is shown to have two phases 
both of which possess a rich soUton structure. The two phases differ in a topological 
term which is induced in one phase, but not in the other. This so-called 6'-term is shown 
to change the statistics of certain solitons. 

One-dimensional metals are frequently discussed in terms of the Hubbard model 
defined by the lattice Hamiltonian | 



i?H = -t^{c]„Cj+i^^ +li.c.) + U^nj^nji. (2.52) 



3 



It contains two parameters, U and t, representing the Coulomb interaction between a 
pair of valence electrons on the same atom, and the matrix element for the hopping 



of an electron to a neighboring atom, respectively [106|. When U is chosen positive, 
the Coulomb interaction is repulsive. In (2.52), the operator cj^ creates an electron of 



spin a{=], [) at the location of the jth atom, and rija = c)-^Cja is the electron-number 
operator at site j. The sum is over all lattice sites. It is assumed that the Coulomb 
interaction is highly screened, with a screening length of the order of a lattice spacing 
a, so that we can restrict ourselves to the on-site Coulomb repulsion. In addition, it is 
assumed that the electrons are well localized, so that the tight-binding approximation 
is vaUd. 



The Hubbard Hamiltonian ( [2.52| ) may be rewritten in the equivalent form 



i/H = -<5](4c,+i,. +h.c.) - hjY, {Sf)\ (2.53) 



3 



where we dropped an irrelevant constant. Here, 5" — ^Cj^a^^Cjr, with <t" {a = 
1, 2, 3 ) the Pauli matrices, is the electron spin operator at site j. The representation 



( |2.53[ ) which explicitly exhibits the spin operator is most convenient to derive the ef- 
fective theory. 

In the large-C/ limit, the Hubbard model at half-filling, where the number of elec- 
trons equals the number of sites, can be mapped onto the spin-i Heisenberg model 
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H^jysrsr+i, (2.54) 



with coupling constant J = At^ /U > 0. Because J is positive, spins on neighboring 
sites favor antiferromagnetic coupHng. This model, solved by Bethe, is known to be 



gapless [107|. In the opposite limit of weak coupling U << t, the Coulomb repulsion 
is a small perturbation and we expect the model to behave qualitatively the same as 
a free-electron gas. These two arguments taken together make it plausible that the 
one-dimensional Hubbard model at half-filling is gapless for all values of the coupling 
constant, implying that it exhibits metallic behavior 

Although the Hubbard chain is meant to describe a one-dimensional metal, it com- 
pletely ignores the electron-phonon interaction. This coupling, particularly in one di- 



mension, can however be very important. A famous argument due to Peierls |108] 
shows that a one-dimensional metal with a partially filled conduction band is unsta- 
ble towards a periodic distortion of the linear lattice which opens an energy gap at the 
Fermi points. Because of this, Peierls concluded that such a one-dimensional system 
at the absolute zero of temperature would probably never have metallic properties and 



would instead be an insulator. Frohlich [109| showed that a periodic lattice distor- 



tion is accompanied by a modulation of the electron number density with the same 



periodicity — a so-called charge density wave [110]. 

Being a fundamental property of one-dimensional metals, we wish to extend the 
Hubbard model so as to account for the Peierls instability. At half-filling, the periodic 
distortion has a period of twice the lattice spacing. The modulation of the electron- 
number density resulting from this instability is given by 

i/p = A^(-l)%,,, (2.55) 

jo- 

where A is the so-called Peierls energy gap, it is typical of the order of 100 K. By 



adding this term to the Hamiltonian ( 2.53 ) we account for a crucial aspect of the one- 
dimensional electron-phonon interaction in the Hubbard model. 

We examine the model in the limit of low energy and small momentum by going to 
the continuum. This is faciUtated by introducing the operators a^, Oo with engineering 
dimension i via 

C2j i^-' \/2a ac(2j) 

C2J+1 i^'+^V2^ao{2j + l), (2.56) 

for even and odd sites, respectively. Here, as well as in the following, we suppress 
spin indices. Summation over these hidden indices is always implied. The powers of 



i in (2.56) represent the factor exp (ikFXi), with xi — 2ja for even sites and xi = 
{2j + l)a for odd sites, respectively, where kp — 7r/2a is the Fermi momentum at 
half-filling. In terms of the new operators, the hopping term of the Hubbard model 



becomes 



(c]c,-+i + h.c.)-. -2itaJ2{ K(2j) - al(2j + 2)] ao(2j + 1) 
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+at(2j + l)[a,(2j + 2)-ae(2j)]}. 



(2.57) 



With the derivative of the operator Oo defined by 

ao(2j + 1) - ao(2j - 1) 



diGo = lim 



2a 



(2.58) 



and a similar definition for dia^, we obtain for this term in the continuum limit 



t^(c]cj+i +h.c.) ^ 2ta / i}''a{-idi)i;. (2.59) 



Here, we introduced the multiplet containing the even-site and odd-site species of 
fermions, 



(2.60) 



with a = Ti the first Pauli matrix which should not be confused with the spin matrices 
a". In deriving ( 2.59 ) we replaced the summation over the lattice sites with an integral: 
2aJ2j ^ J dxi. The right-hand side of ( 2.59 ) is of the form of a massless Dirac 
Hamiltonian. It describes Bloch electrons with a gapless dispersion relation that is 
linear in the crystal momentum ki 

E{h)^VF\h\, (2.61) 

where vp — 2ta is the Fermi velocity. A one-dimensional model described by this 
dispersion relation would exhibit metallic properties. As expected, the term (2.55) 
representing the Peierls instability drastically changes this behavior In the continuum 
limit, this term becomes a mass term in the Dirac theory 

A^(-l)^n, ^ 2aA^ [at(2j)ae(2j) - at(2j + l)a,{2j + 1)] 



(2.62) 



with f3 ~ being the diagonal Pauli matrix. As a result, the gapless spectrum (2.61) 
changes into a spectrum with an energy gap 



and the metallic properties of the free model are lost. 



(2.63) 



We treat the Coulomb interaction of the Hubbard model (2.53) by replacing it with 



where the operator field d" satisfies the constraint equation 



(2.64) 



(2.65) 
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and ((i") = 1 at every site j. In the low-energy, small-momentum limit, the modulus 
Mj can be taken a s a c onstant M. The factor (—1)-', alternating sign from site to 
site, is included in ( 2.64 ) because we are interested in the antiferromagnetic properties 
of the model. In the mean-field approximation, the operator is considered to be 
a classical field. It then describes the staggered magnetization and thus becomes the 
order parameter of the Neel state. In this approximation, the continuum limit of the 
relevant term of the Coulomb interaction becomes 



(2.66) 



where we introduced the effecti ve co uplin g cons tant T, — — | The term (|2.66|) to 



gether with the two other terms ( 2.59| ) and ( 2.62 ) describe the Neel state of the extended 
Hubbard model in the mean-field approximation. The corresponding Lagrangian reads 



where the Fermi velocity up is set to unity, and 



(2.67) 



(2.68) 



are the Dirac matrices. Moreover, ^ = 9^7/^- Although mean-field theory should in 
general not be trusted in lower dimensions, it in this case seems to capture the essential 
physics. 

The effective theory we are seeking is obtained by integrating out the electron 
fields. (We tacitly switched here from the operator formalism to the functional-integral 
approach.) In our discussion, A is considered to be a constant. If one is interested in 
the dynamics of the charge density wave, which we are not, A is assumed to have a 
spacetime-dependent phase, A(a;) = A exp[«(/3(a;)]. 

Since the Lagrangian (2.67) is quadratic in the electron fields, these are readily 
integrated out. The ensuing effective action S'ofr reads 



S^ff = -iTr h-i{^ - A - S d"cr"), 



(2.69) 



which we again evaluate in a derivative expansion [^. On our way, we encounter 
integrals which diverge in the ultraviolet. To handle these, we use dimensional regular- 
ization and generalize the integrals to arbitrary spacetime dimensions D. 

The lowest-order term is of the form r{d")'^, where r is a constant that tends to 
infinity when the dimensional-regularization parameter e = 1 — is taken to zero. 



This term merely renormalizes the first term at the right-hand side of ( 2.64 ). Neglecting 



for the mom ent a possible topological term, we obtain as next order in derivatives the 
contribution [ 111 \ 

Cos = ^ {d^d'^f , (2.70) 

where l/g^ = S^/GttA^. This is a kinetic term for the staggered magnetization in- 
duced by quantum effects. Note that l/g^ > 0, which is required for stability. The 
0(3) nonlinear sigma model (2.70) in 1h-1 spacetime dimensions is known to have only 
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one phase, viz. the so-called quantum disordered phase, with a finite correlation length 
[112]. Due to strong infrared interactions of the "Goldstone modes", they acquire an 



energy gap for all values of the coupling constant g. For this reason one frequently 



refers to the (iH-l)-dimensional state described by the effective theory (2.70), with the 
constraint [dP)^ = 1, as a short-rangedNeel state. In higher spacetime dimensions the 
model has a phase transition between an ordered, weak-coupling phase {g < g^r) and 
a disordered, s trong -coupling {g > ga) phase, where g^ is the critical coupling. 

The form ( 2.69 ) of the action is not conveni ent to derive the topological term. To 
obtain this term we instead follow Jaroszewics [113| and introduce, as we have done 



in the previous section, the decomposition i}} = sx, with s{x) the spacetime-dependent 
SU(2) spin rotation matrix defined in (2. 1 ). The effect of this unitary transformation is 
to rotate the staggered magnetization in every point of spa cetime into a fixed direction. 
With the help of this decomposition, the Lagrangian ( 2.67 ) can be cast in the equivalent 
form 

C = x{^~^-^-^a^ -$)X. (2.71) 

where = —is'^d^s = B'j^a" is the 2x2 matrix-valued field we encountered before 
in the context of superfluid "^He-a. As was the case there, this Lagrangian is invariant 
under the local spin gauge transformation 



^(a;) ^ e— f-)'^' x(x), Blix) ^ B^x) + d^a{x), 



(2.72) 



with spin gauge field B^. This local gauge freedom derives from the circumstance that 



the transformed rotation matrix s' ~ s exp{iaa^) also satisfies (2.1), i.e., s'V"(i"s 



a . That is to say, determines the spin rotation matrix s only up to a multiplicative 
U(l) factor exp(ia(T^). In terms of _B^, the kinetic term {d^id")^ is given by 



(2.73) 



where the summation over the spin index only involves the components a = 1, 2. The 
kinetic term is thus seen to be independent of the spin gauge field B^^. The term (Sj^) , 
which would give the gauge field a mass, is not generated as it would violate local spin 
gauge symmetry. On the other hand, the topological term that in principle can arise in 
this context depends only on SjJ 



2tt 



(2.74) 



It contains, as is frequently the case with topological terms, the antisymmetric Levi- 
Civita symbol and it is linear in derivatives. It can alternatively be written in terms of 
the 0(3) nonlinear sigma field as [82|: 



Ce = —e^.e'^l'^d'-d^dPd,d\ 
an 

where we recognize the winding number of the map 

d^ix) : -> 



(2.75) 



(2.76) 
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of compactified spacetime into the internal two-sphere parameterized by d", cf. 



(2.26) 



Since it is independent of B}^ and i?^, these components may, as far as deriving the 



topological term is concerned, be neglected in the Lagrangian (2.71). This leads us to 
consider the following Lagrangian 



4- (A + E) 



Xe 



(2.77) 



where we made explicit the two fermion species Xc and Xo, which are seen to have 
different masses and to couple to the local spin gauge field in an opposite way. This 
is because the Neel state is characterized by a staggered magnetization, flipping sign 
when we go from an even to an odd site. 



To investigate whether the topological term (2.74) is induced by the fermions, we 
could proceed as before and apply the derivative expansion. But in 1+1 dimensions we 
have another tool available, namely that of bosonization, which we shall use instead. 
Since (2.77) contains only a U(l) symmetry, we can employ Abelian bosonization. One 
could also bosonize the original theory (2.67), as was first done by Wen and Zee [pT]], 
but this involves the more intricate non- Abelian bosonization. The Abelian bosoniza- 
tion rules — on which we shall comment in the next chapter — read 



XcXc 

Xo7a»Xo 



K co8{v' 47r(/)c) 



(2.78) 



with K a for our purposes irrelevant positive constant and 4>c a real Bose field. With 
these and similar rules for Xo^ the Lagrangian (2.77) can be represented in bosonized 
form as 



where we introduced the fields 



1 



V, 



(2.79) 



(2.80) 



0± = (0e ± </'o)/V2, 

and where the potential V is given by 

V = 2AA'cos(\/2^(?!)+) cos(V2^0-) - 2T.K sm{V2^(t)+) sin(%/2^0_). (2.81) 



Note that in (2.79 ) only the field 4>- couples to the local spin gauge field. The potential 
V has a minimum 

_ r -2|A|if for |A|>|I]| 
1/min - I _2|S|X for |A| < |S| 



at 



for A > S > 0; 



(2.83) 
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4VI f-0<A<S, 



(2.84) 



and similar values for A and S negative. With these constant (/)-values, the third term 



in the Lagrangian ( 2.79 ) takes the form of t he top ological term ( 2.74 ). We find that the 
^-parameter, which is an angle, is given by [1 11 1 



= 7r[sgn(A + S)-sgn(A-S)], 



(2.85) 



valid for arbitrary sign of A and E. 

More specific, 9 = 0(mod27r) when |A| > and 9 — 7r(mod27r) when 
|A| < That is, for |A| > the effective theory is simply the 0(3) nonlinear 
sigma model (2.70), with coupling constant squared > Qn. As we remarked before, 
the state it describes is the short-ranged Neel state with finite correlation length. 

For I A| < |S]|, on the other hand, (2.70) is to be augmented by the topological term 
In terms of the field d° 



(2.74) 



the effective theory now reads 



8tt 



(2.86) 



with = ±7r. The 6'-term changes the physical content of the 0(3) nonlinear sigma 
model in a dramatic way. This term has been extensively studied in the context of the 
antiferromagnetic Heisenberg spin chain which can also be mapped onto the nonlinear 



sigma model [114|. Whereas the half-integer spin chain is known to be gapless, the in- 
teger spin chain has a gap in the excitation spectrum — the so-called Haldane gap. The 
two cases differ in the topological term which is induced with a coefficient 6 — 2-Ka 
(mod27r), where a is the spin. For integer spins the coefficient is zero and the ex- 
citation spectrum has an energy gap. For half-integer spins, on the other hand, the 
coefficie nt is nonzero and it can be shown that as a consequence the spectrum becomes 
gapless 1 115|. Using these results, we obtain the phase diagram depicted in Fig. 2.3 



The standard Hubb ard chain is recovered by setting A = 0. It is therefore effectively 
described by ( 2.70 ), implying that the model has a gapless spectrum and thus metalhc 
properties for all values of the Coulomb interaction U . The disordered phase of the ef- 
fective theory, describing the insulating state of the extended Hubbard model, becomes 
accessible when the Peierls energy gap A is taken sufficiently large compared to S, 
measuring the strength of the Coulomb interaction. 

The fermions of the original theory appear in the bosonized theory as solitons to 
which we next turn. An example of a soliton in the regime < A < S, where the 
groundstates are given by ( 2.84 ), is the configuration 0+(a;i) = = f{xi), 

where the function /(xi) interpolates between the two different ground states 



(2.87) 



This soliton carries one unit of fe rmio n-number charge N . This follows because ac- 
cording to the bosonization rules ( 2.78 ), the fermion-number current is represented 

as 



di,<j) 



(2.88) 
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Figure 2.3: Phase diagram of the extended Hubbard model (2.67), where E represents 
the on-site Coulomb interaction and A is the Peierls energy gap. 



implying that 



N 



jo = -\/ - [0+(oo) - (/)+(-oo)] = 1. 



(2.89) 



Similar solitons of unit fermion-number charge can be constructed in the other sectors 
of the phase diagram. 

It can be easily demonstrated that these solitons carry in addition a spin i, so that 
they have the same quantum numbers as the original fermions [^. Using arguments 
similar to those of the previous section, one can show that the spin current, given by 



3u = ^5^"7mV', 



(2.90) 



is again represented in the effective theory by the isospin current of the 0(3) non- 
linear sigma model: 

ja fOiPiif) r2 91) 



Also Eq. ( 2.42 ) is valid, so that 



I dC 1 /2 „ , 
2aM = -2V^^^^- 



This gives as spin for the soliton under consideration 



1 



(2.92) 



(2.93) 



It may therefore be thought of as representing the fundamental fermion of the origi- 
nal theory. Similar conclusions hold for the solitons in the other sectors of the phase 
diagram. 



(glAmstex 



2.2 Peierls Instability 



63 



Up to this point we have only considered the part of the bosonized theory involving 
the scalar fields <j)± . We now turn to the nonlinear sigma field d" and will argue that 
twists in this field may give certain solitons an anomalous spin. For definiteness we 
shall assume in the following that the space dimension is compactified into a circle 



-L < xi < L, with xi = —L and xi = L labeling the same point. The model (2.67) 



admits two possible boundary conditions for d", namely periodic d°'{—L) — d°'{L) 



and antiperiodic d'^{—L) = — ones |81 1. That is, the two phases split up into an 
even and odd sector. 

We are interested in the contribution jg ^ to the spin current stemming from the 
topological term 

The corresponding charge Sg is zero in the even sector, and 

Se = / Je.o = f - d'^i-L)] = f ^"(-i), (2.95) 

Jxi 47r In 

in the odd sector of the theory. It follows that Sg is nonzero only in the gapless phase 
where 9 = ±7r, and with antiperiodic boundary conditions imposed on d". A soli ton in 
this part of the theory acquires consequently an extra contribution ag to the spin given 
by 

cjg=±\. (2.96) 
Since the solitons have a canonical spin of i, the topological term in the odd sector 



of the gapless phase transmutes these fermions into bosons |81 1. This ability of Cg to 
change statistics it shares with the Chern-Simons term in (2+l)-dimensional theories. 
We next consider the critical lines A = ±E. At these lines, either the component 



Xo or Xo in (2.77) is gapless, and the above analysis does not apply. Let us first focus 



on the critical line A = E. The bosonized form of (2.77) then reads 



\{dt.^^f + \{^^.4>of - (A + E)Xcos(V4^</)e) 

\e^Mo-^o)df,Bl, (2.97) 

where the gaplessness of the field Xo is reflected in the absence of a cosine term for 
0o. As a result, the integration over this field is a simple Gaussian which can easily be 
carried out to give a gauge-invariant mass term for the local spin gauge field 




dadu 



Bt (2.98) 



Q2 

This mass generation by gapless fermions in two spacetime dimensions is the famous 



Sch winger mechanism [116|. The contribution of Xc to the effective theory may be 



approximated by substituting the value 0c = 0, or ^/tt /2 in (2.97), corresponding to 



the minimum of the potential (A -I- Y,)K cos(v47r(/ic) for A = E < 0, or A = S > 0. 
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It follows that only in the latter case a topological term with 6 = tt is generated. 
The Schwinger mechanism operates similarly at the critical line A — — S. Here, the 
topological term (with 9 — n) persists only at the lower half of the Une, defined by 

A = -S > 0. 



2.3 Statistics-Changing Phase Transition 



The (l+l)-dimensional model ( ^.71 ) w e dis c ussed in the previous section has similar 



characteristics in 2+1 dimensions |117 118, 119, 12C]. While the Dirac Hamiltonian 



naturally appears in (l+l)-dimensional systems where it describes the fermionic exci- 
tations around the two Fermi points [|2|, this is not the case in 2+1 dimensions. How- 
ever, there exist two- dimensional systems — so-called semimetals — that have point-like 
Fermi surfaces [121 ]. In these materials, the valence and conduction band intersect in 
discrete points. Near such a degeneracy of two energy levels, the Hamiltonian describ- 
ing the full theory may be approximated by a 2 x 2 Dirac Ha miltonian descri bing just 
the two-level subsystem in the vicinity of the diabolic points [122, 123, 124]. In this 



section, we briefly discuss some of the salient features of the (2+l)-dimensional model. 

For our purposes, it suffice to consider the simplified, Abelian version of the model 
given by ( 2.77 ). We have chosen a two-dimensional representation of the Dirac algebra 



7o = /3 = T3, 7i=*T2, 72 = iTi, (2.99) 

{iti, li^} = 25^.., g^Lu = diag(l, -1, -1), (2.100) 
with T the Pauli matrices. With the help of the derivative expansi on [p6| ], one can easily 



compute the induced fermion-number current . It is given by 1 1 1 9 1 



1 - ^G 



where is the dual field — €i_i^\d^B'^ and 



,?=i[sgn(A + I])-sgn(A-S)] 



(2.101) 



(2.102) 



The first term here is the contribution from the even-site spinor Xc, while the second 
term is the contribution from the odd-site spinor Xo- The relative minus sign reflects 
the opposite charge of Xc and Xo with respect to the local spin gauge symmet ry. W e 
encountered the dual field G^ previously in the context of superfluid ■^He-a [see ( 2.29 )] , 
where it was shown that it is proportional to the topological current of the 0(3) 
nonlinear sigma model. The corresponding charge Q is the winding number of the 
map (2.26). 

We thus arrive at the conclusion that the induced fermion-number current is 
proportional to the topological current of the effective theory. 



(2.103) 



The proportionality constant d is nonzero for | E | > | A | , where sgn( A + S) — sgn( A - 
E) = 2sgn(S), and zero for |E| < |A|. 
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Let us first discuss the case |E| > |A|. The situation of two different conserved 
currents which become proportional to each other is typical for ordered states where 
the residual symmetry links up different groups of the original symmetry. For the 



model under consideration, ( ^.103| ) implies that the fermion-number symmetry U(l)^ 
generated by the charge N = J^ja, and the flux symmetry U(l)* generated by the flux 



Go, 



(2.104) 



are spontaneously broken in the following manner: 

U(l)^ xU(l)* dU(1)^-^*. 



(2.105) 



To identify the Abelian Goldstone mode associated with this spontaneous symmetry 
breakdown, we also calculate the induced "Maxwell" term. 



4 ' 



(2.106) 



with Gfj_i, = d gBf — d^B^. To the one-loop order we find for the coefficient 11 
appearing here 1 1 19 1 



n 



1 

8^ 



lA-Sl 



(2.107) 



If we imagine integrating over the local spin gauge field, the partition function can be 
written as 



DS^exp 



nc 



(2.108) 



where we omitted a gauge-fixing factor for the spin gauge field and coupled the current 
to a background field A^. This coupling allows us, by diff erentiating with respect 
to Afj_, to compute the induced fermion-number current ( 2.103 ). 

To make the Goldstone mode explicit we employ a procedure frequently used to 
derive dual theories and introduce the change of variables B^ in the functional 

integral. Since the dual field fulfills the Bianchi identity 9^(5^ = 0, we introduce 
a Lagrange multiplier ip in the functional integral to impose this constraint: 



ipd^G^ 



(2.109) 



Performing the Gaussian integral over the dual field G^i, we obtain an expression for 
the partition function in terms of a gapless scalar field: 



Z 



exp 



2n 



(9,,^ + -dA^ 



(2.110) 



This is the Goldstone field we were seeking. It should be noted that in 2h-1 dimensions, 
both a massless vector field and a real scalar field represent one degree of freedom. 
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This can be easily understood by noting that a photon has only one transverse degree 
of freedom in two space dimensions. In terms of the Goldstone field Lp, the fermion- 
number current becomes 

]^.^~^0^'P + 'dA^) (2.111) 



as in BCS theory. From this we conclude that the (2+l)-dimensional model (2.71) 
exhibits superconductivity. It should be noted that the mechanism leading to this, with 
the spin gauge field playing a decisive role, is entirely different from that in classic 
superconductors. 

On comparing the expression ( [2.1 1 1[ ) for the fermion-number current with the one 



obtained directly from (2.108), we conclude that the dual field is related to the 



Goldstone field via 

G^ = -^d^4>. (2.112) 

We next turn to the case jSj < |A|. When one of the two spinors Xc or Xo be- 
come gapless, i.e., A + S = OorA — E^O, the induced fermion-number current 



discontinuously drops to zero, as was first noticed by Chen and Wilczek 1 125 1. So, for 
|S| < I A| the argument that the fermion-number current and the topological current 
are proportional, which led to the conclusion that we have the spontaneous symmetry 



breaking (2.105), is invalid and the symmetries U(l)^ and U(l)* remain unbroken 



here. The discontinuities in the induced charges at the lines |S| = |A| merely reflect 



that these are critical lines [1 19| 



The restoration of symmetry when one passes from the broken phase to the sym- 
metric phase should be accompanied by a loss of the gapless mode. To see that this is 
indeed the case we note that besides the Maxwell term, also a Chern-Simons terms is 
generated here, 

Ce ^ ^e^^^xBld.Bl, (2.113) 

with 9 = — i7r[sgn(A + E) + sgn(A — E)]. The first and second term here arise from 
the even-site and odd-site spinors, Xc and Xo, respectively. Observe that in the ordered, 
gapless phase where |E| > |A|, the Chern-Simons term is zero. The Euler-Lagrange 
equation for in the presence of a Chern-Simons term and with the background field 
set to zero yields d^G^y oc G^, or in dual form d^Gi, ~ di,G^ cx G^^. Upon 
taking the divergence of this last equation and using the previous one as well as the 



Bianchi identity d^G^ = 0, one finds 1 126] 



a^G,, cxG^. (2.114) 



This equation together with (2.112) shows that Lp has become massive. From this we 
conclude that due to the presence of the Chern-Simons term in the symmetric phase, 
the gapless mode of the ordered phase vanishes, as it should be. 

A last point we like to mention is that the phase transition we have been discussing 
here is a statistics-changing phase transition. This is easily understood using the result 
[ 9l|] th at a Chern-Simons term imparts spin to a skyrmion. The presence of this term 



(2.113) in the unbroken phase turns a skyrmion, which is an ordinary boson in the 



broken phase where no Chern-Simons term is generated, into a fermion. Whence, on 
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crossing the phase boundary |E| = |A|, skyrmions undergo a spin transmutation and 
change their statistics from fermionic to bosonic or vice versa. 



2.4 Fluxons 

In this section, we examine a gas of electrons confined to a plane. We present a simple 
method to calculate the quantum numbers induced by a uniform magnetic background 



field [123, 127|. We will in particular focus on the quantum numbers carried by a 
fluxon — a point-like object carrying one unit of magnetic flux that can be pictured as 
being obtained by piercing the spatial plane with a magnetic vortex. It is found that 
it carries both the charge and spin of a fermion. Unlike what one expects, the spin of 
the fluxon does not originate from a Chern-Simons term, but from a so-called mixed 
Chern-Simons term involving two different gauge fields, viz. the electromagnetic gauge 
field and spin gauge field which describes the spin degrees of freedom. 
Let us consider the Lagrangian 

C = + ^^- Hp)i^ + hi:'<\a'^i:, V = J| ) , (2.115) 

governing the dynamics of the Pauli spinor field il>, with Grassmann components ^/j^ 
and ill describing the electrons of spin-| and | and chemical potential /i « ep = 
kp/2m which is well approximated by the Fermi energy. We introduced an external 
source b coupled to the spin density, to compute the induced spin. The Pauli Hamilto- 
nian 

= :^ip - eA)^ - gofiB^a^H + eAo, (2.116) 
Zm 

with fiB — e/2m the Bohr magneton and go = 2 the free electron ^-factor, contains 
a Zeeman term which couples the electron spins to the background magnetic field H. 
Usually this term is omitted. The reason is that in realistic systems the ^-factor is much 
larger two — the value for a free electron. In strong magnetic fields relevant to, say, 
the QHE, the energy levels of spin- J, electrons are too high to be occupied so that the 
system is spin polarized, and the electron's spin is irrelevant to the problem. Since we 
consider a free electron gas, we include the Zeeman term. We describe the uniform 
magnetic field H by the vector potential Aq = Ai ~ 0; A2 = Hxi. The energy 
eigenvalues of the Pauli Hamiltonian are the famous Landau levels 

ij±^M(;+i)_£^,^, (2.117) 

with (7± = ±i for spin-t and spin-| electrons, respectively. 

Integrating out the fermionic degrees of freedom, one finds the one-loop effective 
action: 

S^a^ [ C,fi^-iTT\n{po-Hp+fi+^b<j^). (2.118) 

J X 

The key observation to evaluate this expression is that with our gauge choice, the theory 
is translational invariant. Although a translation by a distance Ai in the a; 1 -direction 
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changes the vector potential, A2 ^ A2 + HAi, this change can be canceled by a 
gauge transformation Af^ + d^a, with gauge parameter a = HA1X2, so that 

the theory is invariant under the combined symmetry. As a result, the theory effectively 
reduces to a (0+l)-dimensional theory, and each Landau level is infinitely degenerate. 
The number of degenerate states per unit area is given by \eH\/2Tr for each level. With 
these observations, we can write Ccs as 

^off = I (^0 - + M + 

+ ln(A:o -£^;,- 5^)]- (2.119) 

The presence of only an integral over the loop energy and no momentum integrals 
reflects that this is effectively a theory in zero space dimensions. 

The induced fermion number is obtained by differentiating £cff with respect to the 



chemical potential, as can be inferred from the original Lagrangian (2.1 15), and setting 
the external source h to zero, 

dLcS 

Jo 



(2.120) 

fc=0 



d^JL 

To evaluate the resulting energy integral, we employ the integral 

' -0{O, (2.121) 

ZTTl fco + ? — 

containing, as usual in nonrelativistic calculations [Q, an additional convergence fac- 
tor exp{ikori). The function 9{(,) at the right-hand side is the Heaviside unit step func- 
tion. The value for the induced fermion-number density thus obtained is 

io = ^a++L), (2.122) 
where l± is the number of filled Landau levels for spin-f and spin- J, electrons, 

1^ = e(^ + 1), (2.123) 



with 5(x) the integer-part function denoting the largest integer less than x, and 

eH 

l^±=t-i+ cr± (2.124) 

m 

their effective chemical potentials. Implicit in this framework is the assumption that the 
Fermi energies of the spin-f and the spin-| electrons lie between two Landau levels. 



so that the integer-part function is well defined. We see from ( 2.122 ) that the so-called 
filling factor i/j^, defined as 

.. = p|^ = ^. + /-, (2.125) 

takes on integer values only. This was to be expected for an ideal electron gas at zero 
temperature; given a value of the Fermi energy ep, a Landau level below the Fermi 
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surface is filled, while a level above it is empty. (When the Fermi energy and that of a 
Landau level coincide, the value of the integer-part function S is ambiguous.) 

If, in addition to a magnetic field, there is also an electric field E in the plane, a 
Hall current is induced perpendicular to the two fields. The current carried by the filled 



Landau levels is obtained by multiplying the induced density (2.122) with the drift 
velocity E/H. In this way one finds: 



J2 



sgn{eH) 
2^ 



(2.126) 



where, without loss of generality, the electric fi eld is c hosen in the xi direction: Aq = 
—Exi. The induced fermion-number density ( 2.122 ) and Hall current ( 2.126 ) corre- 
spond to a Chern-Simons term. 



in the effective theory, with 



sgn(eH) 
27r 



{1+ + 1-). 



(2.127) 



(2.128) 



Because of the factor sgn(eH), which changes sign under a parity transformation, this 
Chern-Simons term is invariant under such transformations. 

We next turn to the magnetic properties of the two-dimensional electron gas in a 
constant magnetic field. The induced spin density calculated from the effective La- 
grangian via 



db 



b=0 



(2.129) 



(2.130) 



turns out to be independent of the filled Landau levels l±, viz. 

eH 
47r 

This follows from the symmetry in the spectrum i?;^^ ~ and E^^ = 

E^^ if eH < 0. The spin magnetic moment, or magnetization M, is obtained from 
( 2.130 ) by multiplying s with twice the Bohr magneton ^b, 



-H. 



(2.131) 



(2.132) 



47rm 

This leads to the text-book result for the magnetic spin susceptibility xp 

_ dM _ e2 
dH 47rm 

with v{Q) = m/2-K the density of states per spin degree of freedom in two space 
dimensions. To see how the spin contribution ( ^.132 ) to the magnetic susceptibility 
compares to the orbital contribution, we evaluate the feg -integral in the Lagrangian 
(|2.119[) with & = to obtain 



\eH\ 
27r 



i=0 <;=± 



E^i^-Ef) 



(2.133) 
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The summation over the Landau levels I is easily carried out with the result for small 
fields 



C. 



cff 



An ^ 

<;=± 



{eHf 
Am 



[eHl 



[(2af -1], 



(2.134) 



where o- = and /i± is given by (2.124). The first term at the right-hand side of 
(2.134), which is independent of the magnetic field, is minus the energy density of a 
free electron gas particle contribution. 



2tt 



2m 



2m 



(2.135) 



The factor 2 at the right-hand side accounts for the spin-t and spin-| fermions, while 
the step function shows that only energy levels below the Fermi level contribute. The 
energy density is negative because the levels are measured relative to the Fermi energy. 



The second term in (2.134) yields the low-field susceptibility 



X= (-1)2^+12^*1 K0)[(2a)2-1)], (d = 2). (2.136) 
We have cast it in a form valid for a = 0,^, 1. The term involving the factor (2cr)2 is 



the spin contribution which reduces to (2.132) for a 
the three-dimensional expression 



i . Let us compare (2.136) with 



(2.137) 



where i/(0) = mfcF/27r^ now denotes the three-dimensional density of states per spin 
degree of freedom at the Fermi sphere. We see that the ratio of orbital to spin contri- 
bution to X is different in the two cases. In addition, whereas a three-dimensional elec- 
tron gas is paramagnetic (x > 0) because of the dominance of the spin contribution, 
the two-dimensional gas is not magnetic (x = 0) at small fields since the diamagnetic 
orbital and paramagnetic spin contributions to x cancel. 

The induced currents we just calculated may be used to compute the induced quan- 
tum numbers carried by certain scalar field configurations. The case of interest to us is 
the fluxon. This is a point-like object carrying one unit of magnetic flux 27r/e which 
can be described by a magnetic field 



= — <5(x). 

e 



(2.138) 



According to (2.130), a single fluxon carries a spin a-(^ = h and, since for small fields 



Jo 



TT 



\eHl 
2tt 



(2.139) 



where the first term at the right-hand side is the two-dimensional fermion-number den- 
sity fim/ir = kp/2TT, also one unit of fermion-number charge. That is, in a nonrel- 
ativistic electron gas, a fluxon carries the quantum numbers of a fermion. However, 
the close connection between the spin of a fluxon and induced Chern-Simons term for 
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arbitrary large fields that exists in a relativistic context, is lost. This can be traced back 
to the fact that in nonrelativistic theories, the electron spin is an independent degree 
of freedom. Below, we point out that the spin of the fluxon does not derive from the 
ordinary Chern-Simons term, but from a so-called mixed Chern-Simons term. Such a 
term is absent in a relativistic context. 

It is interesting to note that whereas the spin of an antifluxon, which can be de- 
scribed by the magnetic field 



He. 



2-K 



<5(x), 



(2.140) 



is — its fermion-number charge is the same as that for a fluxon because of the abso- 



lute values appearing in (2.139). This is what one expects in a nonrelativistic theory. 

Because fluxons carry the quantum numbers of fermions, the exclusion principle 
forbids two fluxons to be in the same state. This is important when calculating, for 
example, the orbital angular momentum L of a state with A^^ fluxons since they have 



to be put in successive orbital angular momentum states. We find for L \ 128] 



(2.141) 



In this way, the total angular momentum J = S + L, with S = crg,iVg, the total spin 
carried by the fluxons, becomes 



J = S + L = a^Nl^^Nl. 



(2.142) 



We next investigate the origin of the induced spin density ( 2.130 ) we f ound i n the 
nonrelativistic electron gas. To this end we slightly generalize the theory ( Z.115 ) and 
consider the Lagrangian 



Po 



eAn 



1 

2™ 



(P - eA)2 



(2.143) 



It differs from ( 2.1 15 ) in that the spin source term is omitted and the magnetic field in 
the Zeeman term is allowed to point in any direction in spin space labeled by the index 
a — 1,2, 3. It is convenient to consider a magnetic field with fixed magnitude, but a 
spacetime-dependent orientation 



H"{x) ^ Hd"{x), 



(2.144) 



where H is constant, and d" i s a uni t vector in spin space. The gauge field A^^ appearing 
in the extended Lagrangian ( 2.143 ) still gives the magnetic field perpendicular to the 
plane V x A = iJ. As we did various times before, we make the decomposition 



tjjix) = s{x)x{x) 



(2.145) 



with s{x) the spacetime-dependent SU(2) matrix first introduced in (2.1). In terms of 
these new variables the Lagrangian (2.143) becomes 



Po - eAo - Bq + fi 



1 
2m 



(p - eA - B) 



X 



eH 
2m 



xV^x, (2.146) 



(glAmstex 



2.4 Fluxons 



72 



where — —is^d^s is the field (2.4). The theory takes again the form of a gauge 
theory with gauge field i?^. The spin-density operator 



becomes in these new variables [see (2.42)] 



^0 =i?"Via^X = -^i?"^ — 



aS^ 



(2.147) 



(2.148) 



In deriving the first equation we employed the identity (2.40), relating the SU(2) ma- 
trices in the j = ^ representation, 



s = exp {i^e^a") , 
to those in the adjoint representation (j = 1), 

R = exp J,",j) . 
The matrix elements of the generators in the latter representation are 



.07 



(2.149) 



(2.150) 



(2.151) 



The projection of the spin density Jq onto the spin-quantization axis, i.e., the di- 
rection of the applied magnetic field [see Eq. (2.5C)], 



1 dC 



2dBl' 



(2.152) 



only involves the spin gauge field B^^ . So, when calculating the induced spin density 
we may set the fields B}^ and i?^ to zero and consider the simpler theory 



2m 



(2.153) 



where the effective Fermi energies for the spin-f and spin- J, electrons are given in 
(2.124) and eA^ — eA^ ± BjJ. Both components X] ^nd xi induce a Chern-Simons 
term, so that in total we have 



\{e+ + e^)€,,,^{e^A,,d,Ax + BlB^Bl) 



(2.154) 



)^nu\B^^di,A\, 



where the last term involving two different vector potentials is a mixed Chern-Simons 
term. The coefficients are given by, cf. (2.128) 



^ sgn(eg) 



(2.155) 
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assuming that \eH\ > ^\eijdiBj\, so that the sign of eH is not changed by spin gauge 
contributions. The integers l± are the number of filled Landau levels for spin-f and 



spin-| electrons given by (2.123). Since — = sgn(ei/), we obtain for the induced 
spin density s precisely the result (2.130) we found in the preceding section. 



2 dB^ 



_ eH 
B?=o An ' 



(2.156) 



The present derivation shows that the induced spin i n the n onrelativistic electron gas 
originates not from the standard Chern-Simons term (2.127), but from a mixed Chern- 
Simons term involving the electromagnetic and spin gauge potential. 

The first term in (2.154) is the standard Chern-Simons term, the combination 6+ + 



9- precisely reproduces the result (2.128) and is related to the induced fermion-number 
density (2.122). 
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Chapter 3 

Dual Theories 



In this chapter we shall discuss dual theories of a superfluid "^He film and of a BCS 
superconductor in two and three space dimensions. The first use of a duality transfor- 
mation in contemporary physics is generally attributed to Kramers and Wannier who 
applied such a transformation to the Ising model on a square lattice [1291. The models 
we consider in this chapter all possess vortex solutions. The stability of these defects, 
which appear in the original formulations as singular objects, is guaranteed by a non- 
trivial topology. The surrounding hyper sphere of a vortex is a circle [r = 1 in Eq. 
(2.25)], so that the relevant homotopy group is the fundamental homotopy group. In 
d = 2, a vortex is a point object (e^ = 0), while in d = 3 it is a line defect (ed — 1). A 
loop circling any of these vortices cannot be deformed to a point without encountering 
a singularity. 

The duality transformations of the models discussed here are aimed at obtaining a 
field theoretic description of the vortices they contain. Whereas these excitations are of 
topological nature in the original formulations, they become the elementary excitations 
of the dual theory. 



3.1 Superfluid ^He Film 

In this section, we shall derive the dual theory describing a ''He film. It is well known 
that this system undergoes a Kosterlitz-Thouless phase transition at a temperature well 
below the bulk transition temperature. The superfluid low-temperature state is char- 
acterized by tightly bound vortex-antivortex pairs which at the Kosterlitz-Thouless 
temperature unbind and thereby disorder the superfluid state. The disordered state, 
at temperatures still below the bulk transition temperature, consists of a plasma of un- 
bound vortices. We shall see in what way the dual theory, which gives a field theoretic 
description of the vortices, accounts for these phenomena. 

The phase transition is an equilibrium transition, we can accordingly ignore any 
time dependence. The important fluctuations here, at temperatures below the bulk 
transition temperature, are phase fluctuations so that we can consider the London limit 
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= 




= 27r 



Figure 3.1: Vortex of unit winding number The dashed line denotes the cut in the 
spatial plane along which ip makes a jump of 2tt. 



and take as Hamiltonian 

H = ^Ps^l (3.1) 

where ps is the superfluid mass density which we assume to be constant and Vg is the 
superfluid velocity 

V, = -(V(^- yP). (3.2) 

TO 

Following Kleinert ^1%, we included a vortex gauge field (p^ to account for possible 
vortices in the system. A vortex in two space dimensions is, as we mentioned in the 
introduction, a point-Uke object. It is characterized by the winding number of the map 

(/.(x) : Si ^ Si (3.3) 

of a circle S^. around the vortex into the internal circle parameterized by (p. When 



we circle an elementary vortex once, (p changes by 2n (see Fig. 3.1). In itself this 
is harmless because (p is a compact field with periodicity 2tt, i.e., Lp and ip + 2tt are 
identified. However, the jumps taking place along a cut in the spatial plane would turn 
the superfluid velocity into a discontinuous field. This is physically unacceptable. The 
field ip^ is introduced with the purpose to compensate for these jumps. Its curl yields 
the vortex density riv consisting of delta functions at the vortex positions x" 

V X v?P(x) = -27rnv(x) = -27r^ w;"J(x - x"), (3.4) 
leading to the vorticity 

V X Vs = 27r— . (3.5) 

TO 

The integer Wa is the winding number of the vortex located at x" . We shall restrict our- 
selves to vortices of unit winding number, so that Wa — ±1 for a vortex and anti vortex, 
respectively. 

The canonical partition function describing the equilibrium configuration of 7V+ 
vortices and iV_ antivortices in a superfluid ^He film is given by 

^--iN^n£/D,exp(-/3/^H), (3.6) 



with Ti. the Hamiltonian (3.1) and N = iV+ + the total number of vortices and 



antivortices. The factors N+l and A^_! arise because the vortices and antivortices are 
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indistinguishable, and Yla /x° denotes the integration over the positions of the vortices. 
The functional integral over ip is Gaussian and therefore easily carried out, with the 
result 



1 



N+\N-\ 



n 



exp 



^ WaWp In 



a, 13 



(3.7) 



The constant a, with the dimension of a length, in the argument of the logarithm is 
included for dimensional reasons. Physically, it represents the vortex core diameter. 



Apart from an irrelevant normalization factor, Eq. (3.7) is the canonical partition func- 

±q, where 



tion of a two-dimensional Coulomb gas with charges qa 



qwa 



q = y/2TTps/T 



(3.8) 



Let us rewrite the sum in the exponent appearing in ( |3.7| ) as 



a,/3 



a, 13 



In 



- ln(0) 



+ ln(0) {T.1c}j 



(3.9) 



where we isolated the self-interaction in the last term at the right-hand side. Since 
ln(0) — —oo, the charges must add up to zero so as to obtain a nonzero partition 
function. From now on we will therefore assume overall charge neutrality, 9" — 0, 
so that iV+ = iV_ = N/2, where N must be an even integer. To regularize the 
remaining divergence, we replace ln(0) with an undetermined, negative constant — c. 
The exponent of (3.7) thus becomes 



2^' 



a, 13 



|x" -x/^l 



2 ^ 



In 



(3e,N, (3.10) 



where = cq^ /2 physically represents the core energy, i.e., the energy required to 
create a single vortex. In deriving this we used the identity X^Q^t/j 9q9/3 = ^ X]q 5'q ^ 
—Nq'^ which follows from charge neutrality. Having dealt with the self-interaction, we 



limit the integrations Yla /x^ (3.7) over the location of the vortices to those regions 



where they are more than a distance a apart, |x" — x*^^! > a. The grand-canonical 
partition function of the system can now be cast in the form 



E 

N=Q 



-,N 



[{N/2) 



112 



n 



exp 



^ <Z„<Z0 In (Ix^-x'^' I) 



a^l3 



(3.11) 



where z = exp(— /3ec) is the fugacity. We suppressed an irrelevant dimensionful factor 
a^iPi xhe sys tem is kn own to undergo a phase transition at the Kosterlitz- 

Thouless temperature [130, 131 1 



1 2 7^ 

4 2 



(3.12) 
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triggered by the unbinding of vortex-antivortex pairs. It follows from this equation 
that the two-dimensional superfluid mass density Ps{T), which varies from sample 
to sample, terminates on a line with universal slope as T approaches the Kosterlitz- 
Thouless temperature from below [132]. 

To derive the dual theory we note that ln(|x|) is the inverse of the Laplace operator 



i-V2ln(|x|)=5(x). 



(3.13) 



This allows us to represent the exponential function in (3.11) as a functional integral 
over an auxiliary field 0: 



exp 



^g„g/3ln(|x«-x^|) 



D0exp 



(3.14) 



where p(x) — (7q,5(x— x") is the charge density. In this way, the partition function 
becomes 



Z 



^ \{nww n 



D0exp 



1 



(3.15) 



In a mean-field treatment, the functional integral over the auxiliary field introduced in 
(3.14) is approximated by the saddle point determined by the field equation 



(3.16) 



When we introduce the scalar variable <i> := iTcf), this equation becomes formally 
Gauss' law, with <I> the electrostatic scalar potential. The auxiliary field introduces in 
(3.14) may therefore be thought of as representing the scalar potential of the equivalent 
Coulomb gas 1 17|. 

On account of charge neutrality, we have the identity 



N 



[{N/2) 



112 



N 

n 



(3.17) 



where we recall that N is an even number. The factor A^!/[(A^/2)!]^ is the number of 
charge-neutral ter ms co ntained in the binomial expansion of the left-hand side. The 
partition function ( 3.15 ) may thus be written as [n7 1 



E 



(2z)_ 
iV! 



N 



N=0 

D(/)exp 



D0exp 
1 



{V(t)f - 2zcos(#) 



N 



(3.18) 
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where in the final form we recognize the sine-Gordon model. This is the dual theory we 



were seeking. Contrary to the original formulation (3.6), which contains the vortices 
as singular objects, the dual formulation has no singularities. To see how the vortices 
and the Kosterlitz-Thouless phase transition are represented in the dual theory we note 
that the field equation of the auxiliary field now reads 

iTV^fj) = 2Trzq (e*?"^ - e-*^'^) . (3.19) 



On comparison with the previous field equation ( 3.16[ ), it follows that the right-hand 



side represents the charge density of the Coulomb gas. In terms of the scalar potential 
Eq. ( 3.19| ) becomes the Poisson-Boltzmann equation 



= -2TTq {z e-'^"* - z e'^^*) , (3.20) 

describing, at least for temperatures above the Kosterlitz-Thouless temperature, a plas- 
ma of positive and negative charges with density n±, 

ze^'^"'^, (3.21) 

respectively. The fugacity z is the density at zero scalar potential. (It is to recalled that 
we suppress factors of a denoting the diameter of the vortex cores.) Equation ( 3.20[ ) is 



a self-consistent equation for the scalar potential $ giving the spatial distribution of the 



charges via (3.21 ). It follows from this argument that the interaction term 2z cos{q(j)) 
of the sine-Gordon model represents a plasma of vortices. 

The renormalization group applied to the sine-Gordon model reveals that at the 



temperature phase of tightly bound neutral pairs and a high-temperature plasma phase 



Kosterlitz-Thouless temperature Tkt — there is a phase transition between a low- 



of unbound vortices [133|. In the low-temperature phase, the (renormalized) fugac- 
ity scales to zero in the large-scale limit so that the interaction term, representing the 
plasma of unbound vortices, is suppressed. The long-distance behavior of the low- 
temperature phase is therefore well described by the free theory (V(/>)^/47r/3, repre- 
senting the gapless Kosterlitz-Thouless mode. This is the superfluid state. The expec- 
tation value of a single vortex vanishes because in this gapless state its energy diverges 
in the infrared. 

An important characteristic of a charged plasma is that it has no gapless excita- 
tions, the photon being transmuted into a massive plasmon. To see this we assume that 
g<I> << T, so that sinli(/3(7<i>) « /3(7$. In this approximation, the Poisson-Boltzmann 



equation (3.20) can be linearized to give 

(V^ - m^)$ = 0, ml= Airfizq^. (3.22) 

This shows us that, in contradistinction to the low-temperature phase, in the high- 
temperature phase, the scalar potential describes a massive mode — the plasmon. In 
other words, the Kosterlitz-Thouless mode acquires an energy gap mu- Since it pro- 
vides the high-temperature phase with an infrared cutoff, isolated vortices have a finite 
energy now and accordingly a finite probability to be created. This Debeye mechanism 
of mass generation for the photon should be distinguished from the Higgs mechanism 
which operates in superconductors (see below) and also generates a photon mass. 
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Another property of a charged plasma is that it screens charges. This so-called 
Debeye screening may be illustrated by adding an external charge to the system. The 
linearized Poisson-Boltzmann equation (3.22) then becomes 

(V^ - m2,)$(x) = -2^qo'5(x), (3.23) 

with qo the external charge which we have placed at the origin. The solution of this 
equation is giv en by $(x) = qoKo{mT:,\x.\) with Kq a modified Bessel function. The 
mass term in ( 3.23 ) is (2tt times) the charge density induced by the external charge, 
i.e.. 



And(x) = --!-go™D^o(™D|x|). 
ZTT 



(3.24) 



By integrating this density over the entire system, we see that the total induced charge 
/ Pind = —90 completely screens the external charge — at least in the linear approxi- 
mation we are using here. The inverse of the plasmon mass is the screening length — the 
so-called Debeye screening length. 

To see that the sine-Gordon model gives a dual description of a *He film we cast 
the field equation (3.16) in the form 

72 , 



-mgV X Vg 



(3.25) 



where we employed Eq. (3.5). On integrating this equation, we obtain up to an irrele- 
vant integration constant 



iTdi(t> = -qeijidjip - (p^). 



(3.26) 



This relation, involving the antisymmetric Levi-Civita symbol, is a typical one between 
dual variables. It also nicely illustrates that although the dual variable is a regular 
field, it nevertheless contains the information about the vortices which in the original 
formulation are described via the singular vortex gauge field <p^. 

Given this observation it is straightforward to calculate the current-current correla- 
tion function {gi{'k)gj{—'k)), with 



the mass current. We find 

(g,(k)g,(-k)) 



1 



where the average is to be taken with respect to the partition function 

1 



D0exp 



47r/3 



(V0)^ 



(3.27) 



(3.28) 



(3.29) 



which is obtained from ( 3.18[ ) by setting the interaction term to zero. We obtain in this 
way the standard expression for a superfluid 

1 



(3.(k)g,(-k)> 



f3 k2 



(^■ijk kikj 



(3.30) 



The 1 /k^ reflects the gaplessness of the (/)-field in the low-temperature phase, while the 
combination J^jk^ — kikj arises because the current is divergent free, V • g(x) = 0, or 
k • g(k) = 0. 
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3.2 Two-dimensional Superconductor 



We now turn to the dual description of a superconducting film. We thereto minimally 
couple the model of the preceding section to a magnetic field described by the magnetic 
vector potential A. For the time being we ignore vortices by setting the vortex gauge 
field if^ to zero. The partition function of the system then reads 

Z = j^'P j DAS(A) exp (^^(3 J , (3.31) 

where S(A) is a gauge-fixing factor for the gauge field A, and H is the Hamiltonian 



n 



(3.32) 



with 



Vs = — (V(y5 - 2eA). 



(3.33) 



The double charge 2e stands for the charge of the Cooper pairs which are formed at the 



bulk transition temperature. The functional integral over (f in (3.31 ) is easily carried 
out with the result 



Z = 



J 



DAS(A)exp<^-^ 



(V X A)-^ + m-^Ai Si-i 



(3.34) 



A, 



where the last term, with m\ — Ae^ps/m^, is a gauge-invariant, albeit nonlocal mass 
term for the gauge field generated by the Higgs mechanism. The number of degrees of 
freedom does not change in the process. This can be seen by noting that a gapless gauge 
field in two dimensions represents no physical degrees of freedom. (In Minkowski 
spacetime, this is easily understood by recognizing that in 1 + 1 dimensions there is 
no transverse direction.) Before the Higgs mechanism took place, the system therefore 
contains only a single physical degree of freedom described by if. This equals the 



number of degrees of freedom contained in (3.34) 



We next introduce an auxiliary field h to linearize the first term in (3.34), 



exp 



(V X A)^ 



— 1 D/iexp 



1 



h{V X A) 



(3.35) 



and integrate out the gauge-field fluctuations [with a gauge-fixing term (1/2q;)(V • 
A)^]. The result is a manifestly gauge-invariant expression for the partition function in 
terms of a massive scalar field h, representing the single degree of freedom contained 
in the theory: 



Z= I D/iexp^-^ 



(3.36) 



To understand the physical significance of this field, we note from (3.35) that it satisfies 
the field equation 

h = i(3V X A. (3.37) 
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That is, the fluctuating field h represents the local magnetic induction, which is a scalar 
in two space dimensions. Equation ( 3.36| ) shows that the magnetic field has a finite 
penetration depth Al = In contrast to the original description where the func- 

tional integral runs over the gauge potential, the integration variable in (3.36) is the 
physical field. 

We next include vortices. The penetration depth Al provides the system with an 
infrared cutoff so that a single magnetic vortex in the charged theory has a finite energy. 
Vortices can therefore be thermally activated. This is different from the superfluid 
phase of the neutral model, where the absence of an infrared cutoff permits only tightly 
bound vortex-antivortex pairs to exist. We expect, accordingly, the superconducting 
phase to describe a plasma of vortices, each carrying one magnetic flux quantum ±7r/e. 
The partition function now reads 



Z = 



V — 



Hj ^ J^fJ DAS(A) exp (^-13 J 7?j (3.38) 



where z is the fugacity, i.e., the Boltzmann factor associated with the vortex core en- 
ergy. The velocity appearing in the Hamiltonian (3.32) now includes the vortex gauge 
field 



1 



(V(p-2eA- (p^). 



(3.39) 



The vortex gauge field (p^ can be shifted from the first to the second term in the 
Hamiltonian (3.32) by applying the transformation A — > A — (p^ /2e. This results in 
flie shift 

VxA^VxA-B^, (3.40) 



with the plastic field 



5" = -*o5I^"a'5(x-x") 



(3.41) 



representing the magnetic flux density. Here, $o — vr/e is the elementary flux quan- 
tum. Repeating the steps of the previous paragraph we now obtain instead of (3.36) 



Z 



E 

N±=0 

xexp 



nTinJ. 



n 



Dh 



1 



(3.42) 



where h represents the physical local magnetic induction h 

h = if3{V X A- B^) ^ i(3h. 



(3.43) 



The field equation for h obtained from ( |3.42[ ) yields for the magnetic induction: 

~V'^h + m\h = m\B^ , (3.44) 
which is the familiar equation in the presence of magnetic vortices. 
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The last term in (3.42) shows that the charge g with which a magnetic vortex cou- 



ples to the fluctuating /i-field is the product of an elementary flux quantum (contained 
in the definition of B^) and the inverse penetration depth niA — 1/Al, 



(3.45) 



For small fugacities the summation indices and iV_ can be restricted to the values 



0, 1 and we arrive at the partition function of the massive sine-Gordon model [134| 



Dh exp 



1 

2P 



1 



2zcos ( ^oh 



(3.46) 



This is the dual formulation of a two-dimensional superconductor The magnetic vor- 
tices of unit winding number Wa = ±1 turned the otherwise free theory (3.36) into an 
interacting one. 



The final form (3.46) demonstrates the rationales for going over to a dual theory. 
First, it is a formulation directly in terms of a physical field representing the local 
magnetic induction. There is no redundancy in this description and therefore no gauge 
invariance. Second, the magnetic vortices are accounted for in a nonsingular fashion. 
This is different from the original formulation of the two-dimensional superconductor 
where the local magnetic induction is the curl of an unphysical gauge potential A, and 
where the magnetic vortices appear as singular objects. 

Up to this point we have discussed a genuine two-dimensional superconductor As 
a model to describe superconducting films this is, however, not adequate. The reason 
is that the magnetic interaction between the vortices takes place mostly not through 
the film but through free space surrounding the film where the photon is gapless. This 
situation is markedly different from a superfluid film. The interaction between the 
vortices there is mediated by the Kosterlitz-Thouless mode which is confined to the 
film. A genuine two-dimensional theory therefore gives a satisfactory description of a 
superfluid film. 

To account for the fact that the magnetic induction is not confined to the film and 
can roam in outer space, the field equation (3.44) is modified in the following way 
pl|l36|] 



(3.47) 



Here, 1/A^ = dm\, with d denoting the thickness of the superconducting film, is 
an inverse length scale, denotes the coordinates in the plane, h the component of 
the induction field perpendicular to the film, and ddix^) is a smeared delta function of 
thickness d along the xs-axis 



= for 1 > d/2 
+ for 1 < d/2 



(3.48) 



The reason for including the smeared delta function at the right-hand side of (3.47) is 



that the vortices are confined to the film. The delta function in the second-term at the 
left-hand side is included because this term is generated by screening currents which 
are also confined to the film. 
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To be definite, we consider a single magnetic vortex located at the origin. The 



induction field found from (3.47) reads 



with Jq the 0th Bessel function of the first kind. At small distances from the vortex 
core (X±q >> 1) 

fe(x±,0)- . (3.50) 
47rA± x_L 



while far away (X±q << 1) 



h{^^,0)^^^. (3.51) 



This last equation shows that the field does not exponentially decay as would be the 
case in a genuine two-dimensional system. The reason for the long range is that most 
of the magnetic interaction takes place in free space outside the film where the photon 
is gapless. If, as is often the case, the length Aj^ = l/dm\ is much larger than the 
sample size, it can be effectively set to infinity. In this limit, the effect of the mag- 



netic interaction, as can be seen from (3.50), diminishes and the vortices behave as 



in a superfluid film. One therefore expects a superconducting film to also undergo a 
Kosterlitz-Thouless transition at some temperature Tkt characterized by an unbinding 
of vortex-antivortex pairs. The first experiment to study this possibility was carried out 



in Ref. 1137|. Because the transition temperature Tkt is well below the bulk temper- 
ature Tc where the Cooper pairs form, the energy gap of the fermions remains finite 
at the critical point |]138|]. This prediction has been corroborated by experiments per- 



formed by Hebard and Palaanen on superconducting films [139]. For temperatures 
Tkt < T < Tc, there is a plasma of magnetic vortices which disorder the supercon- 
ducting state. At Tkt vortices and antivortices bind into pairs and algebraic long-range 
order sets in. 



3.3 Bosonization 



The two dual models ( 3.36 ) and ( 3.46 ) we encountered in the preceding section de- 
scribing a superconducting film without and with vortices included are reminiscent 
of the bosonized massl ess and massive Schwinger model, respectively. The massless 
Schwinger model [116] describes gapless fermions interacting with an electromagnetic 
field in one space and one time dimension. It is defined by the Lagrangian 



lp2 

4 /J,l^ 7 



(3.52) 



where ^= A^7^ with 7^ the Dirac matrices (2.68), ^ is a two-component Grassmann 
field describing the fermions, and ^ = Since the theory is bilinear in tp, the 

functional integral over the Grassmann fields in the partition function 



D'^/JD^P / DA^E{A^) exp 



C 



(3.53) 
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with 5(A^) a gauge-fixing factor, can be easily carried out. It yields a functional 
determinant 



DipDip exp 



Det - e4) 



(3.54) 



which can be evaluated in closed form using, for example, the derivative expansion 



[|140p. One finds 

Det (/^ — e^) — exp 



52 



The partition function (3.53 ) can now be written in a form. 



(3.55) 



(3.56) 



X exp • 



'-Ff. 



-m\A 



92 



A. 



which is the Minkowski analog of the partition function (3.34) describing a supercon 



ducting film without vortices. Equation (3.56) shows us that the photon has acquired a 



mass niA, given by m\ = e^/tt, where it should be borne in mind that in one spatial 
dimension the electric charge e has mass dimension one. This mass generation due to 
gapless fermions is called the Schwinger mechanism. 

We have seen that the partition function ( 3.34 ) could be equivalently represented 
by the massive scalar field theory ( 3.36 ). T he sa me holds for the Schwinger model. 
Rather tha n line arizing the Maxwell term in ( 3.56 ), which would be the analog of what 
we did in ( 3.35 ), we represent the term obtained in ( 3.55 ) as a functional integral over 
aBose field [|141|: 



DijjD'ip exp 

00 exp ■{ i 



1 



(3.57) 



We recognize here two of the bosonization rules ( 2.78 ) we employed before, viz.: 

—;=e^„du<p, (3.58) 



where — ijj'^f_iip is the electromagnetic current. The integral over the gauge field can 
now be easily carried out, say in the Lorentz gauge 9^ = 0, to yield the bosonized 
form of the massless Schwinger model 



Z = / D0cxp 



(3.59) 
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This shows that the model is equivalent to a massive scalar theory. It is the Minkowski 
analog of the dual theory (3.36) of a superconducting film without taking into account 
vortices. 

To understand the physical origin of the scalar field and why it is massive, it 
should be noted that the massless Schwinger model has besides the local U(l) gauge 
symmetry 

(3.60) 



also a global U( 1 ) chiral symmetry. Under chiral transformations 

where A is the transformation parameter and 75 denotes the matrix 75 
corresponding Noether current is the axial current 



(3.61) 
17071. The 



j'jl = #7m75V'- (3-62) 

At the quantum level, the chiral U(l) symmetry is spontaneously broken by a fermion 
condensate [142| 



Zn 



with 7 ; 



0.577216 Euler's constant. By virtue of the identity 
the axial and electromagnetic current are related via 



With the bosonization rule (3.58) we then obtain the correspondence 



(3.63) 



(3.64) 



(3.65) 



(3.66) 



identifying as the Goldstone field of the broken U(l) chiral symmetry. This can, 
however, not be the end of the story [143|. Goldstone modes are by definition gapless, 
but Eq. (3.59) shows that </) has acquired a mass, meaning that the chiral symmetry 
is no longer an exact symmetry of the bosonized theory. The axial current, which is 
classically conserved, has an anomaly at the quantum level as can be seen by invoking 
the field equation obtained from ( 3.59| ), 



(3.67) 



The destruction of the chiral symmetry at the quantum level can also be seen directly in 
the fermionic formulation, where it can be attributed to the effect of one-loop graphs. 
The anomaly (3.67) follows there from integrating out the fermionic degrees of free- 
dom. Using, for example, the derivative expansion [140|, one finds 



niA 



da A,. 



(3.68) 
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On comparing the two expressions for the anomaly, we see that the Bose field i 
sents the dual field strength F := e^iydf^Ai,, 



1 

ruA 



-F. 



I repre- 



(3.69) 



This is the Minkowski analog of the relation ( 3.37| ) we found in the context of a super- 
conducting film. 

We next consider possible vortex contributions. Let us go back to the functional 
determinant appearing in (3.55) and evaluate it in the presence of a vortex. To this end 
we go over to Euclidean spacetime where a vortex becomes an instanton. There exists 



a powerful index theorem [144| relating the zero eigenvalues of the massless Dirac 
operator — e^to the so-called Pontryagin index 



(3.70) 



An 



which is the winding number of the vortex. The theorem states that the index of the 
Dirac operator, defined by the number n+ of zero-modes with positive chirality minus 
the number n_ of zero-modes with negative chirality, is given by the Pontryagin index. 



(3.71) 



In the presence of vortices Q ^ Q, so that there are always zero-modes. Since the 
determinant of an operator is the product of its eigenvalues, the determinant appearing 
in (3.55) is zero when vortices are included. Hence, in the massless Schwinger model, 
vortices do not contribute to the partition function. The suppression of instanton con- 
tributions by gapless fermions is a common feature [{77|]. 

The situation is different in the massive model, defined by the Lagrangian 



1 p2 

4 fJ.!^' 



(3.72) 



First of all, the mass term —mipip in the Lagrangian explicitly breaks the chiral symme- 
try. So, even if there were no anomaly, we expect a massive bosonized theory. Second, 
because the massive Dirac operator contains no zero-modes, the index theorem is of 



no relevance here and vortices should be included in (3.59). The bosonized partition 



function then becomes, cf. (3.42) 



N±=0 



^N++N- 

n+inJ. 



n 



D^ 



(3.73) 



2 ™2.2]^^/^P^ 



where z = exp(iS'c), with Sc the vortex core action. The plastic field F^ describes the 
vortices, 

F^ ^ -'S^o'^WcSix - x"), (3.74) 
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where the sum is over all vortex locations and Wa is the winding number of the ath 
vortex. The elementary flux quantum in the Schwinger model is $0 = 27r/e, which 
is twice the value in a superconductor because the electric charge of a Cooper pair 
is twice that of an electron. We shall consider only vortices of unit winding number 
Wa = ±1. For z small we can restrict ourselves to configurations having zero or one 



vortex. The partition function then becomes the Minkowski analog of (3.46) 



Z = y"D0exp (^z^ji [(a^0)2-m^02] _2^^os(g0)|^ ^ (3.75) 



with g — $o"iyi as in (3.45). From this we see that the inclusion of the mass term in 
the Schwinger model resulted in the cosine term in the bosonized theory. That is, we 
have the correspondence 

m'ip'tp 22: cos (^ViTTcjyj , (3.76) 

where we used that 

g = ^orriA = \/47r (3.77) 



in the Schwinger model. Surprisingly, (3.76) is the remaining bosonization rule of 



(2.78). Instanton contributions, which were suppressed in the massless model, generate 



the interaction term in the bosonized form of the massive Schwinger model. We note 



that even in the absence of the anomaly, implying that the mass term in (3.75 ) would 
be absent, the bosonized theory would still be massive. This follows from noting that 
the expansion of the cosine term contains a term quadratic in cj). 



3.4 Ginzburg-Landau Model 

In this section, we give arguments as to why the Ginzburg-Landau model of a three- 
dimensional superconductor is not well suited to estabhsh the critical properties of a 
type-II superconductor at zero external field. 
The model is defined by the Hamiltonian 

n=\{W- 2zeA)(/.|' + + X\cj>\^ + i(v x A)^ + ^(V • A)^, (3.78) 

where the coefficients 2e and are the electric charge and mass of the complex 
0-field, while A is a coupling constant. The mass term changes sign at the critical 
temperature Tc. We have added a gauge-fixing term with parameter a. To acquire 



a physical understanding of what the Hamiltonian (3.78) describes, we bring to the 
attention the well-know fact jl7| ] that a |</)|'' -theory gives a field theoretic description of 
a gas of closed loops with contact repulsion. This equivalence is based on Feynman's 



observation [145| that the Green function 

r pikx 

of the free theory with positive mass term, 

7^0 = |V0|2 + m2|0|2, (3.80) 
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can be expressed as a path integral. (In contrast to the previous two chapters, we in this 
and the following chapter denote the free Green function by G without an index 0.) To 
this end, the Schwinger proper-time method is invoked to write the right-hand side of 
(3.79) as an integral over the proper time r 1 146|: 

G(x) = / dre' 



g»k.Xg-rk 



dr e 



(-) 

\4ttt J 



3/2 



where we used the identity 



1 
a 



dr e 



(3.81) 



(3.82) 



Accordi ng to the path-integral formulation of quantum mechanics [ 147|, the right-hand 
side of ( 3.81 ) can be represented as a sum over all paths of a relativistic particle with 
mass running from at imaginary time to x at time r: 



G(x) = / dr 



x(r)— X 



(o)=o 



X>x(r')e" 



-So 



The (Euclidean) action 



Jo 



(3.83) 



(3.84) 



with x(r) — dx(r)/dr, is an elusive representation of (m^ times) the length L of the 
path 

m^L^m^j drVx2(r'), (3.85) 

with now being interpreted as the line tensio n. T o establish this connection we 
consider the canonical momentum obtained from (3.85) 



(3.86) 



which is seen to satisfy the constraint 

p2 = ml. (3.87) 

To incorporate this constraint we go over to the Hamilton formalism and write instead 
of (pSl) 



m^L — > / dr' [p • x — a{p 
Jo 



' ml] 



(3.88) 



where the constraint is implemented with the help of the Lagrange multiplier a. The 
path integral is now over phase space. Because the action is invariant under reparame- 
terization, i.e., under the transformation r t{t), this multiplier can be given a fixed 
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value; a = 1, for example. The integral over p is Gaussian and immediately yields the 



action (3.84) 



The partition function (with fields and coupling constants rescaled such that no 
explicit temperature dependence appears in the Boltzmann factor) 



= / D0*D<^ exp[- I Ho], (3.89) 

involves only closed paths: 

,rT^„.^„2 , ^^^^ 



In(Zo) = - ln[Det(p2 + ml)] = -Tr ln(p2 + ml) (3.90) 



k 



with the action ( [3.84[ ). 

The 101^-interaction in the Ginzburg-Landau model results in an additional term 

Sx = -X f dr'dr" S [x(t') - x(r")] (3.91) 

"'0 

in the action, which gives an extra weight each time two loops — one parameterized by 
r' and one by t" — intersect. Physically, it represents a repulsive contact interaction 
between loops. Finally, the coupling of the field (j> to the magnetic vector potential A 
via the electric current 

je = -2ei0* V 0- 2(2e)2A|(/i|2 (3.92) 
with a charge 2e, results in the additional term 

Se = 2ie [ dr' X • A[x(t')], (3.93) 
Jo 

showing that the loops described by the Ginzburg-Landau model are electric current 
loops. 

On entering the superconducting phase, characterized by a negative mass term 
(m| < 0), the electric current loops proliferate and \(p\ develops a vacuum expecta- 
tion value. In the London limit, the fluctuations in the modulus of cf) are neglected and 
the field is written as 

0(x) = 4=6'"^''^ (3-94) 
v2 

with (j) a constant. Formally, this corresponds to taking the limit —oo and 

A ^ oo in the Hamiltonian such that 



101' = -^ (3.95) 
is finite. Physically, denotes the mass density ps of the superconducting electrons. 
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Let us first define an order parameter describing the superconducting state. The 
field (f> as it stands can not play this role because it is not gauge invariant. But it can be 
used to construct a gauge-invariant Mandelstam-like operator in the following way 

W^(L,)=exp|i (^(z) + J A(x)-J(x) |, (3.96) 

where tp is the phase of cj) and J describes an external current line originating in z, 

V- J(x) = 2e(5(x-z), (3.97) 

or 



Ji{x) = 2e y" dy,5{x ~ y). 



(3.98) 



Here, denotes a path running from z to infinity. The second term in ( |3.96| ) is incor- 
porated to render the operator gauge invariant. Indeed, under a gauge transformation 



A(x) -> A(x) + Vq!(x), ip{z) (p{z) + 2ea(z), 
the operator W{L^) is invariant 



W[L^) W{L^) exp 



ia{-i) + i I Va(x) • J(x) 

X 



(3.99) 



(3.100) 



where in the last step we performed an integration by parts. The gauge invariance of 
H^(Lz) can be made more explicit by writing it in the equivalent form 



exp 



2^ 



{Vip - 2eA) • J 



(3.101) 



To show that this operator indeed is the order parameter of the superconducting 
state, we calculate the correlation function 

{W{L^)W*{L-,)) = J DADp W{L,)W*{L^) exp (-/^) > (3-102) 

obtained by introducing besides an external current line originating in z also one ter- 
minating at z. The Hamiltonian H is the Ginzburg-Landau Hamiltonian (3.78) in the 
London limit. Since both integrations are Gaussian, they are easily carried out with the 
result [|148[ 



{WiL^)W*{L^)) 
1 



(3.103) 



^p(x) G(x - y) p{y) + J,(x) G(x - y) J,(y) 



where p(x) = 2e[(5(x — z) — (5(x — z)] denotes the electric current source and sink, 
and 



G(x) 



k2 



47r Ixl 



(3.104) 
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is the scalar Green function, with 



rriA = 2e|(^| 



(3.105) 



the photon mass. 

In the superconducting phase, the electric current lines proliferate and carry no en- 
ergy; they are of no physical significance, only the endpoints are. The external current 
J can then be written as a gradient of a potential U, 



J = -WU, 

with V^?7(x) = —p{x). In this way, the correlation function becomes 

{W{z)W*{i)) = exp 



TT-T / PWG'o(x-y)p(y) 



(3.106) 



(3.107) 



where Go (x — y) is the gapless scalar Green function obtained from the massive one 
(3.104) by taking the limit niA — > 0. For x = y we have a diverging self-interaction 
which is irrelevant and can be eliminated by defining a renormalized operator Wr via 



W,{z) = H/(z)exp 



1 



rG(0) 



We then find for this operator 



(W^r(z)W;(z)) =exp 



(3.108) 



(3.109) 



where L^z. is a path connecting the source at z with the sink at z, and {L^z \ is the length 
of the path. For large separation this correlation function tends to one, implying that 
the operator Wr{z) develops an expectation value in the superconducting phase. 

It is important to realize that this operator is gauge invariant. Frequently, the super- 
conducting phase is referred to as a phase of spontaneously broken gauge symmetry 
[149]. This should, however, not be taken too literally. A celebrated theorem due to 
Elitzur 1 150 1 states that a local gauge symmetry can never be spontaneously broken. 
An operator that by developing a nonzero vacuum expectation value would sponta- 
neously break the gauge symmetry does not exist. To break the symmetry, the operator 
must transform under the gauge group, but then it cannot develop an expectation value 
because only gauge-invariant objects can. 

In the normal phase, electric current lines carry energy and therefore become phys- 
ical. The form (3.106) is then no longer applicable and we have to use Eq. (3.98) in- 
stead. We shall assume that the current lines have a certain width of the order l/|m^|. 



where the mass will be identified later on in this chapter (see page 106). To pre- 
vent infrared divergences, we give the magnetic vector potential A a small mass fi. The 
correlation function ( 3. 103) then becomes after renormalization 



{Wr{Lz)W;{Lz)) = cxp(-Mw|izz|)exp 



1 /2e\^ e-^l^-^l 



47r \ /i 



(3.110) 
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Figure 3.2: One-loop mass correction. The straight and wiggly lines represent the </>- 
and A-field Green functions, respectively. 



where now L^^ denotes the shortest path connecting the source at z with the sink at z, 
while |Lzz| denotes the length and Mw the line tension of the current line. 



An 



(3.111) 



We see that in the limit |Lzz| oo, the expectation value vanishes. Although the 
operator W distinguishes the superconducting phase from the normal phase, having a 
zero expectation value in the normal phase and a nonzero one in the superconducting 
phase, it is not an order parameter in the strict sense of Landau in that it makes no 
symmetry statement. 

A few remarks are in order. First, the line tension AIw diverges when the arti- 
ficially introduced photon mass fi is taken to zero. This means that the current line 
connecting the source and sink becomes infinitely heavy. Isolated sources and sinks 
can therefore not exist in the normal phase. They are confined in neutral configurations 
of tightly bound pairs. In Minkowski spacetime, where an electric current line becomes 



the worldline of a charged particle, the energy (3.111) denotes the infrared-diverging 



selfenergy [ 15 1 1 of that particle obtained from evaluating the Feynman graph depicted 



in Fig. 3.2. Second, the alert reader mi ght wo nder how the result (3.109) obtained in 
the normal phase is connected with Eq. ( 3.110 ) obtained in the superconducting phase. 
In particular, the line tension AIw of the current line seems n ot to v anish at the critical 
temperature, implying that the first exponential function in ( 3.110 ) would survive the 
transition to the superconducting phase, and that the operator W would vanish there. 
This is in apparent contradiction with previous statements that in the superconducting 
phase, current lines carry no energy and proliferate. This issue will be addressed be- 
low (see page |10^ ) after we applied a duality transformation to the Ginzburg-Landau 
model. 

We continue to investigate the phase transition by taking into account fluctuations. 
The one-loop effective potential is readily calculated in the London limit (3.94). In the 
gauge V • A = 0, one obtains 



Veff(<^) 



ln(A:'^ + m\) 



fc<A 



(3.112) 



We regularized the integral by introducing an ultraviolet cutoff A. An irrelevant ultra- 
violet divergence of the type A'^ log A is ignored, while a divergence of the type A|^p 
is absorbed into a renormalization of the mass parameter, as usual. The effective poten- 



tial Veff (0) has to be added to the tree potential Vq = 



It results in 



a 1 st-order phase transition taking place at a temperature Ti above the temperature 



where the mass term of the Ginzburg-Landau model changes sign (see Fig. 3.3). The 
(positive) value of m? at the 1 st-order transition point, obtained by equating the sum 
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Figure 3.3: Effective potential up to one-loop order. 



of the tree and the one-loop potential and also its first derivative with respect to to 



zero, is given by [152| 



(3.113) 



2 32 

At small values of the so-called Ginzburg-Landau parameter kql^ which is defined as 
the ratio of the two mass scales in the theory. 



iriA 



Al 



(3.114) 



the reliability of this pertu rb ative result has been verified numerically [ 153 1 . It has been 
argued by Kleinert 1 154 1 that for larger values this result is, however, not meaningful 
anymore. At some critical value of the Ginzburg-Landau parameter kql, the Ist-order 
transition goes over into a 2nd-order transition. 

Let us investigate the superconducting-to-normal transition further by applying 
renormalization-group analysis to the Ginzburg-Landau model. This was first done 
by Halperin, Lubensky, and Ma [155| in an e(= 4 — (i)-expansion around the upper 
critical dimension (1 = 4. We shall instead use the fixed-dimension approach of Parisi 
[156, 157]. The model (|3.78[) involves only a single complex field. However, for rea- 



sons that will become clear when we proceed, it is often generalized to contain 2n real 
components: 

/ 01 + 402 

^ ' (3.115) 







V2 



V 



The quantities appearing in ( 3.78 ) should be given an index to indicate that they are 
bare quantities which will change in the process of renormalization. The renormahzed 
fields and parameters are related to the bare ones via 



A = Z 



-1/2 



AO, 



e = Z / eo, 



, ^-1/2 ,0 



(3.116) 



We shall work at criticality and to that end set the mass to zero. To avoid infrared 
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divergences, diagrams are evaluated at finite external momentum k. The renormaliza- 
tion factors can be used to define the following renormalization-group functions 



= i^S- ^Za] 

OK 



70 = 1^-^ ln(^0) 



(3.117) 



where the subscript is to indicate that the bare coupling constants eo, Aq are kept fixed. 
The dimensional parameter k plays the role of renormalization group scale parameter. 
The critical points and A* are determined by the zeros of the renormalization-group 
beta functions. 



OK 



/3a 



K — \ 

OK 



(3.118) 



where = e^K'' 



* and A ~ Xk'^~'^ are the dimensionless coupling constants. When 
evaluated at the critical point, the renormalization-group functions 7^ and 7^ yield the 
anomalous dimension rjA and 77^ of the A- and i/i-field, respectively. 

The various renormalization factors are readily evaluated in ar bitrary dimension 
2 < c? < 4, using the fixed-dimension approach of Parisi [156, 157]. To the one-loop 
order we find 



Za 



Zd, 



c{d) 



Zx = 



d-1 

l + c(d)[(d-l)-a(d-3)]e2 

1 + c{d) \{2n + 8)A - 2a{d - 3)6^ + ^d{d - l)e^/X 



where c(d) stands for the one-loop integral 



kk2(k + q)2 



q2 = l 



T{2~d/2)T^{d/2-l) 
(47r)^/2r(d - 2) ' 



(3.119) 
(3.120) 
(3.121) 



(3.122) 



In deriving this we used analytic regularization to handle the ultraviolet divergences. 
With the one-loop expressions for the renormalization factors, we obtain for the beta 
functions 



= (d-4) 



— n 



c{d) 
d^l' 



(3.123) 



Px = {d- 4) {a - c{d) (271 + 8)A^ - 2{d - l)Xe^ + \d{d - l)e^ } . 

(3.124) 

Note that while Zx and depend on the gauge-fixing parameter a, the function /3a 
extracted from them is independent thereof. We expect this cancelation of gauge- 
dependence in Px to persist in all orders of perturbation theory. Although the beta 
functions are gauge independent, the locations of the fixed points are not universal. 
They depend on the regularization chosen as well as on the renor malization prescrip- 
tion. In the Hmit d ^ 4 — e our formulas reduce to the know form [155, 158[. 
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Figure 3.4: Schematic representation of the /3e2 -function. 



Remarkably, [3^2 is independent of the self-coupHng A. By itself it gives rise to two 
fixed points (see Fig. 3.4): 





{d-l)/nc{d). 



(3.125) 



The neutral fixed point is ultraviolet stable, while the c harged fixed point is infrared 
stable in the e^-direction. Since also Za as given in Eq. ( 3.119 ) and thus are inde- 
pendent of A, the anomalous dimension rjA of the gauge field can be calculated with- 
out any knowledge of the location of the fixed points along t he A- direction. For the 
infrared-stable (IR) fixed point one finds at the one-loop order [ 155 1 



r]A=lA\ 



= A-d. 



(3.126) 



It was pointed out by Herbut and Tesanovic [159| that this is an exact result which 
follows directly from the definition of the /3e2 -function in terms of the renormalization 
factor Za 



Pe2 = 

OK 



^[(d-4)+7A] 



(3.127) 



That is to say, there are no higher-order corrections to the one-loop result ( 3.126 ). We 
have checked this explicitly to the two-loop order in the €{— 4 — (i)-expansion. The 
anomalous dimension is such that the dimension c?^ of the gauge field, which by naive 
power counting equals \{d— 2), is exactly unity. 



c?A = 1 for 2 < d < 4. 



(3.128) 



It is as if we were in the upper critical dimension d = 4. 

Since the result is independent of the number 2n of field components contained in 
the theory, this conclusion can be confirmed in the 1 / n-expansion. To leading order in 
1/n, the inverse gauge field propagator reads in the gauge V • A = and for 2 < < 4 



n„(q) 



c{d) 



(3.129) 



with Pij (q) the transverse projection operator 



(3.130) 
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The exponent d— 2 of the momentum in (3.129) equals d— 2d a, so that we again arrive 



at the conclusion that = 1 in 2 < d < 4. 

It implies that near the charged fixed point, the electric charge scales with the cor- 
relation length ^ in a way one expects from naive power counting 

= ZacI ^ ^''-^ (3.131) 

since Za scales as 

The exact result d^ = 1 has far reaching consequences as it implies that the term 
(V X A)^ in the Hamiltonian, which by naive power counting is marginal, has dimen- 
sion 4. It is therefore an irrelevant operator for all d below the upper critical dimension 
and can be omitted. Formally, this corresponds to taking the limit e — > cx) in the 
Ginzburg-Landau model as can be seen after a rescaling of the gauge field A ^ A/e. 
Without the term (V x A)^, the Ginzburg-Landau model resembles the CP"^^ model, 

/:cp = l(V-^A)0|^ (3.132) 

which has the additional constraint 

\(t>\^^l. (3.133) 

As regards to critical exponents this difference is however irrelevant. Indeed, the well- 
studied 0(n) </>^-model and its nonlinear version, the 0(n) nonlinear sigma model 



where the fields satisfy a condition like (3.133), are generally accepted to be in the 



same universality class [160|. Whereas the linear model is usually investigated in an 
e-expansion around the upper critical dimension (d = 4), the nonlinear model is often 
treated in an expansion around the lower critical dimension (d = 2), so a direct com- 
parison is not possible. But the 1/n-expansion, being the same for both theories and 
applicable in arbitrary dimension 2 < d < 4, can be used to bridge these two regions. 
The results obtained in this way for d close to the upper and lower critical dimension 
are identical to those obtained in the respective e-expansions for large n. 

The 1/n-expansion has also been employed to argue that the Ginzburg-Landau 
model with 2rt-components and the CP"^^ model, too, are in the same universality 



class [161, 162]. Our conclusion that below the upper critical dimension, the term 
(V X A)^ in the Ginzburg-Landau model is irrelevant gives additional support to this 
conjecture. 

For a superconductor, having one complex field {n — 1), this implies that the 



relevant model is the CP° model. Although the e(= d — 2)-expansion [161] predicts a 



nontrivial 2nd-order transition, this model has no degrees of freedom. Physically, this 



means that the charged degrees of freedom are irrelevant at the phase transition [163|. 
In itself this is no reason to reject the Ginzburg-Landau model as a basis for studying the 
phase transition — the "unphysical" limit n ^ of the 0(n) model describes the long- 
distance behavior of a flexible polymer But what is more worrying is that vortices, 
which we saw to be of paramount importance in two dimensions, are not accounted 
for in the discussion so far In addition, there are technical problems with using the 
Ginzburg-Landau model to study the critical properties of the phase transition. 
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2nc(d) 




Figure 3.5: The critical number of components 2nc as a function of the dimensionality 
d below which the one-loop calculation predicts a Ist-order transition. 



Namely, the two beta functions yield no IR fixed point unless the number of field 
components is taken large enough. Below a critical number ndd), determined by the 
condition that 



Ad = 2y/n^ - 2{d - 2){d - l)2(d + l)n - 4(d - l)^{d + 1) (3.134) 
be real, the system has no IR fixed point and the one-loop calculation predicts a Ist- 



order transition. In Fig. 3.5 this number is plotted as function of the dimensionality d. 
We see that the number decreases when the dimensionality is reduced. More precisely, 

2nc(4) = 12(15 + 4Vl5) w 365.9 (3.135) 
2nc(3) = 16(2 + \/6) « 71.2 (3.136) 
2nc(2) = 4V3w6.9, (3.137) 



where the d = 4-result is due to Halperin, Lubensky, and Ma [155|. The figure also 
shows that for a fixed number of field components, 2nd-order behavior is favored when 
d ^ 2, while in the opposite limit, d ^ 4, Ist-order behavior is favored. However, 
even close to the lower critical dimension, where the CP"~^ model predicts a 2nd- 



order phase transition all the way down to n = 1 [161], the Ginzburg-Landau model 
still requires more than 2\/3 « 3.5 complex fields in order to have an IR fixed point. 
The status of the results obtained from the Ginzburg-Landau model in low-order per- 
turbation theory is therefore not clear. 

Given the two-dimensional results, showing the importance of vortices, it is natural 
to consider the dual theory of the three-dimensional Ginzburg-Landau model as an 
alternative to study the phase transition. As before, the dual formulation focuses on 
the vortices in the system, which in three dimensions are line defects. It will again 
turn out that the dual formulation is one directly in terms of physical variables and for 
that reason involves only a global, not a local, symmetry. We will see that the phase 
transition of a type-II superconductor at zero external field becomes more tractable in 
the dual formulation than in the original Ginzburg-Landau formulation. 
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3.5 Dual Ginzburg-Landau Model 



In this section, we apply a duality transformation to the Ginzburg-Landau model. The 
basic idea of this approach originates from three-dimensional lattice studies carried out 
two decades ago [164, 165, 166, 167|. These studies were instigated by the success of 
the Kosterlitz-Thouless theory of the phase transition in a superfluid film [13C, 131 1, 



which we discussed in Sec. 3.1. The three-dimensional lattice studies of the neutral 



[164] and charged a;?/-model 1 165} , 166] were aimed at obtaining a dual description 



in terms of vortices. Following a suggestion by Helfrich and Miiller [168], Dasgupta 



and Halperin [169] carried out a lattice simulation of the superconductor-to-normal 



phase transition in the dual formulation. Their study revealed that for small values of 
the electric charge e, which implies large values for the Ginzburg-Landau parameter 
( [3.114| ), the transition was 2nd-order with xy-exponents. A detailed account of these 
matters as well as an extensive list of references to the literature can be found in Ref. 

Another development underscoring the im port ance of the dual approach to the 
Ginzburg-Landau model was initiated in Ref. ]17C] . Th e basic observation was that 
since a local gauge symmetry can never be broken [150], a local gauge description of 
a phase transition is not feasible. It was argued that the three-dimensional Ginzburg- 
Landau theory contains, in addition to the local gauge symmetry, another global U(l) 
symmetry. When considered in 2 + 1 space-time dimensions with a Minkowski metric, 
this symmetry is generated by the magnetic flux operator. It was demonstrated that this 
symmetry is broken in the normal phase, while it is unbroken in the superconducting 
phase. A genuine order parameter in the sense of Landau was given, and it was shown 
that the massless photon of the normal phase is the Goldstone mode associated with 
the broken flux symmetry. 

To see in which regime magnetic vortices or Abrikosov flux tubes are important, 
we introduce a magnetic monopole at some point z inside the system. A monopole is 
a source of magnetic flux. Due to the Meissner effect, the flux lines emanating from 
the monopole are squeezed into a flux tube. In this way we have created a magnetic 
vortex at zero external field. Electrodynamics in the presence of a monopole was first 
described by Dirac [171] who argued that the combination 



V X A(x) 



BP(x) 



is the physical local magnetic induction h. The subtracted plastic field 



(3.138) 



(3.139) 



with $0 = Tr/e the magnetic flux quantum, removes the field of the so-called Dirac 
string running along some path L^, from the location z of the monopole to infinity. On 
account of Stokes' theorem, the plastic field satisfies the equation 

V • BP(x) = $o^(x- z). (3.140) 

In the presence of the monopole, the Ginzburg-Landau Hamiltonian becomes in the 
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London limit and after integrating out the phase field Lp 



= i(V X A-B 



Pn2 



A, 



1 

2^ 



(V-A)2, (3.141) 



where we gave Ti, the superscript P to indicate the presence of the monopole. Let us 
first verify Dirac's assertion that the combination (3.138) describes a point monopole 
with its Dirac string removed and consider the field equation for the gauge field. 



^,(x)= / G,,(x-y) [VxBP(y)], 



The gauge-field Green function Gy appearing here is 

6ij (^kikj^ /]^" 



«j(x) 



k2 



a 



(3.142) 



(3.143) 



The local magnetic induction corresponding to the classical solution given in (3.142) 
is 

V X A(x) - BP(x) = $oVG(x - z) - m\ f G(x - y)B^(y), (3.144) 

Jy 



where G is the Green function ( 3. 104 ). The first term at the right-hand side corresponds 
to a screened Coulomb force generated by the monopole. The last term, which is 
only present in the superconducting phase where the Meissner effect is operating and 
niA 0, describes the magnetic vortex. If we calculate from the right-hand side of 
(3.144) the magnetic flux through a plane perpendicular to the Dirac string, we find 
that precisely one flux quantum pierces the surface in the negative direction 



y d^x, <i>oa,G(x -z)-m\J G(x - y)Bf (y) 



(3.145) 



Here, d Xi denotes an element of the surface orthogonal to the Dirac string. Equation 



( 3.145 ) confirms Dirac's assertion that the magnetic flux emanating from a monopole 
must be supplied by an infinitesimally thin string of magnetic dipoles and that in order 
to obtain the true local magnetic induction of a genuine point monopole, this string 
has to be subtracted. Whence, the B^-term at the left-hand side of (3.144). While the 



Dirac string is immaterial in the normal phase, the last term at the right-hand side of 
(3.144) shows that it acquires physical relevance in the superconducting phase where 
it serves as the core of the Abrikosov flux tube [172|. 



The energy Ey of this configuration is obtained by substituting the solution (3. 142) 
back into the Hamiltonian ( |3.141 ). It is divergent in the ultraviolet because in the 
London limit, where the mass \m^\ of the ^-field is taken to be infinite, the vortices 
are considered to be ideal lines. For a finite mass, a vortex core has a typical width of 
the order of the coherence length <^ — l/jm^l. This mass therefore provides a natural 
ultraviolet cutoff to the theory. Omitting the (diverging) monopole self-interaction, one 



finds [172] 



dx, / dyi G(x - y) = My 



(3.146) 
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with g = ^QiTiA a combination we also encountered in the two-dimensional charged 
models, |Lz| the (infinite) length of the flux tube, and [ [79[ ] 

My = In f ^) = In(^GL) (3.147) 
Stt \ m-^ J At: 



the line tension. The value kql — l/\/2 for the Ginzburg-Landau parameter ( 3.1 14 ) 
separates the type-11 regime (kgl > 1/V2), where isolated vortices can exist, from 
the type-1 regime (kql < l/'v/Z), where a partial penetration of an external field is 
impossible. Remembering that was the Dirac string, we see from ( 3.146| ) that it 
indeed becomes the core of a vortex in the superconducting phase. 

To identify the regime where vortices are important, we convert the line tension 
My into a magnetic field via 

i/ei = (3.148) 

and compare it with the critical field 

= ^^^5^01, (3.149) 

obtained by equating the tree potential 

Vo = im2|0|2 + iA|^|4 (3.150) 

to —^H^. The critical field physically denotes the smallest value of an external 
field needed to excite vortices. The field H^, on the other hand, is a measure of the 
condensation energy which in turn sets the energy scale. It physically denotes the value 
of an external field at which a type-1 superconductor can no longer resist the magnetic 
pressure and reverts to the normal phase characterized by a perfect penetration of the 
field. The ratio of the two fields 

H^^IH^ (3.151) 

He V2 KGL 

shows that for increasing kgl, vortices become easier to excite. In other words, the 
deeper one enters the type-ll regime, the more important vortex excitations become. 



Since these have been ignored in the calculation of the effective potential (3.112), the 
prediction of a Ist-order transition is only reliable at small kql and breaks down at 
larger values of the Ginzburg-Landau parameter. 

The above construct of inserting a monopole into the system can be used to define a 
so-called disorder parameter Jl7[]. In contrast to the order parameter (3.96 ), a disorder 



parameter should develop an expectation value in the normal, not in the superconduct- 
ing phase. The operator V{L^) describing the monopole with its emerging flux tube is 
easily obtained by noting that in the functional-integral approach, a given field config- 
uration is weigh ted with a Boltzmann factor exp (— Ti^), where the Hamiltonian is 



given by (3.141). From this we infer that the explicit form of the operator is 



^^(Lz) -expjy" [(Vx A).BP-i(BP)']|. (3.152) 
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We are interested in the correlation function where V*{L^) describes 

an additional antimonopole brought into the system at z, with being the accom- 
panying Dirac string running from infinity to z. Since all the integrals involved are 
Gaussian, this expectation value can be evaluated directly. However, we proceed in an 
indirect way to reveal some aspects of the nature of the dual theory and first linearize 
the functional integral over the gauge field by introducing an auxiliary field h. In the 
gauge V • A = 0, which corresponds to setting a = 0, we find 



{V{L^)V*{L^)) = 
j DADhexp 



(3.153) 



ih- (V X A-B^) - '-^A^ 



where now V • B^(x) = (l>o[(5(x— z) — (5(x— z)]. To appreciate the physical relevance 
of the auxiliary field, let us consider its field equation 



h = i(V X A-RP) = ih.. 



(3.154) 



It tells us that apart from a factor i, h can be thought of as representing the local 
magnetic induction h. 

The integral over the vector potential is easily carried out by substituting the field 
equation for A, 



A = X h, 



m , 



back into (3.153), with the result 



{V{L^)V*{L-^)) 



(3.155) 



(3.156) 



Dh exp < — 



1 



-(V X h)2 + 



h 



This shows that magnetic vortices described by a plastic field B^ couple to the fluctu- 
ating massive vector field h, with a coupling constant given hy g — ^qtua = 27r|^| as 
in two dimensions. As the temperature approaches the critical temperature from below, 
(j) tends to zero, so that the vortices decouple from h. The finite penetration depth in 
the superconducting phase is reflected by th e mass term of the h-field ( 3.156 ). 

After carrying out the integral over h in ( 3.156 ), we obtain for the correlation func- 
tion 



p,n(x)G(x-y) p,„(y) 



(3.157) 



.^Br(x)G(x-y)i?,r(y)]|, 



where Pm{x) ~ <l>o[(5(x — z) — (5(x — z)] is the monopole density. The first term in the 
argument of the exponential function contains a diverging monopole self-interaction for 
X = y. This divergence is irrelevant and can be eliminated by defining a renormalized 
operator 

= exp [i$^G(0)] . (3.158) 
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Figure 3.6: Biot-Savart interaction (wiggly line) between two line elements dxi and 
dyi of a magnetic vortex (straight line). 



The second term in the argument is the most important one for our purposes. It repre- 
sents a Biot-Savart interaction between two line elements dxi and di/i of the magnetic 



vortex (see Fig. 3.6). For the renormalized operators we find 



{V,{L^)V;{L^)) = exp(~Afv |L,,|)exp 



2 



V 47r 



(3.159) 



where L^z is the flux tube connecting the monopole at z with the antimonopole at z, 
and |Lzz| is its length. Initially, the two Dirac strings may run to any point at infinity. 
Due to the string tension, however, they join on the shortest path L^z between the 

monopole s. 

The result ( t3.159| ) is central to our line of arguments. It shows thatthe correla- 
tion function {Vr{L2.)V* {L^)) behaves differently in the two phases [170, 173]. In 
the superconducting phase, where because of the Higgs mechanism niA ^ 0, the first 
factor, which amounts to a confining linear potential between the monopole and an- 
timonopole, dominates. As a result, the correlation function decays exponentially for 
distances larger than 1/Afy: 



(3.160) 



This behavior is typical for an operator in a phase without gapless excitations. On 
the other hand, in the high-temperature phase, where niA — 0, the confinement fac- 
tor in the correlation function (3.159) disappears, while the argument of the second 
exponential turns into a pure Coulomb potential. The correlation function remains, 
consequently, finite for large distances: 



{V,iLz)V;iL,)) 



1. 



(3.161) 



By the cluster property of correlation functions this implies that the operator describ- 
ing the finite vortex develops a vacuum expectation value. This signals a proliferation 
of magnetic vortices. Indeed, according to ( 3.147| ) the line tension My of a vortex 
vanishes at the transition point, where — > 0. It should be noted that it is the high- 
temperature phase and not the superconducting phase where K(Lz) develops an ex- 
pectation value. Hence, it is indeed a disorder parameter 

It is interesting to consider the limit tua in (3.156), where the magnetic vortex 
decouples from the massive vector field. This limit yields the constraint V x h = 
which can be solved by setting h = V7. The correlation function then takes the simple 
form 

1 



(k(z)f;(z)> 



D7 exp 



+ i Hp 



(3.162) 
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In the absence of monopoles, the theory reduces to that of a free gapless mode 7 that 
may be thought of as representing the magnetic scalar potential. This follows from 



combining the physical interpretation of the vector field h (3.154) with the equation 
h = V7. Specifically, 

V7==i(V X A-RP). (3.163) 

Inserting the explicit form for the monopole density, which is given by Pm{^) = 
$o[(5(x — z) — (5(x — z)], we see that in terms of the field 7, the correlation func- 
tion reads 

' ' ~ (3.164) 



This demonstr ates th at the operator Vr(iz) describing the fi nite vortex, which was 
introduced in ( 3.152 ) via the singular plastic field ( 3.139 ), is now represented as 
an ordinary field. Since we are in the normal phase, where develops a nonzero 
expectation value, the presence of the phase 7 indicates that this expectation value 
breaks a global U(l) symmetry, with 7 the ensuing Goldstone field. This point will be 
further clarified below. 



We note that Eq. (3.164) reveals also that in the normal phase, the Dirac string 
looses its physical relevance, the right-hand side depending only on the end points z 
and z, not on the path L^z connecting these points. The notion of a magnetic vortex is 
of no relevance in this phase because the vortices proliferate and carry no energy. There 
is also no nontrivial topology to assure their stability. This is the reason for omitting 



any reference to vortex lines in the argument of V in Eqs. ( 3. 162 ) and ( 3. 164 ) 



We are now in a position to write down the dual theory of a three-dimensional 
superconductor It features a grand-canonical ensemble of fluctuating closed magnetic 
vortices of arbitrary shape and length having a repulsive contact interaction. The loop 
gas can be described by a complex disorder j?/'!'* -theory. In addition, as our study of a 
single external vortex revealed, the magnetic vortices couple to the fluctuating vector 
field h w ith a coupUng constant g. Hence, the dual theory is given by [[1741, |l75l |17|, 



163, 1481 



Z = / DhDV'*DV' exp - / H 



with the Hamiltonian 

= :r^(V X h)2 + ih2 + |(v - i^oh-W 



2m\ 



(3.165) 



(3.166) 



where the field ip is minimally coupled to the massive vector field h. Equation (3. 16<: ) 
is a description of the superconducting state in terms of physical variables: the field n 
describes the local magnetic induction, whereas 1/; accounts for the loop gas of mag- 
netic vortices. There are no other physical entities present in a superconductor The 
dual theory has no local gauge symmetry because the vector field h is massive. 



Although ( 3 . 1 65 ) was derived starting from the London limit, it is also relevant near 



the phase transition. The point is that integrating out the size fluctuations of the scalar 
field (f> would only generate higher-order interaction terms. But these modifications do 
not alter the critical behavior of the theory. 
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The line tension My (3.147) appears in the dual theory as a one-loop on-shell mass 



correction to the mass stemming from the graph depicted in Fig. 3.6, which we 
now interpret as a Feynman graph. The straight and wiggly lines represent the ijj- and 
h-field correlation functions, respectively. We have used dashed lines to distinguish 
the Feynman graphs of the dual theory from those in the Ginzburg-Landau model. 

A measure for the interaction strength of a massive vector field in three dimensions 
is given by the dimensionless parameter equal to the square of the coupling constant 
multiplied by the range of the interaction. For the dual theory this factor is /ruA ~ 
mA/e^, which is the inverse of the strength of the electromagnetic gauge field A in the 
superconducting phase. This is a common feature of theories which are dual to each 
other 

Another notable property of the dual theory is that in the limit e ^ it changes 
into a local gauge theory [|l7|]. 



Tit/ 



-(Vxh)2 + |(V-z5h)Vl' 



(3.167) 



as can be checked by rescaling the dual field h in the Hamiltonian (3.166). By tak- 
ing the limit e —> in the Ginzburg-Landau model, we obtain a |0| ^-theory, which 
describes a superfluid, with a decoupled gauge field. 

We next investigate what happens with the dual theory as we approach the critical 
temperature. Remember that (f) and therefore tends to zero as T approaches the 
critical temperature from below. From the first term in the Hamiltonian ( 3.166 ) it again 
follows that V X h — > in this Umit, so that we can write once more h = V7, and 



(3.166) becomes 



(3.168) 



This equation shows that 7, representing the magnetic scalar potential, cannot be dis- 
tinguished from the phase of the disorder field. Indeed, let ^qO be this phase. Then, 
the canonical transformation 9 ^ + j absorbs the scalar potential into the phase of 
Tp; the first term in ( 3.168| ) decouples from the theory and yields a trivial contribution 
to the partition function. In this way, the dual theory reduces to a pure |^|''-theory 



(3.169) 



At the transition temperature, the magnetic vortices proliferate and the field ip develops 
an expectation value. The transition is triggered by a change in sign of In the 



London limit, the Hamiltonian (3.169) then takes the simple form 



(3.170) 



with I -ip I the expectation value of the disorder field, | t/J | / \/2 = {\ip\), and where we now 
represented the phase of ip by $o7 to bring out the fact that 7 describes the magnetic 
scalar potential. From a symmetry point of view [170|, 7 is the Goldstone mode of 
the spontaneously broken global U(l) symmetry of the |?/;|'* -theory. Whereas in the 
Ginzburg-Landau formulation a magnetic vortex is described by a singular plastic field 
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B^, in the dual formulation it is represented by the Noether current associated with the 
global U(l) symmetry, 

jg = V V + 2g2h|Vr (3.171) 

This follows from comparing the terms coupling linearly to the fluctuating h field. In 
the normal phase, where ?/; develops a vacuum expectation value, the Noether current 
becomes in the London limit j = V7. This is the usual relation between the current of 
a spontaneously broken symmetry and the ensuing Goldstone field. 

As we will demonstrate next, ^ has the value \^\ = l/$o |17C] of an inverse flux 
quantum, so that with our normalization choice of the phase of the ^-field, Eq. (3. 17C ) 
takes the canonical form. Let us introduce a closed vortex L in the dual theory (3.17C) 
by minimally coupling the magnetic scalar potential to a vortex gauge field (p^ and 
consider the expectation value of this configuration 



{W{L)) = / D7exp 



1 

We linearize the theory by introducing an auxiliary vector field b via 
1, 



exp 



(3.172) 



(3.173) 



Db exp < — 



ib2+z|^|b- ($oV7~c^f 



This amounts to a duality transformation again. The integral over 7 now yields the 
constraint V • b = 0, demanding b to be the rotation of a vector field, b = V x A. 
This gives 



{W{L)) = yDAcxpj-^ 



(3.174) 



where Ji(x) := 2tt\iJ;\ ^^dj/i(5(x — y) describes the closed vortex. It is natural to 
interpret the fluctuating gapless gauge field A as the magnetic vector potential, and 
the closed vortex as an electric current lo op. T his justifies the use of the symbol W 
in (3.172) which was first introduced in (3.96). (The analogy between vortices and 



electric currents was first pointed out by von Helmholtz.) In this way, we get back the 
Ginzburg-Landau model with just one external electric current loop. The full model 
(3.78) is recovered when we consider a loop gas of these defects and describe it by a 
|(/)|'*-theory, provided we make the identification [170| 



2tt\ 



2e, 



(3.175) 



or I?/; I — l/<i>o- It links the expectation value 'tp of the disorder field ijj describing 
the magnetic vortex loops to the coupling constant 2e of the Ginzburg-Landau model. 
It is the exact analog of the relation between the expectation value cf) of the 0-field 
describing the electric current loops and the coupling constant g — ^arriA of the dual 
theory: 

2n\$\=g. (3.176) 
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Since ip vanishes as T approaches the critical temperature from above, the cou- 
pling constant 2e goes to zero at the critical point and the electric current loops de- 
couple from A. Moreover, Eq. ( |3.111[ ), giving the tension of electric current lines in 
the normal phase, shows that indeed these lines carry no energy in the superconduct- 
ing phase where they proliferate. This solves the apparent contradiction we mentioned 
on page Precisely the same things happened with the magnetic vortices: the cou- 
pling constant g of the dual theory vanishes as the critical temperature is approached 
from below, so that the magnetic vortices decouple from the local magnetic induction 
and proliferate. We thus have a complete duality here between magnetic vortices and 
electric current loops in complementary phases. 

At this stage we can also clarify the physi cal rel evance of the ultraviolet cutoff | 
introduced in the calculation of the tension ( 3.111 ) of a current loop. This calculation 
was performed in the London limit where current lines are considered to be infinitely 
thin. Outside the London limit, vortices of the dual theory have a typical width of the 
order of the coherence length l/|m^|, which therefore provides a natural ultraviolet 
cutoff. 

The line tension ( |3.1 1 1[ ) we previously calculated in the framework of the Ginzburg- 
Landau model can also be obtained in the dual theory directly from (3.172). The in- 
tegration over 7 can be carried out by substituting the field equation of the Goldstone 
field 

1 



7(x) = - 



2$n 



Go(x-y)V.v^(y), 



(3.177) 



where Go is the Green function (3.104) with ijia — 0. This yields an expression 



{W{L))^exp{-^J^J^[W X cpP(x)],Go(x-y) [V x ^.^(y)], \ , (3.178) 



very similar to the expression for the disorder parameter obtained in Eq. (3.157) with 
the monopole density set to zero since 



[V X ¥5^(x)], = 27ry d?/^(5(x - y) 



(3.179) 



This immediately gives 



{W{L)) 



_ ^-Mn-\L\ 



(3.180) 



with \L\ the vortex length and Mw the line tension (3.111) 
The field equation we derive from (3.174), 



V X h = a 



(3.181) 



is one of the two basic equations of magnetostatics. It should be noted, however, that 
an additional factor of i shows up here. As a result, the Biot-Savart law yields opposite 
signs from what is usually the case. Two parallel currents repel instead of attract each 
other so that a single current loop prefers to crumple. It follows that a state where 
these current loops proliferate has zero magnetization — as should be the case in the 
superconducting phase of the Ginzburg-Landau model. 
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Figure 3.7: Two different surfaces S and S' spanning the loop L. The flux through the 
surfaces differ by a factor Att. 



Because of the imaginary current, the local magnetic induction generated by a cur- 
rent loop L is also purely imaginary as follows from Ampere's law, 



h(x) = i — 



dy X 



X y 



2$o 



(3.182) 



where 57 (x) is the solid angle that the loop subtends at x. The same result can also be 
directly derived from the dual theory. Rewriting the field equation for 7 obtained from 
(3.172), we can relate the Goldstone field to the solid angle in the following way Jl7[ ] 



7(x) = 

' 2$o 



.2 (x-y)i 1 ^, . 

x-y-^ 2$o 



(3.183) 



where d is an element of the surface S spanned by the current loop L. Together 
with the observation that the local magnetic induction can, apart from a factor i, be 
identified with the gradient of the phase variable 7 



h = —iS/'-f, 



(3.184) 



[see (3.163)], this yields the previous result (3.182) 



As a side remark we note that the magnetic moment density, or magnetization, 
is represented in the dual theory by tp^. This follows from (3.172) showing that an 
electric current loop couples to the magnetic field V7 via tp^. 

With this physical identification of the Goldstone mode it follows that, in the normal 
state, the disorder parameter 14 (x) essentiaUy measures the solid angle 57(x) 



K-(x) 



(3.185) 



Since the operator V{x.) was constructed by putting a magnetic monopole at x, we 
see that ^51(x) is the magnetic flux emanated by the monopole through the electric 
current loop L. As a last remark we note that one can chose two topologically different 
surfaces spanning the loop L (see Fig. 3.7). Both lead, however, to the same phase 



factor because 51 differs only by a factor of An. 

Because the dual Ginzburg-Landau model involves only physical degrees of free- 
dom and a global U(l) symmetry, it is ideally suited as basis to study the critical prop- 
erties of the superconducting-to-normal phase transition |176|. As has been argued 
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above, at the mean-field level the massive induction field decouples from the theory as 
the critical temperature is approached from below. The resulting theory is a li/)!"* -theory 
which is known to undergo a 2nd-order phase transition with a;j/-exponents. This can 
also be seen by integrating out the massive induction field in the low-temperature phase. 



This only leads to changes in the coefficients of the |?/;| -theory. Explicitly [148], 

= + {ml - g'^) IVP +(u- ^) IV'I^. (3.186) 

In deriving this we used dimensional regularization and (irrelevant) higher-order terms 
are omitted. We note that all contributions stemming from the vector field h vanish 
in the limit where T approaches Tc from below, so that H^^cS reduces to a pure IV'I''- 
theory in this limit. 

A one-loop renormalization-group analysis of the dual theory carried out in Ref. 



[ 1 76 ] also led to the conclusion that the disorder field behaves as in a pure | ^fj | ^-theory 
with xy-exponents. In addition, the exponents of the magnetic induction field where 
shown to retain their Gaussian values because it decouples from the theory. The xy- 



behavior of the heat capacity was established experimentally in [177|, while experi 



ments on YBaCuO [178| seem to confirm the prediction of the Gaussian value for the 
divergence of the magnetic penetration depth. The confirmation relies, however, on a 
delicate finite-size analysis, so that there is still room for a reinterpretation of the data 
and the possibility of suggesting a different critical behavior as has been done by Herbut 



[179] who concluded that the penetration depth has an xy-exponent. One recent lattice 



simulation directly of the Ginzburg-Landau model corroborated the Gaussian value for 



[181]. 



the divergence of the penetration depth [180], while an other found an xy-exponent 
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Chapter 4 

Quantum Phase Transitions 



In this chapter, we discuss various continuous phase transitions at the absolute zero 
of temperature — so-called quantum phase transitions. Unlike in their classical coun- 
terparts taking place at finite temperature and being an equilibrium phenomenon, time 
plays an important role in quantum phase transitions. Put differently, whereas the crit- 
ical behavior of classical 2nd-order phase transitions is governed by thermal fluctu- 
ations, that of 2nd-order quantum transitions is controlled by quantum fluctuations. 
These transitions, which have attracted much attention in recent years (for an introduc- 



tory review, see Ref. [182]), are triggered by varying not the temperature, but some 
other parameter in the system, like the applied magnetic field or the amount of dis- 
order. The natural language to describe these transitions is quantum field theory. In 
addition to a diverging correlation length ^, quantum phase transitions also have a di- 
verging correlation time They indicate, respectively, the distance and time period 
over which the order parameter characterizing the transition fluctuates coherently. The 
way the diverging correlation time relates to the diverging correlation length, 

(4.1) 

defines the so-called dynamic exponent z. It is a measure for the asymmetry between 
the time and space directions and tells us how long it takes for information to prop- 
agate across a distance ^. The traditional scaling theory of classical 2nd-order phase 



transitions is easily extended to include the time dimension [ 1 83 1 because relation (4.1) 
implies the presence of only one independent diverging scale. The critical behavior of 
a phase transition at finite temperature is still controlled by the quantum critical point 
provided T < l/£,t- This is what makes quantum phase transitions experimentally 
accessible. 

We start in the next section discussing the so-called superfluid-to-Mott-insulating 



phase transition in the pure case, while in Sec. 4-. 2 we include (quenched) impurities. 



In Sec. 4.3 we discuss the effective theory describing the fractional quantized Hall 
effect and argue that, in principle, it can be employed to describe the quantum phase 
transitions in quantum Hall systems. In Sec. we then apply renormalization-group 



analysis to this theory. In Sec. 4.5 we discuss scaling and hyperscaling theory applied 
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to the systems under study and in Sec. 4.6 we discuss various experiments probing 
quantum phase transitions. 

4.1 Repulsively Interacting Bosons 

The first quantum phase transition we wish to investigate is the superfluid-to-Mott- 



insulating transition of repulsively interacting bosons in the abs ence o f impurities [ 1 84 1 . 
The transition is described by the nonrelativistic |0|''-theory (1.164), which becomes 



critical at the absolute zero of temperature at some (positive) value /ic of the renor- 
malized chemical potential. The Mott insulating phase is destroyed and makes place 
for the superfluid phase as increases. Whereas in the superfluid phase the single- 
particle (Bogoliubov) spectrum is gapless and the system compressible, the single- 
particle spectrum of the insulating phase has an energy gap and the compressibility k 
vanishes here. 

The nature of the insulating phase can be best understood by putting the theory on 
a lattice. The lattice model is defined by the Hamiltonian 

Hn = -t^(a]aj+i + li.c.) + ^(-A^L^ij + Ufi^j), (4.2) 

3 3 

where the sum J^j is over all lattice sites. The operator creates a boson at site j 
and hj ~ a}-aj is the particle number operator at that site; t is the hopping parameter, 
U the interparticle repulsion, and /il is the chemical potential on the lattice. The zero- 



temperature phase diagram is as follows |184|. In the limit t/U 0, each site is 



occupied by an integer number n of bosons which minimizes the on-site energy (see 



Fig. 4.1) 



e{n) = -fii,n + Un^. (4.3) 

It follows that within the interval 2n— 1 < ii-l/U < 2n+l, each site is occupied by ex- 
actly n bosons. When the chemical potential is negative, n — 0. The intervals become 
smaller when t/U increases. Within such an interval, where the particles are pinned to 
the lattice sites, the single-particle spectrum has an energy gap, and the system is in the 
insulating phase with zero compressibility, k = n^^dn/dfii, = 0. Outside these inter- 
vals, the particles delocalize and can hop through the lattice. Being at zero temperature, 
the delocalized bosons condense in a superfluid state. The single-particle spectrum is 
gapless here and the system compressible (k ^ 0). 

As t/U increases, the gap in the single-particle spectrum as well as the width of the 
intervals decrease and eventually vanish at some critical value tc- For values t > ic of 



the hopping parameter, the superfluid phase is the only phase present (see Fig. 4.2 ). The 



continuum model (1.164), with renormalized chemical potential jj, > describes the 



condensed delocalized lattice bosons which are present when the density deviates from 



integer values (see Fig. 4.1). In the limit /i fi^ from above, the number of delocalized 
bosons decreases and eventually becomes zero at the phase boundary /i = /ic between 
the superfluid and insulating phases. 

Various quantum phase transitions belong to the universality class defined by the 
zero-density transition of repulsively interacting bosons. For example, itinerant quan- 
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3 
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1 



6 f^dU 



Figure 4. 1 : Schematic representation of the average number n of particles per site as 
function of the chemical potential /iL at some finite value of the hopping parameter 

t < ie [|191|]. 




superfluid 



tc/U 



t/U 



Figure 4.2: Schematic representa tion of the phase diagram of the lattice model (4.2) at 
the absolute zero of temperature [ 1 84 ] . 
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turn antiferromagnets [185, 184 ^^'^^ lower-dimensional (clean) supercon- 



ductors belong to this universality class. As we have seen in Sec. 1.4, Cooper pairs 
become tightly bound composite particles in the strong-coupling limit, which are de- 
scribed by the nonrelativistic |0|^-theory with a weak repulsive interaction. For fi > 
^c, the field cj) now describes the condensed delocalized Cooper pairs. When the chemi- 
cal potential decreases, the condensate diminishes, and the system again becomes insu- 
lating for /i < /ic 1 138 1 . By continuity, we expect also the superconductor-to-insulator 
transition of a (clean) weakly interacting BCS superconductor to be in this universal- 
ity class. The restriction to lower dimensions is necessary for two different reasons. 
First, only for d < 2 the penetration depth is sufficiently large [see, for example, be- 
low (^)], so th at it is appropriate to work in the limit Al ^ oo with no fluctuating 
gauge field [188|. Second, in lower dimensions, the energy gap which the fermionic 
excitations face remains finite at the critical point, so that it is appropriate to ignore 
these degrees of freedom. Moreover, since also the coherence length remains finite at 
the critical point, the Cooper pairs look like point particles on the scale of the diverging 
correlation length associated with the phase fluctuations, even in the weak-coupUng 
hmit I pSt . 

The nonrelativistic |0|'*-theory is also of importance for the description of the frac- 



tional quantized Hall effect (FQHE) (see Sec. |4-.3[ ). As function of the applied magnetic 
field, this two-dimensional system undergoes a zero-temperature transition between a 
so-called quantum Hall liquid, where the Hall conductance is quantized in odd frac- 
tions of e^/27r, or, reinstalling Planck's constant, e^/ft,, and an insulating phase. Here, 
the nonrelativistic |(/)|^-theory describes — after coupling to a Chern-Simons term — the 
original electrons contained in the system bound to an odd number of flux quanta. The 
Hall liquid corresponds to the phase with fj, > Hc, while the other phase again describes 
the insulating phase. In this picture, the Hall liquid is characterized by a condensate of 
composite particles. 

It should be noted however that in most of these applications of the nonrelativistic 
1 01 ''-theory mentioned here, impurities play an important role; this will be the subject 
of the succeeding section. 

The critical properties of the zer o-den sity transition of the nonrelativistic \(j>\'^- 
theory were first studied by Uzunov [ 189 |. To facilitate the discussion let us make 
use of the fact that in nonrelativistic theories the mass is — as far as critical phenomena 
concerned — an irrelevant parameter which can be transformed away. This transforma- 
tion changes, however, the scaling dimensions of the (/)-field and the coupling constant 
which is of relevance to the renormalization-group theory. The engineering dimensions 
are 

[x] = -l, [t] = -2, [/io]-2, [\o]=2-d, [0] = id, (4.4) 

with d the number of space dimensions. In two space dimensions the coupling constant 
Ao is dimensionless, showing that the |0|^-term is a marginal operator, and d^ = 2 the 
upper critical space dimension. Uzunov showed that below the upper critical dimension 
there appears a non-Gaussian IR fixed point. He computed the corresponding critical 
exponents to all orders in perturbation theory and showed them to have Gaussian val- 
ues, 1/ = i, z — 2, J] — 0. Here, ly characterizes the divergence of the correlation 
length, z is the dynamic exponent, and 77 is the correlation-function exponent which 
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Figure 4.3: A closed oriented loop. 

X+ XX-XXK 

- xxxx - 

Figure 4.4: Ring diagrams renormalizing the vertex function of the neutral |0|^-theory. 



determines the anomalous dimension of the field cj). The unexpected conclusion that a 
non-Gaussian fixed point has nevertheless Gaussian exponents is rooted in the analytic 
structure of the nonrelativistic propagator at zero (bare) chemical potential (^o = 0): 



G(fc) 



ze 



ikon 



ko 



IT] 



(4.5) 



where, as before, 77 is a small positive constant that has to be taken to zero after the 
loop integrations over the energies have been carried out. By setting jiq = 0, we 
are considering the system at criticality. The rule fco ^ + (4-.5) expresses 
the fact that in this nonrelativistic theory particles propagate only forward in time. In 
Feynman diagrams involving loops with more than one propagator, the integrals over 
the loop energy are convergent and can be evaluated by contour integration with the 
contour closed in either the upper or the lower half plane. If a diagram contains a loop 
which has all its poles in the same half plane, it consequently vanishes. Pictorially, 
such a loop has all its arrows, representing the Green functions contained in the loop, 
oriented in a clockwise or anticlockwise direction [ |190[ (see Fig. 4.3 ). We will refer 
to them as closed oriented loops. Owing to this property most diagrams are zero. 
In particular, all self-energy diagrams vanish. The only surviving ones are the so- 
called ring diagrams which renormalize the vertex (see Fig. 4.4). Because this class of 
diagrams constitute a geometric series, the one-loop result is already exact. The vertex 
renormalization leads to a non-Gaussian fixed point in d < 2, while the vanishing of all 
the self-energy diagrams asserts that the exponents characterizing the transition are not 
affected by quantum fluctuations and retain their Gaussian values [189]. These results 



(glAmstex 



4.1 Repulsively Interacting Bosons 



114 



Table 4. 1 : The upper critical space dimension dc of a nonrelativistic (NR) and a rela- 
tivistic (R) quantum theory with a 10^'" interaction term. 



k 


4(NR) 


4(R) 


2 


2 


3 


3 


1 


2 


oo 





1 



have been confirm ed b y numerical simulations in d = 1 [191| and also by general 



scaling arguments [ |184| | 



We have seen that dc = 2 is the upper critical dimension of the nonrelativistic 
I I ^-theory. Dimensional analysis shows that for an interaction term of the form 



Cint = -ffol-^r" (4.6) 



the upper critical dimension is 



4 = (4.7) 

The two important physical cases are dc = 2, k = 2 and dc = 1, fc = 3, while dc ^ 
when k ^ oo. For space dimensions d > 2 only the quadratic term, is relevant 
so that here the critical behavior is well described by a Gaussian theory. 

In the corresponding relativistic theory, the scaling dimensions of t and x are, of 
course, equal [t] = [x] = — 1 and [0] = ^(d — 1). This leads to different upper critical 
(space) dimensions, viz., 

instead of (p~7|). The two important physical cases are here dc = 3, fc = 2 and dc = 2, 
fc = 3, while dc ^ 1 when k ^ oo. On comparison with the nonrelativistic results, we 
see that the nonrelativistic theory has an upper critical space dimension which is one 



lower than that of the corresponding relativistic theory (see Table 4-.1 ). Heuristically, 
this can be understood by noting that in a nonrelativistic context the time dimension 
cou nts d ouble in that it has a scaling dimension twice that of a space dimension [see 
Eq. ( |4-.4[ )], thereby increasing the effective spacetime dimensionality by one. 

From this analysis it follows that for a given number of space dimensions the crit- 
ical properties of a nonrelativistic theory are unrelated to those of the corresponding 
relativistic extension. 

In closing this section we recall that in a one-dimensional relativistic theory — 
corresponding to the lowest upper critical dimension (dc = 1) — a continuous sym- 
metry cannot be spontaneously broken. However, the theory can nevertheless have a 
phase transition of the so-called Kosterlitz-Thouless type. Given the connection be- 
tween the relativistic and nonrelativistic theories discussed above, it seems interesting 
to study the nonrelativistic theory at zero space dimension (d = 0) to see if a similar 
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rich phenomenon as the Kosterlitz-Thouless transition occurs in the quantum theory. 
This may be of relevance to so-called quantum dots. 



4.2 Including Quenched Impurities 



In the preceding section, we saw that in the absence of impurities repulsively inter- 
acting bosons will undergo a 2nd-order quantum phase transition. As was pointed 
out there, this universality class is of relevance to various condensed-matter systems. 
However, in most of these systems, as well in *He in porous media, impurities play 
an essential if not decisive role. For example, the two-dimensional superconductor-to- 
insulator transition investigated by Hebard and Palaanen [192| is driven by impurities. 
This means that, e.g., the correlation length ^ diverges as | A* — Aj"" when the param- 
eter A characterizing the disorder approaches the critical value A*. Hence, a realistic 
description of the critical behavior of these systems should include impurities. 

Some years ago, it has been argued that upon including quenched impurities, the 
quantum critical behavior of the nonrelativistic |0|*-theory becomes unstable [193|. 
Only after introducing an artificial high-energy cutoff, an IR fixed point was found by 



Weichman and Kim 1 1 94 ] . However, as was pointed out by the authors, such a cutoff is 



difficult to justify as it would imply that time is discrete. So, it is widely accepted that 
the random nonrelativistic |0|^-theory has no perturbatively accessible IR fixed point, 
if any at all 1 184 |. The absence of an IR fixed point in the nonrelativistic 1 01 ^-theory 
theory would imply that also the quantum critical behavior of the systems mentioned 
in the preceding section, is unstable with respect to impurity influences. Because of 
its implications for the description of the critical behavior of these systems we have 
revisited the problem. Below, it will be shown that, contrary to general conviction, the 
random nonrelativistic |0|'*-theory does have a new IR fixed point. The calculations 
are performed without introducing a (physically unnatural) high-energy cutoff. 

To account for impurities, we add to the nonrelativistic |0|^-theory (1.164) a term 



[183] 



/:a-^(x)|0|^ 



(4.9) 



with V'(x) a random field. As we have seen in the preceding section, the theory be- 
comes critical in the limit where the bare chemical potential tends to zero, /io — * 0. We 
shall study the random theory in the symmetrical state where the chemical potential is 
negative and the global U(l) symmetry unbroken. We therefore set /io = —ro again, 
with ro > 0. We leave the number of space dimensions d unspeci fied f or the moment. 
The random field tp(x) is assumed to be Gaussian distributed 1 183]: 



P{^p) — exp 



(4.10) 



characterized by the disorder strength Aq. The engineering dimension of the random 
field is the same as that of the chemical potential which is one, = 1, while that of 
the parameter Ao is [Ao] = 2 — c? so that the exponent in (4.10) is dimensionless. The 
quantity 

Z[ij] ^ I D</>*D0 exp (if £ ) , (4.11) 
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where C now stands for the Lagrangian (1.164) with the term (4.9) added, is the zero- 
temperature partition function for a given impurity configuration In the case of 
quenched impurities, the average of an observable 0{(f)* , 4>) is obtained as follows 



(0(r»)= / DV;P(V)(0(<^*»)^, 



(4.12) 



where {0{(t>* , <j>))ip indicates the g rand- canonical average for a given impurity config- 
uration, i.e., taken with respect to ( 4.11 ). In other words, first the ensemble average is 
taken, and only after that the averaging over the random field is carried out. 

Since ^(x) depends only on the d spatial dimensions, the impurities it describes 
should be considered as grains randomly distributed in space. When — as is required 
for the study of quantum critical phenomena — time is included, the static grains trace 
out straight w orldli nes. That is to say, the impurities are line-like. It has been shown by 
Dorogovtsev 1 195 1 that the critical properties of systems with extended defects must 
be studied in a double epsilon expansion, otherwise no IR fixed point is found. The 
method differs from the usual epsilon expansion, in that it also includes an expansion 
in the defect dimensionality ed. To carry out this program in the present context, where 
the defect dimensionality is determined by the dimensionality of time, the theory has 
to be formulated in ed time dimensions. The case of interest is ed = 1, while in the 
opposite limit, ed 0, the random nonrelativistic |0| ''-theory reduces to the classical 
spin model with random (point-like) impurities. Hence, ed is a parameter with which 
quantum fluctuations can be suppressed. An expansion in ed is a way to perturbatively 
include the effect of quantum fluctuations on the critical behavior. Ultimately, we will 
be interested in the case ed = 1. 

To calculate the quantum critical properties of the random theory, which have been 
first studied in [193], we will not employ the replica method [196], but instead follow 
Lubensky 



[ |197[ ]. In this approach, the averaging over impurities is carried out for 
each Feynman diagram separately. The upshot is that only those diagrams are to be 
included which remain connected when Ap , the parameter characterizing the Gaussian 
distribution of the impurities, is set to zero [ 198 ] . To obtain the relevant Feynma n rule s 
of the random theory we average the interaction term (4.9) over the distribution (4.10): 



Dij; P{ip) exp 



exp 



4* 



Ao J d^Hd^H' d'^x x)\^\(l){t', x) 



(4.13) 



The randomness is seen to result in a quartic interaction term which is nonlocal in time. 



The factor i'^'^ appearing in (4.13) arises from the presence of ed time dimensions, each 



of which is accompanied by a factor of i. The Feynman rules of the random theory are 
now easily obtained 



_j-«dgi(i^l+'^2H l-Wjj)77 



UJi + UJ2 + ■ ■ ■ + (^e 



k-^ — ro + irj 



(4.14) 
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where we note that the Lagrangian in time dimensions involves instead of ju st one 
time derivative, a sum of ed derivatives: dt ^ dti+dt2 + - ■ ■ + dt^^ ■ The integral ( 4.17 ) 
has an additional convergence factor exp(zci;?7) for each of the energy integrals. This 
factor, which is typical for nonrelativistic quantum theories is to be included in 
self-energy diagrams containing onl y on e (/)-propagator. 

Following Weichman and Kim 1 194], we evaluate the integrals over loop energies 
assuming that all energies are either positive or negati ve. T his allows us to employ 
Schwinger's propertime representation of propagators 1146], which is based on the 
integral representation of the gamma function, 

— = TTTT / —r^e-^". (4.15) 

The energy integrals we encounter to the one-loop order can be carried out with the 
help of the equations 



1 



(27r)'=d LUi+LU2^ 1-' 



X zt irj 



r(i-£d) 

(27r)^d 



sgn(a;)|a;|''' ^ (^i 



±z sgn(x)7r€d 



(27r)^i wi + u;2 H h 



1X1] 



sin(i7red) 



sgn(x) 



(4.16) 



(4.17) 



where rj is again an infinitesimal positive constant which is to be taken to zero after 
the energy integrals have been carried out. The prime on the integrals is to remind 
the reader that the energy integrals are taken over only two domains with either all 
energies positive or negative . The energy integrals have been carried out by using again 
the integral representation ( 4.15 ) of the gamma function. In doing so, the integrals 
are regularized and — as is always the case with analytic regularizations — irrelevant 
divergences suppressed. 



By differentiation with respect to x, Eq. (4.16) can, for example, be employed to 
calculate integrals involving integrands of the form 1 / {uji+llJ2 + ■ ■ ■+^ea — x+irj)'^. Is 
is easily checked that in the limit 1, where the energy integral can be performed 
with help of contour integration, Eqs. (4.16) and (4.17) reproduce the right results. 
When considering the limit of zero time dimensions (ed 0), it should be remembered 
that the energy integrals were taken over two separate domains with all energies either 
positive or negative. Each of these domains is contracted to a single point in the limit 
ed 0, so that one obtains a result which is twice that obtained by simply purging any 
reference to the time dimensions. 

Before studying the random theory, let us briefly return to the repulsively inter- 
acting bosons in the absence of impurities. In this case, there is no need for an ed- 
expansion and the formalism outlined above should yield results for arbitrary time 
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dimensions < ed < 1, interpolating between the classical and quantum limit. After 



the energy integrals have been performed with the help of Eqs. (4.16) and (4-. 17), the 



standard technique of integrating out a momentum shell can be applied to obtain the 
renormalization-group equations. For the correlation-length exponent v we obtain in 



this way [201] 



1 



1 



2 (to + 4) - (m + 3)ed 



cos^(^7red) 



(4.18) 



Here, e = 4 — 2ed — d is the deviation of the effective spacetime dimensionality from 
4, where it should be noted that in (canonical) nonrelativistic theories, time dimen- 
sions have an engineering dimension twice that of space dimensions. (This property is 
brought out by the Gaussian valu e z = 2 for the dynamic exponent z.) For compari- 
son we have extended the theory ( 1.164 ) to include to complex fields instead of just 
one field. In the classical limit, Eq. ( 4.18 ) give s the well-known one-loop result for a 
classical spin model with 2to real components |ll82 



1 / e TO + 1 

2 \ ^ 2TO + 4 



(4.19) 



while in the quantum limit it gives the result ^ i, as required. 

The exponent (4.18), and also the location of the fixed point, diverges when the 
number of time dimensions becomes ed (m + 4)/(m + 3). Since this value is 
always larger than one, the singularity is outside the physical domain < ed < 1. This 
simple example illustrates the viability of the formalism developed here to generate 
results interpolating between the classical and quantum limit. 

We continue with the random theory. After the energy integrals have been carried 
out, it is again straightforward to derive the renormalization-group equations by inte- 
grating out a momentum shell A/6 < fc < A, where A is a high-momentum cutoff and 
b — exp(Z), with I infinitesimal. Defining the dimensionless variables 



A 



(27r)^< 



■AA" 



(4.20) 



where Kd given in ( 1.249 ) is the area of a unit sphere in d spatial dimensions divided 
by (27r)'', we find J20l|] 



^ = eA-8[r(l-ed) + (TO + 3)r(2-ed)]cos(i7red)A2+6AA 

dA 
"dT 



(e + 2ed)A + 4A2 - 16(m + l)r(2 - Cd) cos(i7red)AA 



df TO -I- 1 cos^(i7red) c ; 

- - 2f + A7T . \^ ^ A-A. 

dl r(ed) smUned) 



(4.21) 



These results are to be trusted only for small values of ed. For illustrative purposes we 
have, however, kept the full ed dependence. The set of equations yields the fixed point 



A* 



1 



e + 6ed 



16cos(i7red)r(l - ed) 2m(l - ed) - 1 



(4.22) 
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A* = 



1 m(l - ed)(2ed - e) + 2ed(4 - 3ed) + ^(2 - e^) 
4 2m(l-ed)-l 



and the critical exponent 



1 

i> = — 
2 



+ 2ed TO + 1 (6ed + e) [cd + cos(7red)] 



16 



16 



2to(1 - ed) - 1 



(4.23) 



The dynamic exponent is given by z — 2 + A*. We see that both A* and A* diverge 
when ed ^ 1 — l/2m. At this point, the fixed point becomes unphysical. The singular- 
ity separates the quantum regime ed ~ 1 from the classical regime ed ~ about which 
perturbation theory is to be carried out. When the equations are expanded to first or- 
der in ed, we recover the IR fixed point found by Weichman and Kim 1 194] using an 
high-energy cutoff: 



A* 



1 e + 6ed 
16 2to - 1 ^ 



A* 



4 



2(to + 4)ed 



2to - 1 



with the critical exponent 



v = — 
2 



1 + - 



1 3me + (5to + 2)2ed 



2to - 1 



(4.24) 



(4.25) 



We thereby provide support for the existence of this fixed point. 

The value of the critical exponent ( 4.25 ) should be compared with that of the classi- 
cal spin model with 2m components in the presence of random impurities of dimension 



1 + - 



1 3TOe + (5to + 2)ed 



2to- 1 



(4.26) 



Taking into account that in a nonrelativistic quantum theory, time dimensions count 
double as compared to space dimensions, we see that both results are equivalent. As 
to the dynamic exponent, we mention that a remarkable simple scaling argument (see 
Sec. 4.5) using the quad ratic terms in the effective theory ( 1.177| ) predicts the exact 
value z = d for ed = 1 [184|. The perturbative result z = 2 + A*, with A* given by 



( |4.24| ), is seen to be far away from this. 

The limit of interest to us, corresponding to ed — 1, is probably difficult to reach 
by low-order perturbation theory for the quantum regime is separated by a singularity 
from the classical regime where perturbation theory applies. Although this might be an 
artifact of the one-loop calculation, it is unlikely that by including a few more loops, 
the quantum regime becomes accessible via the classical regime. We note that the 
singularity moves towards ed = 1 when the number of field components increases 

If the IR fixed point is to be of relevance to the systems mentioned in Sec. 
the impurities have to lead to localization in the superfluid phase. For the model at 



4.1 



hand, albeit in d = 3, this connection has been established by Huang and Meng [199|. 
They showed that in this phase, the impurities give rise to an additional depletion of 
the condensate as well as of the superfluid mass density. (The latter is defined, as we 
discussed in Sec. 1.7, by the response of the system to an externally imposed velocity 
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field as specified by the expression for the momentum density, or mass current g.) 
They found that the depletion of the superfluid mass density is larger than that of the 
condensate, indicating that part of the zero-momentum states belongs to the normal 
fluid rather than to the condensate. They interpreted this as implying that these states 
are trapped by impurities. 

This situation should be contrasted to the one in the absence of impurities, where 
the condensate is depleted only due to the interparticle repulsion. Despite the deple- 



tion, all the particles were nevertheless seen in Sec. 1 .7 to participate in the superflow 
motion at zero temperature. In other words, the normal fluid is dragged along by the 
condensate. This conclusion was based on the observation that the momentum density 
was given by g = mnvg, where n is the total particle number density and Vg the ve- 
locity with which the condensate moves, i.e., the superfluid velocity. (For clarity, we 
have reintroduced the mass parameter m here.) 



In terms of the shifted field ( 1.166| ), the random term (O) reads 



£a = ^(x)(|0p + |0P + + r4>)- (4.27) 

The first two terms lead to an irrelevant change in the chemical potential, so that we 
only have to consider the last two terms, which we can cast in the form 

£a = V'(x)^^$, * = ( ^* ) ■ (4.28) 

The integral over $ is Gaussian in the Bogoliubov approximation and therefore easily 
performed to yield an additional term to the effective action 

Sa = -\( V(x)l'^Go(x-y)l>^(y), (4.29) 

where the propagator Go is the inverse of the matrix A/q introduced in ( |1.170| ) with the 
field U{x) set to zero. Let us first Fourier transform the fields. 



Go{x-y) - / e-''= (^-«) Go(fc) (4.30) 

Jk 

V'(x) = j e*'^-''?/'(k). (4.31) 
The contribution to the effective action then appears in the form 

^A = -^^|^(k)P$^G(0,k)#. (4.32) 



Since the random field is Gaussian distributed [see ( 4.10| )], the average over this field 
representing quenched impurities yields, 

(IVXk)^ ^ lyAo. (4.33) 



The remaining integral over the loop momentum in (4.32) is readily carried out to yield 
in arbitrary space dimensions 

{C^)^ T{l-d/2)(—) |0p(6Ao|0p-/.o)'*/2-iAo. (4.34) 



2 ' ' ' V27r/ 
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Because this is a one-loop result, we may to this order replace the bare parameters with 
the (one-loop) renormalized ones. 



In Sec. 1.6 we saw that due to the interparticle repulsion, not all the particles reside 



in the condensate. We expect that the random field causes an additional depletion of 
the condensate. To obtain this, we differentiate (4.34) with respect to the chemical 
potential. This gives [poojpoil] 



9/1 



-d/2 



2-2-d/2-l 

'-o 



Ao, 



(4.35) 



where uq denotes the number density of particles residing in the condensate. The 
divergence in the limit A ^ for d < 4 signals the collapse of the system when the 
interparticle repulsion is removed. 

We next calculate the mass current g to determine superfluid mass density, i.e., the 
mass density flowing with the superfluid velocity Vg. As we have seen in the preceding 
section, in the absence of impurities and at zero temperature all the particles participate 
in the superflow and move on the average with the velocity Vg. We expect this no longer 
to hold in the presence of impurities. To determine the change in the superfluid mass 



density due to impurities, we replace /io with /ics as defined in (1.236) and ido with 
ido— (u—Vs) ■ (— iV) in the contribution ( 4. 32 ) to the effective action, and differentiate 
it with respect to the externally imposed velocity, — u. We find to linear order in the 
difference u — Vg: 

g = PsVs + PnU, (4.36) 



with the superfluid and normal mass density [201 1 



m \ n riA 

d 



Pn 



-mn^. 



(4.37) 



We see that the normal density is a factor i/d larger than the mass density mn^ 
knocked out of the condensate by the impurities. (For d = 3 this gives the factor | 



first found in Ref. [199].) Apparently, part of the zero-momentum states belongs for 
c? < 4 not to the condensate, but to the normal fluid. Being trapped by the impuri- 
ties, this fraction of the zero-momentum states are localized. This shows that the phe- 
nomenon of localization can be accounted for in the Bogoliubov theory of superfluidity 
by including a random field. 

An other realistic modification of the nonrelativistic |0|*-theory is to include a l/r 
Coulomb repulsion. Using very general scaling arguments (see Sec. 4.5) in the context 
of the effective theory. Fisher, Grinstein, and Girvin [ 188 1 predicted that in the presence 
of impurities the inclusion of this interaction changes the value of the dynamic expo- 
nent to z = 1, which is exact again. This prediction has been confirmed in experiments 
on the superconductor-to-insulator transition in two-dimensional films [ 1 92 ] . The same 
value z = 1 has also been found in the quantum Hall transitions, and — more recently — 
in metal-to-insulating transitions observed in dilute two-dimensional electron systems 
in silicon MOSFETS [ |20^ p03[ ^04] , p05| ] . This suggests that the 1 /r Coulomb interac- 
tion plays an important role in these systems too. However, let us at this stage mention 



the renormalization-group results obtained by Giamarchi and Schulz [206] who stud 



ied fermions with short-range interactions in a random one-dimensional system. They 
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found first of all that the superconductor-to-insulator transition these fermions undergo 
is indeed in the universality class of repulsively interacting bosons in the presence of 
impurities. Moreover, their result for the conductivity is consistent with the value z — 1 
[ ]188| ], implying that already a local repulsive interaction as used in the nonrelativistic 
I I* -theory can lead to this value for the dynamic exponent. 



4.3 CSGL Theory 

We now turn to the fractional quantized Hall effect (FQHE) which is the hallmark of 
a new, intrinsically two-dimensional condensed-matter state — the quantum Hall liq- 
uid. Many aspects of this state are well understood in the framework of the quantum- 



mechanical picture developed by Laughlin [207]. Considerable effort has nevertheless 



been invested in formulating an effective field theory which captures the essential low- 
energy, small-momentum features of the liquid. A similar approach in the context 
of superconductors has proven most successful. Initially, only the phenomenological 
model proposed by Ginzburg and Landau [ [75| ] in 1950 was known here. Most of the 
fundamental properties of the superconducting state such as superconductivity — the 
property that gave this condensed-matter state its name, Meissner effect, magnetic flux 
quantization, Abrikosov flux lattice, and Josephson effect, can be explained by the 
model. The microscopic theory was given almost a decade later by Bardeen, Cooper, 
and Schrieffer [|^. Shortly here after, Gorkov made the connection between the 
two approaches by deriving the Ginzburg-Landau model from the microscopic BCS 
theory, thus giving the phenomenological model the status of an effective field theory. 
A first step towards an effective field theory of the quantum Hall liquid was taken 



by Girvin and MacDonald [208] and has been developed further by Zhang, Hansson 
and Kivelson [209], who also gave an explicit construction starting from a microscopic 
Hamiltonian. Their formulation incorporates time dependence which is important for 
the study of quantum phase transitions. This approach has proven very successful 



(for a review see Ref. [210|). In this and the following section, we shall argue that 
the effective theory — the so-called Chern-Simons-Ginzburg-Landau (CSGL) theory — 
can also be employed to describe the quantum phase transition which a quantum Hall 
system undergoes as the applied magnetic field changes. 

We shall in this section first recall some basic properties of the CSGL theory and 
show how it can be used to describe the field-induced Hall-liquid-to-insulator transition 



of a Hall liquid. In Sec. 4.4 we then apply renormalization group theory to the CSGL 
theory and study its critical properties. 

An important ingredient for obtainin g an effective theory of the FQHE was the 
identification by Girvin and MacDonald [ [208[ | of a bosonic operator (j> which exhibits 
(algebraic) off-diagonal long-range order The long-range order was found to be of a 
type known to exist in two-dimensional bosonic superfluids. They argued that this field 
should be viewed as an order parameter in terms of which the effective field theory 
should be formulated. To account for the incompressibility of the quantum Hall liquid 
they suggested to couple to a statistical gauge field a^. The gapless spectrum of 



the neutral system then changes into one with an energy gap [209], thus rendering the 
charged system incompressible. 
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Girvin and MacDonald assumed that the statistical gauge field is governed solely 
by a Chern-Simons term 

Q 

with V X a the statistical magnetic field and 6 a constant. Because of the absence of 
a kinetic term (the usual Maxwell term), the statistical gauge field does not represent 
a physical degree of freedom. In a relativistic setting, a Maxwell term is usually gen- 
erated by quantum corrections so that the statistical gauge field becomes dynamical 
at the quantum level. The quantum theory then differs qualitatively from the classical 
theory. On the other hand, as we shall see below, this need not be the case in a nonrel- 
ativistic setting. That is to say, the Ansatz of the absence of a Maxwell term is here not 
necessarily obstructed by quantum corrections. 



The CSGL theory is described by the Lagrangian [ 209 1 



C = i^*Do^ - :^|D0|2 + ^o|0|' - Ao|0|" + ^cs- (4.39) 

The covariant derivatives Dq = do + ieAo + ieao and D = V — ieA — iea give 
a minimal coupling to the applied magnetic and electric field described by the gauge 
field Af^ and also to the statistical gauge field. For definiteness we will assume that our 
two-dimensional sample is perpendicular to the applied magnetic field, defining the z- 
direction, and we choose the electric field to point in the x-direction. The charged field 
(f> represents the Girvin-MacDonald order parameter describing the original electrons 
bound to an odd number of flux quanta. To see that it indeed does, let us consider the 
field equation for qq: 

^ -e6iV X a. (4.40) 
The simplest solution of the CSGL Lagrangian is the uniform mean-field solution 

= a = -A, ao = -^o = 0, (4.41) 

where n indicates the constant fermion number density. The statistical gauge field is 
seen to precisely cancel the applied field. The constraint equation ( 4.4C| ) then becomes 



n = eOH, (4.42) 

with H the applied magnetic field. Now, if we choose 9^^ — 2tt{21 + 1), it follows on 
integrating this equation that, as required, with every electron there is associated 21 + 1 
flux quanta: 

where = ^/ *&o> with ^ = J^H the magnetic flux, indicates the number of flux 



quanta. Equation ( 4.42 ) implies an odd-denominator filling factor vh which is defined 
by 
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The coupling constant Aq (> 0) in ( [4-.39 ) is the strength of the repulsive contact 



interaction between the composite particles, and /io is a chemical potential introduced 
to account for a finite number density of composite particles. 

It is well known from anyon physics that the inclusion of the Chern-Simons term 



changes the statistics of the field (f> to which the statistical gauge field is coupled [211 1. 
If one composite particle circles another, it picks up an additional Aharonov-Bohm fac- 
tor, representing the change in statistics. The binding of an odd number of flux quanta 
changes the fermionic character of the electrons into a bosonic one for the composite 
particles, allowing them to Bose condense. The algebraic off-diagonal long-range or- 
der of a quantum Hall liquid can in this picture be understood as resulting from this 
condensation. Conversely, a flux quantum carries 1/(2/ + l)th of an electron's charge 



IgOj, and also 1/(2/ + l)th of an electron's statistics [ |212| ]. 



The defining phenomenological properties of a quantum Hall liquid are easily 



shown to be described by the CSGL theory [209, ElOr. From the lowest-order ex 



pression for the induced electromagnetic current one finds 

eJ^ = 1^ = = ^ = -e^Oe^jOoa, - d,ao) = eHe,,Ej, (4.45) 



with E the applied electric field and where the Lagrangian (4.39) is written as a sum 
C = + -Ccs- It follows that the Hall conductance <Jxy is quantized in odd fractions 
of e^/27r, or, reinstalling Planck's constant, /h. This result can also be understand in 
an intuitive way as follows. Since the composite particles carry a charge e, the appUed 
electric field gives rise to an electric current 

I = e— (4.46) 
at 

in the direction of E, i.e., the x-direction. This is not the end of the story because the 
composite objects carry in addition to electric charge also 21 + \ flux quanta. When the 
Goldstone field if encircles 21 + 1 flux quanta, it picks up a factor 27r for each of them 

V-(^ = 27r(2Z + l). (4.47) 

r 

Now, consider two points across the sample from each other. Let the phase of these 
points initially be equal. As a composite particle moves downstream, and crosses the 
line connecting the two points, the relative phase A</j between them changes by 2tt{21+ 
1). This phase slippage [^ leads to a voltage drop across the sample given by 

yH = -<9oA(^=(2/ + l)$o^, (4.48) 
e at 

where the first equation can be understood by recalling that due to minimal coupling 
do(p dof + cAq. For the Hall resistance we thus obtain the expected value 

P.y='^ = (2l + 1)^. (4.49) 

If the CSGL theory is to describe an incompressible liquid, the spectrum of the 
single-particle excitations must have a gap. Without the coupling to the statistical gauge 
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field, the spectrum is given by the gapless Bogoliubov spectrum (1.174). To obtain the 



single-particle spectrum of the coupled theory, we integrate out the statistical gauge 



field. The integration over oq was shown to yield the constraint (4.40) which in the 
Coulomb gauge V • a = is solved by 

a, = le,A|0|^ (4.50) 

The integration over the remaining components of the statistical gauge field is now 
simply performed by substituting (4.50|) back into the Lagrangian. The only nonzero 



contribution arises from the term — e^|(/>pa^/2TO. The spectrum of the charged system 
acquires as a result an energy gap cUc 



E{k) = ^iu^ + e^k) + 2fioeik), (4.51) 

with LUc = iJio/Wm\Q. To lowest order, the gap equals the cyclotron frequency of a 
free charge e in a magnetic field H 

n eH 

c^c = 7— = — . (4.52) 

The presence of this energy gap results in dissipationless flow with a^x = 0. 

These facts show that the CSGL theory captures the essentials of a quantum Hall 
liquid. Given this success, it is tempting to investigate if the theory can also be em- 
ployed to describe the field-induced Hall-liquid-to-insulator transitions. It should how- 
ever be borne in mind that both the l/|x|-Coulomb potential as well as impurities 
should be incorporated into the theory in order to obtain a realistic description of the 
FQHE. The repulsive Coulomb potential is believed to play a decisive role in the for- 
mation of the the composite particles, while the impurities are responsible for the width 
of the Hall plateaus. As the magnetic field moves away from the magic filling factor, 
magnetic vortices will materialize in the system to make up the difference between 
the applied field and the magic field value. In the presence of impurities, these de- 
fects get pinned and do not contribute to the resistivities, so that both a^x and (Txy are 
unchanged. Only if the difference becomes too large, the system reverts to an other 
quantum Hall state with a different filling factor 

For positive bare chemical potential /ig > the Girvin-MacDonald order parame- 
ter (/) has a nonvanishing expectation value given by |0p = /io/2Ao, implying that the 
composite particles are condensed. When /io 0, the condensate is drained of com- 
posite particles, and at /iq — 0, it vanishes altogether The system becomes critical here 
and reverts to the insulating phase characterized by a negative bare chemical potential. 

In the spirit of Landau, we take a phenomenological approach towards the field- 
induced phase transition of the CSGL theory. And assume that when the applied mag- 
netic field H is close to the upper critical field iJ at which the quantum Hall liquid 
with filling factor i/ff is destroyed, the chemical potential of the composite particles 
depends linearly on H, i.e., /in ~ ~ ^H. This state can of course also be de- 

stroyed by lowering the applied field. If the system is near the lower critical field H^^^ , 
we assume that the chemical potential is instead given by /iq ~ eH — eH:j^ . This is 
the basic postulate of our approach. 
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In the CSGL Lagrangian (4-. 39) we have again transformed the mass m of the non- 
relativistic j^j^'-theory away. In addition to the engineering dimensions (4.4), we have 
for the Chern-Simons field 



[eat] ~ 1, [eao] = 2, [0] = 0. 



(4.53) 



In two space dimensions the coupling constant Xq was seen to be dimensionless, im- 
plying that the |(/)|'*-term is a marginal operator. We see from (4.53) that also the Chern- 
Simons term is a marginal operator. Hence, the CSGL theory contains precisely those 
terms relevant to the description of the Hall-liquid-to-insulator transition in a quantized 
Hall system. 



4.4 Renormalization Group 

We now turn our attention to the critical properties of the CSGL theory. A remaining 
open problem in this context is the nature of the Hall-liquid-to-insulating transitions 
in the theory. Taking place at the absolute zero of temperature, these transitions are 
a pure quantum phenomenon. They are driven either by impurities, a periodic poten- 
tial, or by the applied magnetic field. The impurity-driven phase transition has been 



considered in Ref. [213], while the transition driven by a periodic potential has been 



studied in Refs. [214, 215]. The latter problem was mapped onto a relativistic CSGL 



theory and it was concluded [ 215 1 that the transition induced by a periodic potential is 
Ist-order We in this section shall discuss the transition induced by the magnetic field 
in the absence of both impurities and a periodic potential. We will be working with the 
standard nonrelativistic CSGL theory. As has been eluded to above, there are essential 
differences between relativistic and nonrelativistic Chern-Simons theories. There are 
features special to the nonrelativistic theory which in our view are important for under- 
standing the field-induced transition, and which cannot be implemented in a relativistic 
formulation. 

Experimentally, if the external field is changed so that the filling factor vu moves 
away from an odd-denominator value, the system eventually becomes critical and un- 
dergoes a transition to an insulating phase. We will argue that this feature is encoded 
in the CSGL theory. We shall show that the theory has an IR fixed point, implying that 
the phase transition is 2nd-order, and determine the critical exponents. They turn out 
to be universal and independent of the filling fraction. 

We will be working in the Coulomb gauge. To this end we add a gauge-fixing term 



(V • a)^/2agf to the Lagrangian (4.39), and take the limit agf 0. In this gauge, the 



propagator of the statistical gauge field is purely off-diagonal 



kj 



-wwwwv = Go^ik^ = -G,o(fc) = iae^jri, (4.54) 
i 



with a — 1/9. The two propagators (4.5) and (4.54) of the CSGL theory are causing 
infrared divergences. We regularize these by introducing an infrared cutoff via the sub- 
stitution ^ + 2ro in both propagators. In the scalar propagator, this amounts to 
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-i\ -ie -ie^Sj^j «e(pj+p^) 



Figure 4.5: Feynman rules of the CSGL theory. 




Figure 4.6: A closed oriented loop containing in addition a gauge field propagator. 



taking a negative chemical potential, /iQ — — rp, so that we are in the symmetrical state 
of the theory. Because we consider the theory in its upper critical dimension and at crit- 
icality, the parameter rg is the only dimensional parameter present. It is necessary to 
introduce a dimensional parameter in a theory not having any, not just to avoid infrared 
divergences, but also to be able to carry out the renormalization-group program — this 
parameter playing the role of renormalization-group scale parameter. There are various 
ways to do so, our proposal is one possibility. Another would be to take the renormal- 
ization point at nonzero external momenta, as has been done by, for example, Bergman 



and Lozano [ [216[ |. The critical properties are, of course, independent of the particular 
choice. 

We evaluate Feynman diagrams according to the following computational scheme. 
First, we carry out the discrete sums over repeated indices as well as the integrations 
over loop energies. Then, to handle ultraviolet divergences, the resulting momentum in- 
tegrals are analytically continued to arbitrary space dimensions d. It is to be noted that 
this step is taken after the discrete sums have been carried out and the antisymmetric 
tensor — which makes sense only in d = 2 — has disappeared from the expressions. 
Finally, the counter terms and the critical exponents are determined by employing the 
minimal subtraction scheme in which only those parts of the counter terms are consid- 
ered which diverge in the limit d ^ 2. The Feynman rules for the vertices are given in 



Fig. 4.5 



The reason we have chosen the Coulomb gauge is that in this gauge t he sta tistical 



gauge field propagator does not depend on the energy variable ko [see Eq. ( )4.54[ )]. This 
means that the analytic structure of a given diagram is determined solely by its scalar 
propagators. As a result, the rule that closed oriented loops vanish is not corrupted 



if a loop contains in addition gauge field propagators. (See Fig. |4.6| for an example.) 
Hence, also in the coupled theory most diagrams vanish, making the quantum critical 
behavior of the CSGL theory tractable. 



©Amstex 



4.4 Renormalization Group 



128 



Figure 4.7: A one-loop and a two-loop example of vanishing diagrams which would 
otherwise renormalize the charge e. 



In particular, all the diagrams which would renormalize the charge e vanish. (A 
one-loop and a two-loop example of this class of diagrams are depicted in Fig. 4.7 ) Put 
differently, no Maxwell term for the statistical gauge field is generated at the quantum 
level so that the statistical gauge field propagator (4-.54) remains purely off-diagonal. 
This is an important characteristic of the nonrelativistic CSGL theory not shared by 
the corresponding relativistic extension where a Maxwell term is automatically gener- 
ated at the quantum level. Since the Chern-Simons term is also not renormalized, the 
statistical gauge field propagator (4.54) is not effected at all by quantum corrections. 



As in the neutral case, the boson self-energy E also vanishes at every loop order in 
the charged theory, and the only object that renormalizes is the self-coupling parameter 
Xq. Consequently, if an IR fixed point is found in the CSGL theory, the 2nd-order phase 
transition described by it has Gaussian exponents. 

In order to examine the presence of an IR fixed point, we compute the beta function 
using the scheme outlined in [217|. To this end we evaluate the diagrams: 



1 



2r 



h"^ — r + irj ko — ^k^ — r + irj 
Id, 



with A and r denoting the renormalized parameters, 

r(i-fi/2) 1 



Id 



(4.55) 



(4.56) 



and 



k2 + 2r _ ik2 - r + irj k2 + 2t 



" ' '"U 2 

k2 



k (k2 + 2r)2 



where we introduced the abbreviation 



(4.57) 



^d 



(4.58) 
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XXK 



+ 4 




+ 2 



Figure 4.8: The relevant two-loop diagrams. 



There are three other diagrams of the form ( 4.57 ), one having the right external legs 
crosse d and two having the triangle turned upside-down. To account for thes e, the 
result ( 4.57 ) must be included with an additional factor of 4. The diagrams ( 4.55 ) 
and ( 4.57 ) are found to diverge in the limit d ^ 2. The divergences are cancelled by 
including the one-loop counter term 



X 



(4.59) 



with e = 2 — d the deviation from the upper critical dimension and VI2 — ^ /47r. This 
leads in the usual way to the one-loop beta function in the upper critical dimension 
(c? = 2 in our case) 



1 - A 



d_ 
'dX 



a 



d_ 

da 



1 



(4A2 



(4.60) 



where ai is {i times) the residue of the simple pole in the counter term (4.59). The 



result is in accordance with previous studies by Lozano and Bergman 1216| who used 
a momentum cutoff to regularize the ultraviolet divergences. 

We continue with the two-loop calculation. The relevant diagrams are given in Fig. 
4.8. The factor 4 in the second diagram is to account for the related diagram having 
the 4-vertex to the right rather than to the left and for the diagrams having the triangles 
turned upside-down. The factor 2 in the last diagram accoun ts for the fact that there 
are two vertices which can be replaced with the counter term ( 4.59 ). The first diagram 
is readily shown to yield 



XXK 



(4.61) 



whereas the second gives 




-2iXa' 



1 (k-l).l 



k.l k2 + 2r 12 + 2r (k - 1)2 + 2r ' 



(4.62) 



To evaluate the remaining integrals, we employ the identity 
f k If 1 



/k (k2 + 2r) [(k - 1)2 + 2r] 2 (k2 + 2r) [(k - 1)2 + 2r] 



(4.63) 
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which can be proven by introducing a new integration variable k' = — k + 1 in the 
left-hand side. In this way, we obtain for the second diagram 




= iXa" (I J 



Jd) 



(4.64) 



where Jd is the integral 



k,i k2 + 2r P + 2r (k - 1)2 



2r 



(4.65) 



whi ch re mains finite in the limit d — > 2. Finally, t he las t diagram in Fig. 4.8 is given 
by (4.55) with —iX replaced with the counter term (4.59). 



When the results are added, the sum of the three diagrams in Fig. 4.8 is shown to 
diverge again when the space dimension is given its physical value (d = 2). The theory 
can be rendered finite to this order by introducing the two-loop counter term 



= -16i nj—X{AX^ - a^) 



(4.66) 



Surprisingly, the counter term does not contain a simple pole in e. Since the beta 
function is determined by the residue of the simple pole, it follows that the one-loop 
result (4.6C) is unchanged at the two-loop level. 

The calculation of the beta function was extended to third order in the loop ex- 
pansion by Freedman, Lozano, and Rius [218] using differential regularization, and to 
fourth order by this author [219] applying the more conventional methods explained in 
this section. We found that the third and fourth order diagrams diverge like and 
1 / e^, respectively. In particular, the simple poles in e dropped out so that the one-loop 
beta function is unaffected by three- and four-loop corrections. It is tempting to con- 
jecture that this feature persists to all orders in perturbation theory, implying that — just 
as in the neutral system which corresponds to taking the limit 6 oo — the one-loop 
beta function ( 4.60 ) is exact. 

The beta function is schematically represented in Fig. 4.9. It yields an IR fixed 

2 

point A* = determined by 6* = a~^, or, equivalently, by the filling factor More 
precisely, the strength of the repulsive coupling at the fixed point A* = 7r(2Z + 1) in- 
creases with the number 2/ + 1 (= 1, 3,5,.. .) of flux quanta bound to the electron. 
The presence of the fixed point shows that the CSGL theory undergoes a 2nd-order 
phase transition when the chemical potential of the composite particles tends to zero. 
Since the self-energy E is identically zero, the nontrivial fixed point has nevertheless 
Gaussian exponents, v = ^, z = 2, 7y = 0. It should be noted that only the location of 
the fixed point depends on 9, the critical exponents — which in contrast to the strength 
of the coupling at the fixed point are independent of the regularization and renormal- 
ization scheme — are universal and independent of the filUng factor. This "triviality" is 
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/3(A) 




Figure 4.9: Schematic representation of the /J-function. 



in accord with the experimentally observed universality of the FQHE. A dependence 
of the critical exponents on could from the theoretical point of view hardly be made 
compatible with the hierarchy construction [220| which implies a cascade of phase 
transitions. From this viewpoint the present results are most satisfying: the CSGL the- 
ory is shown to encode a new type of field-induced 2nd-order quantum phase transition 
that is simple enough not to obscure the observed universality of the FQHE. 

Most of the properties found in the nonrelativistic theory are not shared by the 
relativistic version of the CSGL theory. On the contrary, by calculating the effec- 
tive potential, Pryadko and Zhang [215| argued that the relativistic theory undergoes a 
fluctuation-induced Ist-order transition — similar to the one in massless scalar electro- 
dynamics in 3h-1 spacetime dimensions first studied by Coleman and Weinberg [221 1. 
From the standpoint of dimensional analysis the discrepancy in critical behavior is not 
surprising. Because the effective spacetime dimensionality in the nonrelativistic theory 
is increased by one, the scaling dimension of a given operator generally differs in both 
theories. 

As a side remark we mention that the theory ( 4.39 ) with zero che mica l potential 
has also been successfully applied to the problem of any on scattering [ [216| ]. From it, 
the two-particle scattering amplitude, which was known from the work of Aharonov 
and Bohm [ 222 , 2 1 1 1 on the scattering of a charged particle from a magnetic flux tube, 
could be calculated. The result of Aharonov-Bohm could be reproduced provided the 
coupling constant was given the specific value \* — corresponding to the IR fixed 
point. It is sometimes presented as if taking this value is a mere choice. However, the 
scattering amplitudes considered here pertain to low-energy, small-momentum phe- 
nomena, and these are governed by the IR fixed point. In other words, in describ- 
ing these processes, there is no arbitrariness at all; one must fix the value of the self- 
coupling A to correspond to the IR fixed point. 

Because of the necessity to introduce a dimensional parameter to avoid infrared 
divergences and to be able to carr y out the renormalization-group program, t he sc ale 
invariance of the classical theory ( 4.39 ) with /i = 0, first discussed in Ref. [223|, is 
broken at nonfixed points. Only at the fixed point the theory is scale invariant. From 
the discussion above, it follows that the low-energy, small-momentum behavior of the 
system is also scale invariant since it is governed by the fixed point. 

The same remarks apply to the thermod ynam ic properties of an any on g as which 
can also be computed from the Lagrangian ( 4.39 ), as was done in Ref. [224|. These, 
too, are low-energy, small-momentum properties governed by the IR fixed point so that 
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one must set A* = when calculating them. 

The results obtained so far are gratifying. The CSGL theory, which captures the es- 
sential properties of a quantum Hall liquid, was shown to be capable of also describing 
the Hall-liquid-to-insulator transitions — at least in the absence of impurities. However, 
impurities plays an essential part in the FQHE, so a realistic description of the critical 
behavior of a quantum Hall liquid should account for impurities. Including (quenched) 
impurities in the superfluid state of repulsively interacting bosons was shown to lead to 
the phenomenon of localization which is of paramount importance to the FQHE. Since 
upon invoking a statistical gauge field this superfluid state represents a quantum Hall 
liquid, this is in principle a good starting point to study the quantum critical behavior of 
the CSGL model in the presence of impurities. There is, however, a technical problem. 
We have worked in a double epsilon expansion thereby leaving physical spacetime. 
Unfortunately, the Chern-Simons term governing the statistical gauge field of the the- 
ory is defined only in 2 + 1 dimensions. This makes a study of the critical behavior 
of the random CSGL theory in a double epsilon expansion impossible, and a different 
method, possibly nonperturbative in character, is required. 

Luckily, using general scaling arguments and information contained in the effective 
theories considered in this report one is able to acquire some additional understanding 
of the behavior of various quantum phase transitions studied in experiment, although 
precise estimates for critical exponents other than the dynamic one cannot be made on 
this basis. 



4.5 Scaling Theory 

The traditional scaling theory of classical 2nd-order phase transitions, first put for- 



ward by Widom [225], is easily extended to include the time dimension [I83| be 



cause relation ( 4-. 1 ) implies the presence of only one independent diverging scale. Let 
6 — K ~ Kc, with K the parameter that drives the phase transition, measure the dis- 
tance from the critical coupling Kc- A physical observable at the absolute zero of 
temperature 0{kQ, |k|, K) can in the critical region close to the transition be written as 

0{h,k,K)^^''oo{ko^t,\m, iT = 0), (4.67) 

where do is the dimension of the observable O. The right-hand side does not depend 
explicitly on K; only implicitly through ^ and The closer one approaches the critical 
coupling Kc, the larger the correlation length and time become. 

Since a physical system is always at some finite temperature, we have to investigate 



how the scaling law (4-. 67) changes when the temperature becomes nonzero. The easi- 
est way to include temperature in a quantum field theory is to go over to imaginary time 
T = it, with T restricted to the interval < t < /3. The temporal dimension becomes 
thus of finite extend. The critical behavior of a phase transition at finite temperature is, 
as we remarked before, still controlled by the quantum critical point provided £^t < (3. 
If this condition is fulfilled, the system does not s ee the finite extend of the time di- 
mension. Instead of the zero-temperature scaling ( 4.67| ), we now have the finite-size 
scahng 

0(fco,|k|,i^,T)=/3^°/^0(fco/3,|k|/3l/^/?/6), (T ^ 0). (4.68) 
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The distance to the quantum critical point is measured by the ratio ~ \5Y'^ /T. 
Let us continue to derive hyperscaling relations. To this end, we consider the two 



terms in the effective theory (1.177) quadratic in the Goldstone field Lp with m effec- 



tively set to 1 and write it in the most general form [E26| : 



C 



(2) 



(4.69) 



The coefficient ps is the superfluid mass density which in the presence of, for example, 
impurities does not equal mn — even at the absolute zero of temperature. The other 
coefficient, 

3fi 

— - lim noo(0,k), (4.70) 

op |kHO 



with IIoo the (0 0)-component of the full polarization tensor (1.179), involves the full 



compressibility and particle number density. This is because the chemical potential is 
represented in the effective theory by p = —do(p and 



dp^ 



(4.71) 



Equation (4.69) leads to the general expression of the sound velocity 



(4.72) 



at the absolute zero of temperature. 

Let us next assume that the chemical potential is the control parameter, so that 
6 (X p — Pc denotes the distance from the phase transition, and £^ Now, on the 

one hand, the singular part of the free energy density fsing arises from the low-energy, 
long-wavelength fluctuations of the Goldstone field. (Here, we adopted the common 
practice of using the symbol / for the density il/V and of referring to it as the free 
energy density.) The ensemble averages give 



-2z 



(4.73) 



On the other hand, dimensional analysis shows that the singular part of the free energy 
density scales near the transition as 



/si 



-(d+z) 



(4.74) 



Combining these hyperscaling arguments, we arrive at the following conclusions: 



-id~z) 



\6\ 



{d—z)iy 



(4.75) 



The first conclusion is consistent with the universal jump (3.12) predicted by Nelson 



and Kosterlitz [132] which corresponds to taking z — and d — 2. Since ^ ~ l^l""^, 
/sing can also be directly differentiated with respect to the chemical potential to yield 
for the the singular part of the compressibility 



{d+z)u-2 



(4.76) 
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Fisher and Fisher [ [226| continued to argue that there are two alternatives. Either k 

^sin^ 

. . . ^ . 1 1 

Z 



implying zv = 1; or the full compressibility k is constant, implying z — d. The 
former is consistent with the Gaussian values ly ~ ^, z — 2 found by Uzunov |189] 
for the pure case in d < 2. The latter is believed to apply to repulsively interacting 
bosons in a random media. These remarkable simple arguments thus predict the exact 
value z — d for the dynamic exponent in this case. 

The above hyperscaling arguments have been extended by Fisher, Grinstein, and 
Girvin [ I88| to inc lude the l/|x|-Coulomb potential. The quadratic terms in the effec- 
tive theory ( 1.201 ) may be cast in the general form 



C 



(2) 
off 



(4.77) 



where e is the renormalized charge. From ( 1.201 ) we find that to lowest order: 

The re normal ized charge is connected to the (0 0)-component of the full polarization 
tensor ( 1.107 ) via 

(4.79) 



lim 

IkHo noo(0,k) 



A simple hyperscaling argument like the ones given in the preceding paragraph shows 
that near the transition, the renormalized charge scales as 



(4.80) 



They then argued that in the presence of random impurities this charge is expected to 
be finite at the transition so that z = 1. This again is an exact results which replaces 
the value z = dof the neutral system. 

The quantum phase transitions we are considering take place in charged systems 
and are mainly probed by conductivity a or resistivity p measurements. To see how 
the conductivity a rel ates to the superfluid mass density ps, we minimally couple the 
effective theory ( 4.69 ) to an el ectro magnetic gauge field. The only relevant term for 
this purpose is the first one in ( 4.69 ) with Vip replaced by — eA, where we allow 
the superfluid mass density to vary in space and time. The term in the action quadratic 
in A then becomes in the Fourier representation 

Sa = -ie^ [ A(-fco, -k)p,(fco, k)A(fco, k). (4.81) 



The electromagnetic current. 



j(fco,k) 



6Sa 



(5A(-fco,-k) 
obtained from this action can be written as 

j(fco,k) = cr(fco,k)E(fco,k) 



(4.82) 



(4.83) 



(glAmstex 



4.6 Experiments 



135 



with the conductivity 

a{k) = -ie^^ (4.84) 
fco 

essentially given by the superfluid mass density divided by fcp, where it should be 
remembered that the mass m is effectively set to 1 here. By virtue of the scaling 



relation (4.75), it follows that a scales as [227] 



ar^C'"^-^^. (4.85) 

In other words, the scaling dimension of the conductivity and therefore that of the 
resistivity is zero in two space dimensions. 

Let us now consider a quantum phase transition triggered by changing the appUed 
magnetic field, i.e., 5 (x H — He- The critical field scales with ^ as i/c ^ ^o/^'^- In 
fact, this expresses a more fundamental result, namely that the scaling dimension 
of A is one, 

dA = 1, (4.86) 

so that |A| ~ ^"^ From this it in turn follows that E ^ and that 

the scaling dimension d^^, of Aq is z, 

dAo = z, (4.87) 

so that Aq ^ ^ For the DC conductivities in the presence of an external 

electric field we have on account of the general finite-size scaling form ( 4.68[ ) with 
fco = |k| = 0: 

a{H,T,E) = <;{d''yT,d'''-^+^'^/E). (4.88) 

This equation shows that conductivity measurements close to the quantum critical point 
collapse onto a single curve when plotted as function of the dimensionless combina- 
tions S'^^ /T and 5'^^^'^^^ / E. The best collapse of the data determines the values of uz 
and v{z + 1). In other words, the temperature and electric-field dependence determine 
the critical exponents v and z independently. 



4.6 Experiments 



4.6.1 Superconductor- To-Insulator Transition 

The first experiments we wish to discuss are those performed by Hebard and Paalanen 
on superconducting films in the presence of random impurities [192, 139|. It has been 
predicted by Fisher [228| that with increasing applied magnetic field such systems 
undergo a zero-temperature transition into an insulating state. (For a critical review of 



the experimental data available in 1993, see Ref. [229].) 



Let us restrict ourselves for the moment to the TA-plane of the phase diagram by 
setting the applied magnetic field H to zero. For given disorder strength A, the system 
then undergoes a Kosterlitz-Thouless transition induced by the unbinding of magnetic 
vortex pairs at a temperature Tkt well below the bulk transition temperature (see Sec. 
3.2). The Kosterlitz-Thouless temperature is gradually suppressed to zero when the 
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disorder strength approaches criticality A Ac. The transition temperature scales 
with the correlation length ^ | Ac - A|~'' as Tkt ~ 

In the i/A-plane, i.e., at T = 0, the situation is as follows. For given disorder 
strength, there is now at some critical value He of the applied magnetic field a phase 
transition from a superconducting state of pinned vortices and condensed Cooper pairs 
to an insulating state of pinned Cooper pairs and condensed vortices. The condensa- 
tion of vortices disorder the ordered state as happens in classical, finite temperature 
superfluid- and superconductor- to-normal phase transitions [|l7|]. When the disorder 
strength approaches criticality again, is gradually suppressed to zero. The critical 
field scales with ^ as i7c ^ ^o/^"^- Together, the scaling results for Tkt and He imply 
that 

^2/z 



He 



(4.89) 



This relation, linking the critical field to the Kosterlitz-Thouless temperature, provides 
a direct way to measure the dynamic exponent z at the _ff = 0, T = transition. 
This has been done first by Hebard and Paalanen [192, I39|. Their experimental de- 
termination of Tkt and H^ for five different films with varying amounts of impurities 
confirmed the relation (4.89) with 2/ z — 2.04 ± 0.09. The zero-temperature critical 
fields were obtained by plotting dpxx/^T\H versus H at the lowest accessible tem- 
perature and interpolating to the the field where the slope is zero. The resulting value 
z — 0.98 ± .04 is in accordance with Fisher's prediction [2281, 2 = 1, for a random 
system with a l/|x| -Coulomb potential. 

Hebard and Paalanen |192] also investigated the field-induced zero-temperature 
transition. The control parameter is here S oc H — He- When plotted as function of 
\H ~ Hel/T^^""''" they saw their resistivity data collapsing onto two branches; an 
upper branch tending to infinity for the insulating state, and a lower branch bending 
down for the superconducting state. The unknown product vhZh is experimentally 
determined by looking for which value the best scaling behavior is obtained. Further 



experiments carried out by Yazdani and Kapitulnik |23C] studying the electric-field 
dependence of the resistivity also determined the product vh{zh + !)■ The two inde- 
pendent measurements together fix the critical exponents vh and zh separately. From 
their best data, Yazdani and Kapitulnik extracted the values [ 23C ] 



ZH = 1.0±0.1, = 1.36 ±0.05. 



(4.90) 



4.6.2 Quantum-Hall Systems 

We continue to discuss the field-induced quantum phase transitions in quantum Hall 



systems. Since an excellent discussion recently appeared in the literature |182], we 
shall be brief, referring the reader to that review for a more thorough discussion and 
additional references. 

One can image transitions from one Hall liquid to another Hall liquid with a dif- 
ferent (integer or fractional) filling factor, or to the insulating state. Experiments seem 
to suggest that aU the quantum-Hall transitions are in the same universality class. The 
transitions are probed by measuring the conductivities Cxx and axy The scaling of the 



width of the transition regime with temperature as predicted by Eq. ( 4.8S| ) has been 
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corroborated by DC experiments on various transitions between integer quantum-Hall 
states which were all found to yield the value Xjvz — 0.42 ± 0.04 [231 1. Also the 
scaling b ~ has been corroborated by experiment which yielded the value 

j^(z + l) « 4.6 [232]. Together with the previous result obtained from the temperature 
scaling this gives 

zwl, 1/W2.3. (4.91) 

The value of the dynamic exponent strongly suggests that it is a result of the presence 
of the l/|x|-Coulomb potential. The correlation length exponent v is seen to be large. 



4.6.3 2d Electron Systems 

Recently, silicon MOSFET's at extremely low electron number densities has been stud- 



ied [202, 203, 204, 205). Earlier experiments at higher densities seemed to confirm the 



general believe, based on the work by Abrahams etal. [ 227 ] , that such two-dimensional 



electron systems do not undergo a quantum phase transition. In that influential paper, 
it was demonstrated that even weak disorder is sufficient to localize the electrons at the 
absolute zero of temperature thus excluding conducting behavior Electron-electron 
interactions were however not included. As we saw in Sec. 1.6, the l/|x|-Coulo mb in - 
teraction becomes important at low densities and the analysis of Abrahams et al. \ 227] 
no longer applies. 

The recent experiments have revealed a zero-temperature conductor-to-insulator 
transition triggered by a change in the charge carrier density n. That is, the distance to 
the critical point is in these systems measured by (5 c>c n — fic- Like in the quantum- 
Hall systems, these transitions are probed by measuring the resistivity. It scales with 
temperature near the transition according to the scaling form (4.88) with H set to zero. 
For n < fie, where the Coulomb interaction is dominant and fluctuations in the charge 
carrier density are suppressed, the electron system is insulating. On increasing the 
density, these fluctuations intensify and at the critical value fic, the system reverts to a 
conducting phase. By plotting their conductivity data as function of with vz = 

1.6 ± 0.1, Popovfc, Fowler, and Washburn [205] saw it collapse onto two branches; 
the upper branch for the conducting side of the transition, and the lower one for the 
insulating side. A similar collapse with a slightly different value \/vz — 0.83 ± 0.08 
was found in Ref. [203], where also the collapse of the data when plotted as function of 
gl^i/{z+i)u obtained. The best collapse resulted for l/(z + l)v = 0.37 ± 0.01, 
leading to 

z = 0.8 ±0.1, J/ =1.5 ±0.1. (4.92) 

The value for the dynamic exponent is close to the expected value z — 1 for a charged 
syste m with a l/|x|-Coulomb interaction, while that of v is surprisingly close to the 
value ( 4.90 ) found for the superconductor-to-insulator transition. 

A further experimental result for these two-dimensional electron systems worth 
mentioning is the suppression of the conductin g ph ase by an applied magnetic field 
found by Simonian, Kravchenko, and Sarachik [204]. They applied the field parallel 
to the plane of the electrons instead of perpendicular as is done in quantum-Hall mea- 
surements. In this way, the field presumably couples only to the spin of the electrons 
and the complications arising from orbital effects do not arise. At a fixed temperature. 
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a rapid initial raise in the resistivity was found with increasing field. Above a value of 
about 20 kOe, the resistivity saturates. It was pointed out that both the behavior in a 
magnetic filed, as well as in zero field strongly resembles that near the superconductor- 
to-insulator transition discussed above, suggesting that the conducting phase might in 
fact be superconducting. 

4.6.4 Conclusions 

We have seen that general scaling arguments combined with the effective theories, can 
be employed to understand the scaling behavior observed in various quantum phase 
transitions. Most of the experiments seem to confirm the expected value z = 1 for a 
random system with a l/|x|-Coulomb interaction. The number of different universality 
classes present is yet not known. Even if the conductor-to-insulator transition observed 
in silicon MOSFET's at low electron number densities turns out to be in the same 
universality class as the superconductor-to-insulator transition, there are still the field- 
induced transitions in quantum-Hall systems, which have a larger correlation-length 
exponent. 

The paradigm provided by a repulsively interacting Bose gas, seems to be a good 

starting point to describe the various systems. However, high-precision estimates cal- 
culated from this theory with impurities and a l/|x| -Coulomb interaction included are 
presently lacking. 
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